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Abstract 

We present the final part of the proof of the Kepler 
Conjecture. 



1 Overview 

This section describes the structure of the proof of 
the Kepler Conjecture. 

Theorem 1.1. (The Kepler Conjecture) No pack- 
ing of congruent balls in Euclidean three space has 
density greater than that of the face- centered cubic 
packing. 



This density is tt/vT8 ~ 0.74. 
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Figure 1: The face-centered cubic packing 

The proof of this result is scattered through- 
out several papers. (Every article in the bibliog- 
raphy will be required.) Here, we describe the top- 
level outline of the proof and give references to the 
sources of the details of the proof. 

By a packing, we mean an arrangement of balls 
that are non-overlapping in the sense that the inte- 
riors of the balls are pairwise disjoint. Consider a 
packing of congruent balls in Euclidean three space. 
There is no harm in assuming that all the balls 
have unit radius. The density of a packing does 
not decrease when balls are added to the packing. 



Thus, to answer a question about the greatest possi- 
ble density we may add non-overlapping balls until 
there is no room to add further balls. Such a pack- 
ing will be said to be saturated. 

Let A be the set of centers of the balls in a satu- 
rated packing. Our choice of radius for the balls im- 
plies that any two points in A have distance at least 
2 from each other. We call the points of A vertices. 
Let B(x, r) denote the ball in Euclidean three space 
at center x and radius r. Let Six, r, A) be the finite 
density, defined by the ratio of A(x, r, A) to the vol- 
ume of B(x,r), where A(x, r, A) is defined as the 
volume of the intersection with B(x, r) of the union 
of all balls in the packing. Set A(x, r) = A<~)B(x, r). 

Recall that the Voronoi cell Cl(v) around a vertex 
ii £ A is the set of points closer to v than to any 
other ball center. The volume of each Voronoi cell in 
the face-centered cubic packing is \/32. This is also 
the volume of each Voronoi cell in the hexagonal- 
close packing. 

Let a : A — > R be a function. We say that a is 
negligible if there is a constant C\ such that for all 
r > 1, we have 

53 < Cir 2 . 

We say that the function a is fee- compatible if for 
all v € A we have the inequality 



32 < vol(n(u)) + a(v). 

Lemma 1.2. If there exists a negligible fec- 
compatible function a : A — > K for a saturated pack- 
ing A, then there exists a constant C such that for 
all r > 1, we have 



8(x,r,A) < 7r/Vl8 + C/r. 

Proof. The numerator A(x, r, A) of S(x, r, A) is at 
most the product of the volume of a ball 4ir/3 with 
the number |A(ar, r+ 1) | of balls intersecting B(x, r). 
Hence 



A(x,r,A) < \A(x,r + 1)|4tt/3. 



(1) 
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In a saturated packing each Voronoi cell is con- 
tained in a ball of radius 2 centered at the center 
of the cell. The volume of the ball B(x, r + 3) is at 
least the combined volume of Voronoi cells lying en- 
tirely in the ball. This observation, combined with 
fcc-compatibility and negligibility, gives 



The density of the face-centered cubic packing is 
a weighted average of these two ratios 



_ <het_ 
18 " 3 
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32|A(a;,r + 1)| < ^ (a(v) + vol(Q(V))) 

veA(x,r+l) 

< Ci(r + l) 2 + vol 5(2;, r + 3) 

< Ci(r + l) 2 + (1 + 3/r) 3 vol5(x, r) 

(2) 

Divide through by vol B(x,r) and eliminate 
\A(x,r + 1)| between Inequality (1) and Inequality 
(2) to get 

6(x,r,A) < -Wl + 3/r) 3 + Ci- 1 
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The result follows for an appropriately chosen con- 
stant C. □ 

Remark 1.3. We take the precise meaning of the 
Kepler Conjecture to be a bound on the essential 
supremum of the function S(x,r) as r tends to in- 
finity Lemma 1.2 implies that the essential supre- 
mum of S(x, r, A) is bounded above by 7r/\/l8, pro- 
vided a negligible fcc-compatible function can be 
found. The strategy will be to define a negligible 
function, and then to solve an optimization prob- 
lem in finitely many variables to establish that it is 
fcc-compatible. 

The paper [2] defines a compact topological space 
X and a continuous function a on that space. 

The topological space X is directly related to 
packings. If A is a saturated packing, then there is 
a geometric object D(v, A) constructed around each 
vertex v G A. D(v, A) depends on A only through 
the vertices in A at distance at most 4 from v. The 
objects D(v,A) are called decomposition stars, and 
the space of all decomposition stars is precisely X. 

The following constants arise in that paper. 

Let Stet be the packing density of a regular tetra- 
hedron. That is, let 5* be a regular tetrahedron of 
edge length 2. Let B the part of S that lies within 
distance 1 of some vertex. Then Stet is the ratio 
of the volume of B to the volume of S. We have 
Stet — V& arctan(v / 2/5). 

Let 5 oc t be the packing density of a regular octa- 
hedron of edge length 2, again constructed as the 
ratio of the volume of points within distance 1 of a 
vertex to the volume of the octahedron. 



This determines the exact value of S oct in terms of 
Stet- We have S oct w 0.72. 

Let pt = -tt/3 + V2S tet w 0.05537 

The following conjecture is made in [2] 

Conjecture 1.4. The maximum of a on X is the 

constant 8 pt w 0.442989. 

Lemma 1.5. An affirmative answer to Conjec- 
ture 1.4 implies the existence of a negligible fcc- 
compatible function for every saturated packing A. 

Proof. For any saturated packing A define a func- 
tion a : A — > R by 

-cj(D(v, A))/(45 oct )+47r/(3£ oct ) = vol(fi(«))+o(«). 

Negligibility follows from [2, Prop. 3.14 (proof)]. 
The upper bound of 8 pt gives a lower bound 

-8pt/(45 oct ) + 47r/(3<5 oct ) < vol(fi(«)) + a(v). 

The left-hand side of this inequality evaluates to 
V32, and this establishes fcc-compatibility. □ 

Theorem 1.6. Conjecture 1.4 is true. That is, the 
maximum of the function a on the topological space 
X of all decomposition stars is 8 pt. 

Theorem 1.6, Lemma 1.5, and Lemma 1.2 com- 
bine to give a proof of the Kepler Conjecture 1.1 

Let to = 1.255 (2to = 2.51). This is a parameter 
that is used for truncation throughout the series 
of papers on the Kepler Conjecture. The signifi- 
cance of this particular choice of truncation param- 
eter will be explained elsewhere. 

Let U(v, A) be the set of vertices in A at distance 
at most 2to from v. From a decomposition star 
D(v, A) it is possible to recover U(v, A) (at least up 
to Euclidean translation: U 1— » U + y, for y € K 3 ). 
We can completely characterize the decomposition 
stars at which the maximum of a is attained. 

Theorem 1.7. Let D be a decomposition star at 
which the maximum 8pt is attained. Then the set 
U(D) of vectors at distance at most 2t from the 
center has cardinality 12. Up to Euclidean mo- 
tion, U(D) is the kissing arrangement of the 12 balls 
around a central ball in the face-centered cubic pack- 
ing or hexagonal-close packing. 
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1.1 Outline of proofs 

To prove Theorems 1.6 and 1.7, we wish to show 
that there is no counterexample. That is, we wish 
to show that there is no decomposition star D with 
value u(D) > 8 pt. We reason by contradiction, as- 
suming the existence of such a decomposition star. 
With this in mind, we call D a contravening decom- 
position star, if 

a(D) >8pt. 

In much of what follows we will tacitly assume that 
every decomposition star under discussion is a con- 
travening one. Thus, when we say that no decom- 
position stars exist with a given property, it should 
be interpreted as saying that no such contravening 
decomposition stars exist. 

To each contravening decomposition star, we as- 
sociate a (combinatorial) plane graph. A restrictive 
list of properties of plane graphs is described in Sec- 
tion 2.3. Any plane graph satisfying these proper- 
ties is said to be tame. All tame plane graphs have 
been classified. (There are several thousand, up to 
isomorphism.) Theorem 2.1, asserts that the plane 
graph attached to each contravening decomposition 
star is tame. By the classification of such graphs, 
this reduces the proof of the Kepler Conjecture to 
the analysis of the decomposition stars attached to 
the finite explicit list of tame plane graphs. 

A few of the tame plane graphs are of partic- 
ular interest. Every decomposition star attached 
to the face-centered cubic packing gives the same 
plane graph (up to isomorphism). Call it G/ cc . 
Likewise, every decomposition star attached to the 
hexagonal-close packing gives the same plane graph 
Ghcp- Let X cr u be the set of decomposition stars 
D such that the set U(D) of vertices is the kissing 
arrangement of the 12 balls around a central ball in 
the face-centered cubic or hexagonal-close packing. 
There are only finitely many orbits of X cr it under 
the group of Euclidean motions. 




Figure 2: The plane graphs Gf cc and Ghcp 



In [2, Lemma 3.13], the necessary local analysis is 
carried out to prove the following local optimality. 

Lemma 1.8. A decomposition star whose plane 
graph is Gf cc or Ghcp has score at most 8pt, with 
equality precisely when the decomposition star be- 
longs to Xcrit- 

In light of this result, we prove 1.6 and 1.7 by 
proving that any decomposition star whose graph 
is tame and not equal to G/ cc or Ghcp is not con- 
travening 

There is one more tame plane graph that is par- 
ticularly troublesome. It is the graph G pen t ob- 
tained from the pictured configuration of twelve 
balls tangent to a given central ball (Figure 3). 
(Place a ball at the north pole, another at the south 
pole, and then form two pentagonal rings of five 
balls.) This case requires individualized attention. 
S. Ferguson proves in [7] that if D is any decompo- 
sition star with this graph, then o~(D) < 8 pt. 




Figure 3: The plane graph G pen t 



To eliminate the remaining cases, more-or-less 
generic arguments can be used. A linear program is 
attached to each tame graph G. The linear program 
can be viewed as a linear relaxation of the nonlin- 
ear optimization problem of maximizing a over all 
decomposition stars with a given tame graph G. 
Because it is obtained by relaxing the constraints 
on the nonlinear problem, the maximum of the lin- 
ear problem is an upper bound on the maximum 
of the original nonlinear problem. Whenever the 
linear programming maximum is less than 8 pt, it 
can be concluded that there is no contravening de- 
composition star with the given tame graph G. This 
linear programming approach eliminates most tame 
graphs. 

When a single linear program fails to give the 
desired bound, it is broken into a series of linear 
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programming bounds, by branch and bound tech- 
niques. For every tame plane graph G other than 
Ghc P , Gf cc , and G pent , we produce a series of linear 
programs that establish that there is no contraven- 
ing decomposition star with graph G. When every 
face of the plane graph is a triangle or quadrilat- 
eral, this is accomplished in [5]. The general case is 
completed in the final sections of this paper. 

1.2 Organization of the Paper 

This paper has been written in such a way that 
several sections can be skipped without disrupting 
the expository flow. To get a general overview of 
the paper, Sections 1, 2, 4, and 10 may be read. To 
get a somewhat more detailed view of the results 
and methods of proof, Sections 1, 2, 3, 4, 5, 10, 11, 
and the first two appendices may be read. 

The paper has been written so as to be substan- 
tially independent of the other papers in the series. 
An appendix gives a summary of the results of [6] in 
a form that is suited to the purposes of this paper. 

A historical account of the Kepler Conjecture can 
be found at [1]. An introduction to the ideas of the 
proof can be found in [11]. An introduction to the 
algorithms can be found at [9]. Speculation on a 
second-generation design of a proof can be found in 
[9] and [12]. 

2 Tame Plane Graphs 

We give a succinct statement of properties that hold 
of all plane graphs considered in the proof of the 
Kepler conjecture. Such graphs are said to be tame. 
A list of all tame graphs has been generated by 
computer. 

2.1 Basic Definitions 

An n-cycle is a finite set C of cardinality n, together 
with a cyclic permutation s of C. We write s in the 
form v i — ► s(v,C), for v G C. The element s(v,C) 
is called the successor of v (in C). A cycle is a 
n-cycle for some natural number n. By abuse of 
language, we often identify C with the cycle. The 
natural number n is the length of the cycle. 

(We may assume that all vertices of all graphs lie 
in some large finite set of vertices f2, if we wish to 
arrange that only finitely many graphs occur in this 
discussion.) 



Let G be a nonempty finite set of cycles (called 
faces) of length at least 3. The elements of faces 
are called the vertices of G. An unordered pair of 
vertices {v, w} such that one element is the succes- 
sor of the other in some face is called an edge. The 
vertices v and w are then said to be adjacent. The 
set G is a plane graph if three conditions hold. 

1. If an element v has successor w in some face 
F, then there is a unique face (call it s'(F, v)) 
in G for which v is the successor of w. (Thus, 
v = s(w,s'(F,v)), and each edge occurs twice 
with opposite orientation.) 

2. For each vertex v, the function F i— > s'(F,v) 
is a cyclic permutation of the set of faces con- 
taining v. 

3. Euler's formula holds relating the number of 
vertices V, the number of edges E, and the 
number of faces F: 

V -E + F = 2. 

(The set of vertices and edges of a plane graph form 
a planar graph in the usual graph-theoretic sense of 
admitting an embedding into the plane.) 

Let len be the length function on faces. Let tri(v) 
be the number of triangles containing a vertex v. 
Faces of length 3 are called triangles, those of length 
4 are called quadrilaterals, and so forth. An face of 
length at least 5 is called an exceptional face. 

Two plane graphs are properly isomorphic if there 
is a bijection of vertices inducing a bijection of faces. 
For each plane graph, there is an opposite plane 
graph G op obtained by reversing the cyclic order of 
vertices in each face. A plane graph G is isomorphic 
to another if G or G op is properly isomorphic to the 
other. 

The degree of a vertex is the number of faces it 
belongs to. An n-circuit in G is a cycle C in the 
vertex-set of G, such that for every v € C, it forms 
an edge in G with its successor: that is, (v, s(v, C)) 
is an edge of G. In a plane graph G we have a 
combinatorial form of the Jordan curve theorem: 
each n-circuit determines a partition of G into two 
sets of faces. 

The type of a vertex is defined to be a triple of 
non-negative integers (p,q,r), where p is the num- 
ber of triangles containing the vertex, q is the num- 
ber of quadrilaterals containing it, and r is the num- 
ber of exceptional faces. When r = 0, we abbreviate 
the type to the ordered pair (p,q). 
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2.2 Weight Assignments 

We call the constant 14.8, which arises repeatedly 
in this section, the target. 

Define a : N -> R by 



r 14.8 n 
1.4 
1.5 




0,1,2, 
n = 3, 
n = 4, 
otherwise. 



Define b : N x N -> M by 6(p, g) = 14.8, except for 
the values in the following table (with the under- 
standing that x = 14.8): 



p = 
1 
2 
3 
4 
5 
G 



9 = 

x 
x 

X 
X 

4.139 
0.55 
6.339 



1 

x 
x 

8.5 
3.642 
3.781 
11.22 

x 



2 

x 
6.95 
4.756 
8.334 

x 

x 

x 



3 

7.135 
7.135 
12.981 

x 

X 
X 
X 



4 

10.649 

x 
x 

X 
X 

x 
x 



Define c : N 



by 



Define d : N 



d= I 



1 



-1.03 
-2.06 
-3.03 

by 
fo 

2.378 
4.896 
7.414 
9.932 
10.916 
U4.8 



n= 3, 
n = 4, 
n = 5, 
n = 6, 
otherwise. 



n = 3, 
n = 4, 
n = 5, 
n = 6, 
n = 7, 
n = 8, 
otherwise. 



A set V of vertices is called a separated set of 
vertices if the following four conditions hold. 

1. For every vertex in V there is an exceptional 
face containing it. 

2. No two vertices in V are adjacent. 

3. No two vertices in V lie on a common quadri- 
lateral. 

4. Each vertex in V has degree 5. 



A weight assignment of a plane graph G is a 
function w : G — > R taking values in the set of 
non-negative real numbers. A weight assignment is 
admissible if the following properties hold: 

1. If the face F has length n, then w(F) > d{n) 

2. If v has type (p, q), then 

J2w(F)>b(p,q). 

veF 

3. Let V be any set of vertices of type (5,0). If 
the cardinality of V is k < 4, then 

w(F) > 0.55k. 

4. Let V be any separated set of vertices. Then 

J2 (w(F)-d(len(F)))>J2a(tri(v)). 
vnF^ii v<=v 

The sum J2 F w(F) is called the total weight of w. 

2.3 Plane Graph Properties 

We say that a plane graph is tame if it satisfies the 
following conditions. 

1. The length of each face is (at least 3 and) at 
most 8. 

2. Every 3-circuit is a face or the opposite of a 
face. 

3. Every 4-circuit surrounds one of the cases il- 
lustrated in Figure 4. 

4. The degree of every vertex is (at least 2 and) 
at most 6. 

5. If a vertex is contained in an exceptional face, 
then the degree of the vertex is at most 5. 



6. 



J2c(len(F))>8, 



7. There exists an admissible weight assignment 
of total weight less than the target, 14.8. 

It follows from the definitions that the abstract 
vertex-edge graph of G has no loops or multiple 
joins. Also, by construction, every vertex lies in at 
least two faces. Property 6 asserts that the graph 
has at least eight triangles. 
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Figure 4: Tame 4-circuits 



2.4 Classification of Tame Plane 
Graphs 

A list of several thousand plane graphs appears at 
[8]. 

Theorem 2.1. Every tame plane graph is isomor- 
phic to a plane graph in this list. 



3 Proof of Classification 

The results of this section are not needed except in 
the proof of Theorem 2.1. 

Computers are used to generate a list of all tame 
plane graphs and to check them against the archive 
of tame plane graphs. We will describe a finite state 
machine that produces all tame plane graphs. This 
machine is not particularly efficient, and so we also 
include a description of pruning strategies that pre- 
vent a combinatorial explosion of possibilities. 

3.1 Basic Definitions 

A partial plane graph is a plane graph with addi- 
tional data: every face is marked as "complete" or 
"incomplete." We call a face complete or incom- 
plete according to the markings. We require the 
following condition. 

• No two incomplete faces share an edge. 



Each unmarked plane graph is identified with the 
marked plane graph in which every face is complete. 
We represent a partial plane graph graphically by 
deleting one face (the face at infinity) and drawing 
the others and shading those that are complete. 

A patch is a partial plane graph P with two dis- 
tinguished faces F\ and P 2 , such that the following 
hold. 

• Every vertex of P lies in Fi or F2 . 

• The face F2 is the only complete face. 

• Fi and F2 share an edge. 

• Every vertex of F2 that is not in Fi has degree 
2. 

Fi and F2 will be referred to as the distinguished 
incomplete and the distinguished complete faces, re- 
spectively. 

Patches can be used to modify a partial plane 
graph as follows. Let F be an incomplete face of 
length n in a partial plane graph G. Let P be a 
patch whose incomplete distinguished face F\ has 
length n. Replace P with a properly isomorphic 
patch P' in which the image of Pi is equal to F op 
and in which no other vertex of P' is a vertex of G. 
Then 

G' = {F 1 eGUP':F'^ F° p , F' ^ F} 

is a partial plane graph. Intuitively, we cut away 
the faces F and Pi from their plane graphs, and 
glue the holes together along the boundary (Figure 
5). (It is immediate that the Condition 3.1 in the 
definition of partial plane graphs is maintained by 
this process.) There are n distinct proper ways of 
identifying F\ with F op in this construction, and 
we let (fi be this identification. The isomorphism 
class of G' is uniquely determined by the isomor- 
phism class of G, the isomorphism class of P, and 
(fi (ranging over proper bijections (fi : F\ 1— > F op ). 

3.2 A Finite State Machine 

For a fixed N we define a finite state machine as 
follows. The states of the finite state machine are 
isomorphism classes of partial plane graphs G with 
at most N vertices. The transitions from one state 
G to another are isomorphism classes of pairs (P, (fi) 
where P is a patch, and (fi pairs an incomplete face 
of G with the distinguished incomplete face of P. 
However, we exclude a transition (P, (fi) at a state 
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Figure 5: Patching a plane graph 



if the resulting partial plane graphs contains more 
than N vertices. Figure 5 shows two states and a 
transition between them. 

The initial states /„ of the finite state machine 
are defined to be the isomorphism classes of partial 
plane graphs with two faces: 

{(1,2,... ,n),(n,n-l,... ,1)} 

where n < N, one face is complete, and the other is 
incomplete. In other words, they are patches with 
exactly two faces. 

A terminal state of this finite state machine is one 
in which every face is complete. By construction, 
these are (isomorphism classes of) plane graphs 
with at most N vertices. 

Lemma 3.1. Let G be a plane graph with at most 
N vertices. Then its state in the machine is reach- 
able from an initial state through a series of transi- 
tions. 

Proof. Pick an face in G of length n and identify it 
with the complete face in the initial state /„. At 
any stage at state G", we have an identification of 
all of the vertices of the plane graph G' with some 
of the vertices of G, and an identification of all of 
the complete faces of G' with some of the faces of 
G (all faces of G are complete). Pick an incomplete 
face F of G' and an oriented edge along that face. 
We let F' be the complete face of G with that edge, 
with the same orientation on that edge as F. Cre- 
ate a patch with distinguished faces F\ — F op and 
i<2 = F' . (Fi and F2 determine the patch up to 
isomorphism.) It is immediate that the conditions 
defining a patch arc fulfilled. Continue in this way 
until a graph isomorphic to G is reached. □ 

Remark 3.2. It is an elementary matter to gener- 
ate all patches P such that the distinguished faces 



have given lengths n and m. Patching is also en- 
tirely algorithmic, and thus by following all paths 
through the finite state machine, we obtain all plane 
graphs with at most N vertices. 

3.3 Pruning Strategies 

Although we reach all graphs in this manner, it is 
not computationally efficient. We introduce prun- 
ing strategies to increase the efficiency of the search. 
We can terminate our search along a path through 
the finite state machine, if we can determine: 

1. Every terminal graph along that path violates 
one of the defining properties of tameness, or 

2. An isomorphic terminal graph will be reached 
by some other path that will not be terminated 
early. 

Here are some pruning strategies of the first type 
(1). They are immediate consequences of the con- 
ditions of the defining properties of tameness. 

• If the current state contains an incomplete face 
of length 3, then eliminate all transitions, ex- 
cept for the transition that carries the partial 
plane graph to a partial plane graph that is the 
same in all respects, except that the face has 
become complete. 

• If the current state contains an incomplete face 
of length 4, then eliminate all transitions ex- 
cept those that lead to the possibilities of Sec- 
tion 2.3, Property 3, where in Property 3 each 
depicted face is interpreted as being complete. 

• Remove all transitions with patches whose 
complete face has length greater than 8. 

• It is frequently possible to conclude from the 
examination of a partial plane graph that no 
matter what the terminal position, any admis- 
sible weight assignment will give total weight 
greater than the target (14.8). In such cases, 
the all transitions out of the partial plane graph 
can be pruned. 

To take a simple example of the last item, we 
observe that weights are always non-negative, and 
that the weight of a complete face of length n is 
at least d(n). Thus, if there are complete faces 
F± , . . . , Fk of lengths n-± , . . . , rife, then any admis- 
sible weight assignment has total weight at least 
y^_] d(n,i). If this number is at least the target, 
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then no transitions out of that state need be con- 
sidered. 

More generally, we can apply all of the inequal- 
ities in the definition of admissible weight assign- 
ment to the complete portion of the partial plane 
graph to obtain lower bounds. However, we must be 
careful, in applying Property 4 of admissible weight 
assignments, because vertices that are not adjacent 
at an intermediate state may become adjacent in 
the complete graph. Also, vertices that do not lie 
together in a quadrilateral at an intermediate state 
may do so in the complete graph. 

Here are some pruning strategies of the second 
type (2). 

• At a given state it is enough to fix one incom- 
plete face and one edge of that face and then 
to follow only the transitions that patch along 
that face and add a complete face along that 
edge. (This is seen from the proof of Lemma 
3.1.) 

• In leading out from the initial state /„, it is 
enough to follow paths in which every added 
complete face has length at most n. (A graph 
with a face of length to, for to > n, will be also 
be found downstream from I m .) 

• Make a list of all type (p, q) with b(p, q) < 14.8. 
Remove the initial states 1% and 1^ and create 
new initial states 1 PA {I p . q , I p „, etc.) in the 
finite state machine. Define the state I p . q to 
be one consisting of p + q + 1 faces, with p com- 
plete triangles and q complete quadrilaterals 
all meeting at a vertex (and one other incom- 
plete face away from v). (If there is more than 
one way to arrange p triangles and q quadri- 
laterals, create states I PA , I' , I p ' q , for each 
possibility. See Figure 6.) Put a linear order 
on states I p . q . In state transitions downstream 
from J Pj9 disallow any transition that creates a 
vertex of type (p 1 , q'), for any (p\ q') preceding 
{p, q) in the imposed linear order. 




Figure 6: States and i 3 



This last pruning strategy is justified by the fol- 
lowing lemma, which classifies vertices of type (p, q) . 

Lemma 3.3. Let A and B be triangular or quadri- 
lateral faces that have at least 2 vertices in common 
in a tame graph. Then the faces have exactly two 
vertices in common, and an edge is shared by the 
two faces. 

Proof. Exercise. Some of the configurations that 
must be ruled out are shown in Figure 7. Some 
properties that are particularly useful for the exer- 
cise are Properties 2 and 3 of tameness, and Prop- 
erty 2 of admissibility. □ 




Figure 7: Some impossibilities 

Once a terminal position is reached it is checked 
to see whether it satisfies all the properties of tame- 
ness. 

Duplication is removed among isomorphic termi- 
nal plane graphs. It is not an entirely trivial proce- 
dure for the computer to determine whether there 
exists an isomorphism between two plane graphs. 
This is accomplished by computing a numerical in- 
variant of a vertex that depends only on the local 
structure of the vertex. If two plane graphs are 
properly isomorphic then the numerical invariant 
is the same at vertices that correspond under the 
proper isomorphism. If two graphs have the number 
of vertices with the same numerical invariants, they 
become candidates for an isomorphism. All pos- 
sible numerical-invariant preserving bijections are 
attempted until an proper isomorphism is found, 
or until it is found that none exist. If there is no 
proper isomorphism, the same procedure is applied 
to the opposite plane graph to find any possible 
orientation-reversing isomorphism. 

This same isomorphism-producing algorithm is 
used to match each terminal graph with a graph 
in the archive. It is found that each terminal graph 
matches with one in the archive. (The archive was 



originally obtained by running the finite state ma- 
chine and making a list of all the terminal states up 
to isomorphism that satisfy the given conditions.) 

In this way Theorem 2.1 is proved. 

4 From Stars to Graphs 
4.1 Plane Graphs 

A plane graph G is attached to every contravening 
decomposition star as follows. For simplicity, take 
the decomposition star to be attached to a ball of 
the packing, centered at the origin. From the de- 
composition star D, it is possible to determine the 
coordinates of the set U(D) of vertices at distance 
at most 2to from the origin. 

If we draw a geodesic arc on the unit sphere at 
the origin with endpoints at the projections of v\ 
and V2 for every pair of vertices vi, Ua G U(D) such 
that | Vi |, \v2 1, |«i — i>2| < 2to, we obtain a plane 
graph that breaks the unit sphere into regions called 
standard regions. (The arcs do not meet except at 
endpoints [3, Lemma 3.10].) Each standard region 
is defined as the closure in the unit sphere of a con- 
nected component of the unit sphere with all arcs 
removed. 

For a given standard region, we consider the arcs 
forming its boundary together with the arcs that 
are internal to the standard region. We consider the 
points on the unit sphere formed by the endpoints 
of the arcs, together with the projections to the unit 
sphere of vertices in U whose projection lies in the 
interior of the region. 

Remark 4.1. The system of arcs and vertices as- 
sociated with a standard region in a contravening 
example must be a polygon, or one of the following 
configurations of Figure 8 (see [6, Corollary 4.4]). 




Figure 8: Non-polygonal standard regions 



Remark 4.2. Observe that one case is bounded 
by a triangle and a pentagon, and that the others 
are bounded by a polygon. Replacing the triangle- 
pentagon arrangement with the bounding pentagon 
and replacing the others with the bounding poly- 
gon, we obtain a partition of the sphere into simple 
polygons. Each of these polygons is a single stan- 
dard region, except in the triangle-pentagon case 
(Figure 9), which is a union of two standard regions 
(a triangle and a eight-sided region). 




Figure 9: An aggregate forming a pentagon 

Remark 4.3. To simplify further, if we have an ar- 
rangement of six standard regions around a vertex 
formed from five triangles and one pentagon, we re- 
place it with the bounding octagon (or hexagon). 
See Figure 10. (It will be shown in Lemma 6.3 and 
Section 9.2 that there is at most one such configura- 
tion in the standard decomposition of a contraven- 
ing decomposition star, so we will not worry here 
about how to treat the case of two overlapping con- 
figurations of this sort.) 




Figure 10: Degree 6 aggregates 

In summary, we have a plane graph that is ap- 
proximately that given by the standard regions of 
the decomposition star, but simplified to a bound- 
ing polygon when one of the configurations of Re- 
marks 4.2 and 4.3 occur. We refer to the combi- 
nation of standard regions into a single face of the 
graph as aggregation. We call it the plane graph 
G = G(D) attached to a contravening decomposi- 
tion star D. 
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Proposition 5.1 will show the vertex set U is non- 
empty and that the graph G(D) is non-empty. 

When we refer to the plane graph in this manner, 
we mean the combinatorial plane graph in the sense 
of Section 2, as opposed to the embedded metric 
graph on the unit sphere formed from the system 
of geodesic arcs. Given a vertex v in G(D), there is 
a uniquely determined vertex v(D) of U(D) whose 
radial projection to the unit sphere determines v. 
We call v(D) the corner in U(D) over v. 

By definition, the plane graphs in Section 2 do 
not have loops or multiple joins. The edges of G(D) 
are defined by triangles whose sides vary between 
lengths 2 and 2i - The angles of such a triangle 
are strictly less than it. This implies that the edges 
of the metric graph on the unit sphere always have 
arc-length strictly less than tt. In particular, the 
endpoints are never antipodal. A loop on the com- 
binatorial graph corresponds to a edge on the metric 
graph that is a closed geodesic. A multiple join on 
the combinatorial graph corresponds on the metric 
graph to a pair of points joined by multiple mini- 
mal geodesies, that is, a pair of antipodal points on 
the sphere. By the arc-length constraints on edges 
in the metric graph, there are no loops or multiple 
joins in the combinatorial graph G(D). 

In Definition 2.3, a plane graph satisfying a cer- 
tain restrictive set of properties is said to be tame. 
If a plane graph G(D) is associated with a contra- 
vening decomposition star D, we call G(D) a con- 
travening plane graph. 

Theorem 4.4. Every contravening plane graph is 
tame. 

This theorem is one of the main results of this 
paper. Its proof occupies Sections 5 through 9. 

In Theorem 2.1, the tame graphs arc classified up 
to isomorphism. As a corollary, we have an explicit 
list containing all contravening plane graphs. 

5 Contravention is tame 

This section begins the proof of Theorem 4.4. 

To prove Theorem 4.4, it is enough to show that 
each defining property of tameness is satisfied. This 
is the substance of results in the next sections. 
First, we prove the promised non-degeneracy result. 

Proposition 5.1. The construction of Section 4-1 
associates a (nonempty) plane graph with at least 
two faces to every decomposition star D with 
a{D) > 0. 



Proof. First we show that decomposition stars with 
a(D) > have non-empty vertex sets U. (Recall 
that U is the set of vertices of distance at most 2£ 
from the center) . The vertices of U are used in [2] to 
create all of the structural features of the decompo- 
sition star: quasi-regular tetrahedra, quarters, and 
so forth. If U is empty, the decomposition star is a 
solid containing the ball B(to) of radius to, and its 
score is just 

a(D) = vor(D) 

= -AS oct vol{D) + Att/3 
<-4S oct vol{B{t ))+ATr/3 < 0. 

By hypothesis, cr(D) > 0. So U is not empty. 

The results of [2, Eqn 3.11] show that the function 
a can be expressed as a sum of terms indexed by 
the standard regions R. It is proved in [2, Lemma 
3.13 ] that an < 0, unless R is a triangle. Thus, a 
decomposition star with positive score must have at 
least one triangle. Its complement contains a sec- 
ond standard region. Even after we form aggregates 
of distinct standard regions to form the simplified 
plane graph (Remarks 4.2 and 4.3), there certainly 
remain at least two faces. □ 

Proposition 5.2. The plane graph of a contraven- 
ing decomposition star satisfies Property 1 of tame- 
ness: The length of each face is (at least 3 and) at 
most 8. 

Proof. By the construction of the graph, each face 
has at least three edges. The upper bound of 8 faces 
is [6, Corollary 4.4]. □ 

Proposition 5.3. The plane graph of a contraven- 
ing decomposition star satisfies Property 2 of tame- 
ness: Every 3- circuit is a face or the opposite of a 
face. 

Proof. The simplifications of the plane graph in 
Remarks 4.2 and 4.3 do not produce any new 3- 
circuits. (See the accompanying figures.) The result 
is proved in [3, Lemma 3.7]. The proof is purely ge- 
ometric: it expresses a geometric impossibility true 
of all decomposition stars. □ 

Proposition 5.4. Contravening graphs satisfy 
Property 4 of tameness: The degree of every vertex 
is at least 2 and at most 6. 

Proof. The statement that degrees are at least 2 
trivially follows because each vertex lies on at least 
one polygon, with two edges at that vertex. 

The impossibility of a vertex of degree 7 or more 
is found in the proof of [5, Lemma 6.2]. □ 
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6 More Tame Properties 

This section continues in the proof that all contra- 
vening plane graphs are tame. 

6.1 Linear Programs 

To continue with the proof that contravening plane 
graphs are tame, we need to introduce some more 
notation and methods. 

The function a can be written is a sum of func- 
tions on the space of decomposition stars, indexed 
by standard regions. (Sec [2, Eqn 3.11].) Write 
or for the term indexed by the standard region R. 
Then 

*(£>) = 5>fl(r>), 

R 

for every decomposition star. 

Let C = l/(2arctan(V2/5). Let sol(i2) denote 
the solid angle of a region R. We write tr for the 
following modification of or: 

T R {D)=sol{R)Cpt-a R {D) 

and 

t(D) = tr(D) = 4n(pt - a{D). 

Since 4n(pt is a constant, r and a contain the same 
information, but r is often more convenient to work 
with. A contravening decomposition star satisfies 

t(D) < 4ir(pt -8pt= (4ttC - 8)pt. (3) 

The constant (47r£ — 8) pt (and its upper bound 
14.8 pi) will occur repeatedly in the discussion that 
follows. 

Let TZ(F) be the set of standard regions associ- 
ated with a face F in the plane graph. We occasion- 
ally combine the terms a ft associated with the same 
face. The set contains a single standard region, un- 
less multiple standard regions have been aggregated 
according to the constructions of Remarks 4.2 and 
4.3. Set 

CTp = ^2 a R 

Ren(F) 

and 

t f= Yl tr - 

Ren(F) 

In [6, Section 4.4] a natural number n(R) is at- 
tached to the standard region R. Lower bounds on 
Tft(D) are obtained in [6, Theorem 4.4] that depend 



only on n(R). Write t n for the constants in this the- 
orem. The same theorem gives upper bounds s n for 
<Tft(D) that depend only on n(R). 

It is also helpful to bear in mind several bounds 
on angles in the following discussion. Every inter- 
nal angle of every standard region is at least 0.8638 
[5, 10.1.3]. Every internal angle of every standard 
region, except for triangles, is at least 1.153 [5, 4.3]. 

A decomposition star D determines a set of ver- 
tices U(D) that are of distance at most 2t from 
the center of D. Three consecutive vertices p\, P2, 
and P3 of a standard region are determined as the 
projections to the unit sphere of three corners v\, 
1)2, and 1)3, respectively in U(D). If the internal 
angle of a standard region is less than 1.32, then 
it follows by an interval arithmetic calculation that 
\ v i — v 3\ < V&, see [6, Lemma 3.11.4]. 

We are to the point in the proof of the Kepler 
Conjecture where it becomes necessary to use in- 
equalities proved by interval arithmetic in a serious 
way. To use these inequalities systematically, we 
combine inequalities into linear programs and solve 
the linear programs on computer. At first, our use 
of linear programs will be light, but our reliance will 
become progressively strong as this paper develops. 

To start out, we will make use of the inequal- 
ities [5, 10.3.4-8] and [5, 4.1]. These inequalities 
give lower bounds on tr(D) when R is a triangle 
or quadrilateral. To obtain lower bounds through 
linear programming, we take a linear relaxation. 
Specifically, we introduce a linear variable for each 
function tr and a linear variable for each internal 
angle aft. We substitute these linear variables for 
the nonlinear functions tr(D) and nonlinear inter- 
nal angle function into the given inequalities in [5]. 
Under these substitutions, the inequalities become 
linear. Given p triangles and q quadrilaterals at 
a vertex, we have the linear program to minimize 
the sum of the (linear variables associated with) 
Tft(D) subject to the constraint that the (linear 
variables associated with the) angles at the vertex 
sum to at most d. Linear programming yields a 
lower bound tlp(p, q, d) to this minimization prob- 
lem. This gives a lower bound to the corresponding 
constrained sum of nonlinear functions tr. 

Similarly, the inequalities from the same sections 
of [6] yield upper bounds ctlp (p, Q, d) on the sum of 
p + q functions aR, with p standard regions R that 
are triangular, and another q that are quadrilateral. 
These linear programs find their first application in 
the proof of the following proposition. 

Proposition 6.1. The plane graph of a contraven- 
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ing decomposition star satisfies Property 6 of tame- 
ness: 

^c(len(F)) > 8. 

F 

Proof. We will show that 

c{len{F))pt>a F {D) (4) 

Assuming this, the result follows for contravening 
stars D: 

Y, F c(len(F))pt >Ef*f(D) 
= o-(D) > 8pt. 

We consider three cases for Inequality 4. In the 
first case, assume that the face F corresponds to 
exactly one standard region in the decomposition 
star. In this case, Inequality 4 follows directly from 
the bounds [6, Theorem 4.4]. 

In the second case, assume we are in the context 
of a pentagon F formed in Remark 4.2. Then, again 
by [6, 4.4], we have 

<Jf(D) < (c(3) + c(8)) pt < c(5) pt. 

(Just examine the constants c(k).) 

In the third case, we consider the situation of 
Remark 4.3. The six standard regions score at most 
ss +cr LP (5,0,27r- 1.153) < c(8)pt. □ 

Proposition 6.2. Let F be a face of a contraven- 
ing plane graph G(D). Then 

t f {D) > d{len{F))pt. 

Proof. Similar. □ 

Lemma 6.3. Consider the aggregate in Remark 
4.3. 

1. There are at most two vertices (of a standard 
region) of degree six. If there are two, then 
they are adjacent vertices on a pentagon (Fig- 
ure 11 ). 

2. If a vertex of a pentagonal standard region has 
degree six, then the aggregate F of the six faces 
satisfies 

a F (D) < s 8 , 
t f {D) > t s . 

Proof. We begin with the second part of the lemma. 
The sum t f (D) over these six standard regions is 
at least 

h +t L p(5,0,2tt- 1.153) > t s . 




Figure 11: Adjacent vertices of degree 6 



Similarly, 

•55 + cr LP (5, 0, 2tt - 1.153) < s 8 . 

We note that there can be at most one exceptional 
region with a vertex of degree six. Indeed, if there 
are two, then they must both be vertices of the same 
pentagon: 

t s + h > (4ttC - 8) pt. 

Such a second vertex on the octagonal aggregate 
leads to one of the follow constants greater than 
(4-7r£— 8) pt. These same constants show that such a 
second vertex on a hexagonal aggregate must share 
two triangular faces with the first vertex of degree 
six. 

*8 +TLp(4,0,27r- 1.32 -0.8638), or 

t 8 +1.47p£+T LP (4,0,27T- 1.153- 0.8638), 

ts +r LP (5,0,2vr- 1.153). 

(The relevant inequalities involving the constants 
1.47 and 1.32 are found at [6, Lemma 5.12.1].) 

□ 

Proposition 6.4. A contravening plane graph sat- 
isfies Property 5 of tameness: If a vertex is con- 
tained in an exceptional face, then the degree of the 
vertex is at most 5. 

Proof. Assume for a contradiction that v is a ver- 
tex of degree six that lies in an exceptional region. 
There are several cases according to the number k 
of triangular regions at the vertex. 

(k < 2) If there are at least four nontriangular 
regions at the vertex, then the sum of interior angles 
around the vertex is at least 4(1.153) + 2(0.8638) > 
27r, which is impossible. 
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(k = 3) If there arc three nontriangular regions 
at the vertex, then t(D) is at least 2t 4 + t 5 + 
tlp(3, 0, 2tt - 3(1.153)) > (47< - 8) pt. 

(k = 4) If there are two exceptional regions at 
the vertex, then t(D) is at least 2is + tlp(4, 0, 2tt — 
2(1.153)) > (47rC-8)pi. 

If there are two nontriangular regions at the ver- 
tex, then t(D) is at least t 5 +Ti l p(4, 1, 2n- 1.153) > 
(4ttC - 8) pi. 

We are left with the case of five triangular regions 
and one exceptional region. 

When there is an exceptional standard region 
at a vertex of degree six, we claim that the ex- 
ceptional region must be a pentagon. If the re- 
gion is a heptagon or more, then t(D) is at least 
*7 + tlp(5, 0, 2tt - 1.153) > (4ttC - 8) pt. 

If the standard region is a hexagon, then t(D) 
is at least t 6 + tlp(5,0, 2tt — 1.153) > tg. Also, 
s 6 + °Xp(5, 0, 2-7T— 1.153) < sg. The aggregate of the 
six standard regions is a nonagon. The argument of 
[6, Section 4.6] extends to this context to give the 
bound of 8 pt. 

□ 



7 Weight Assignments 

The purpose of this section is to prove the existence 
of a good admissible weight assignment for contra- 
vening plane graphs. 

Theorem 7.1. Every contravening planar graph 
has an admissible weight assignment of total weight 
less than 14.8. 

Given a contravening decomposition star D, we 
define a weight assignment w(F) by 

w(F) = T F (D)/pt. 

Since D contravenes, 

J2 F w(F) =J2 F r F (D)/pt 

= t(D) < (4ttC - 8) pt/pt 
< 14.8. 

The challenge of the theorem will be to prove that 
w(F), when defined by this formula, is admissible. 

The next three lemmas establish that this defini- 
tion of w(F) for contravening plane graphs satisfies 
the first three defining properties of an admissible 
weight assignment. 



Lemma 7.2. Let F be a face of length n in a con- 
travening plane graph. Define w(F) as above. Then 
w(F) > d(n). 

Proof. This is Proposition 6.2. □ 

Lemma 7.3. Let v be a vertex of type (p, q) in a 
contravening plane graph. Define w{F) as above. 
Then 

J2w(F)>b(p,q). 

veF 

Proof. This is [5, Proposition 5.2]. □ 

Lemma 7.4. Let V be any set of vertices of type 
(5,0) in a contravening plane graph. Define w(F) 
as above. Lf the cardinality of V is k < 4, then 

W ( F ) > 0.55fc. 

vr\F ^$ 

Proof. This is [5, Proposition 5.3]. □ 

The following proposition establishes the final 
property that w(F) must satisfy to make it admis- 
sible. 

Proposition 7.5. Let V be any separated set of 
vertices in a contravening plane graph. Define w(F) 
as above. Then 

(w(F) - d(len(F))) > ]T a(tri(v)), 

VnF^V) vGV 

where tri(v) denotes the number of triangles con- 
taining the vertex v. 

The proof will occupy the rest of this section, as 
well as all of Section 8. Since the degree of each 
vertex is five, and there is at least one face that 
is not a triangle at the vertex, the only constants 
tri(v) that arise are 

tri(v) e {0, ... ,4} 

We will prove that in a contravening plane graph 
that the Properties (1) and (4) of a separated set 
are incompatible with the condition tri(v) < 2, for 
some v eV. 

This will allows us to assume that 

tri{v) e {3,4}, 

for all v eV. These cases will be treated in Sections 
7.3 and 7.4. 
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The outline of the proof of Proposition 7.5 is the 
following. Let D be a contravening decomposition 
star. First, assume that every face meeting V cor- 
responds to a single standard region of D. (That 
is, we exclude the aggregates created in Remark 4.2 
and Remark 4.3.) Under this hypothesis, we show 
that tri(v) < 2 is incompatible with the hypothesis. 
This allows us to assume that tri(v) € {3,4}. We 
express the desired inequality as a sum of inequali- 
ties indexed by V, where the inequality indexed by 
v e V depends only on (an explicit description of) 
the function t r {D) for standard regions R corre- 
sponding to a face containing v. These individual 
inequalities are proved first for tri(v) — 3, and then 
for tri(v) = 4. 

Finally in Section 8, we return to the aggregates 
of Remark 4.2 and Remark 4.3 and show that the 
conclusions can be extended to aggregate faces as 
well. 

7.1 Proof that tri(v) > 2 

In this subsection I? is a contravening decomposi- 
tion star with associated graph G(D). Let V be a 
set of vertices that satisfies the conditions of Propo- 
sition 7.5. Let v be a vertex in V such that none 
of its faces is an aggregate in the sense of Remarks 

4.2 and 4.3. 

Lemma 7.6. Under these conditions, for every v <G 
V, tri(v) > 1. 

Proof. If there are p triangles, q quadrilaterals, and 
r other faces, then 

r(D) >E veR MD) 

> rt 5 + T L p(p,q,27r - r(1.153)). 

If there is a vertex w that is not on any of the faces 
containing v, then the sum of tf{D) over the faces 
containing w yield an additional 0.55 pt by [5, 5.2]. 
We run the linear programs for each (p, q, r) and 
find that the bound is always greater than (47r£ — 
8) pt. This implies that D cannot be contravening. 



(P, Q, r) 

(0,5,0) 
(0,9, r > 1) 
(1,4,0) 
(1,3,1) 
(1,2,2) 
(l,«,r>3) 



lower bound 

22.27 pt [5, 5.2] 
t 5 + 4t 4 w 14.41 pt 

17.62 pt [5, 5.2] 
h + 12.58 pt (r LP ) 
2t 5 + 7.53 pt (tlp) 

3^5 + ti 



□ 



Lemma 7.7. Under these same conditions, for ev- 
ery v G V, tri(v) > 2. 

Proof. Assume that tri(v) = 2. We will show that 
this implies that D is not contravening. Let e be 
the number of exceptional faces at v. 

The constants 0.55 pt and 4.52 pt used throughout 
the proof come from [5, 5.3] ; The constants t n comes 
from [6, 4.4]. 

(e = 3): First, assume that there are three ex- 
ceptional faces around vertex v. They must be 
pentagons (2t 5 + f 6 > (4n( — 8)pt). The aggre- 
gate of the five faces is an m-gon (some m < 11). 
If there is a vertex not on this aggregate, use 
3t 5 + 0.55 pt > (4-7rC — 8) pt. So there are at most 
nine triangles away from the aggregate. 

a(D) <9pt+ (3s 5 + 2pt) < 8pt. 

The argument if there is a quad, pentagon, and 
hexagon is the same (£4 + 1§ = 2<5, S4 + sq = 2s§). 

(e = 2): Assume next that there are two pen- 
tagons and a quadrilateral around the vertex. The 
aggregate of the two pentagons, quadrilateral, and 
two triangles is an m-gon (some m < 10). There 
must be a vertex not on the aggregate of five faces, 
for otherwise we have 

a(D) < 8 pt + (2s 5 + 2 pt) < 8 pt. 

The interior angle of one of the pentagons is at 
most 1.32. For otherwise, r LP (2, 1, 2tt - 2(1.32)) + 
2t 5 + 0.55pt> (4nC-8)pt. 

The calculation [6, Lemma 5.12.2] shows that any 
pentagon R with an internal angle less than 1.32 
yields tr(D) > 5.66 pt. If both pentagons have 
an interior angle < 1.32 the lemma follows easily 
from this calculation: 2(5.66) pt + tlp(2, 1, 2ir — 
2(1.153)) + 0.55 > (4ttC - 8) pt. If there is 
one pentagon with angle > 1.32, we then have 
5.66p£+T LP (2, 1, 2tt- 1.153- 1.32) +t 5 + 0.55 pt> 
(4ttC - 8) pt. 

(e = 1): Assume finally that there is one ex- 
ceptional face at the vertex. If it is a hexagon (or 
more), we are done t% + tlp(2,2,27t — 1.153) > 
(47rC — 8) pt. Assume it is a pentagon. The ag- 
gregate of the five faces at the vertex is bounded by 
an m-circuit (some m < 9). If there are no more 
than 11 quasi- regular tetrahedra outside the aggre- 
gate, the score is at most (1 + 2(4.52)) pt + s 5 + 
cr L p(2,2,27r - 1.153) < 8pt ([5, 5.3]). So we may 
assume that there are at least three vertices not on 
the aggregate. 
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If the interior angle of the pentagon is greater 
than 1.32, we have 

r LP (2, 2, 2tt- 1.32) + 3(0.55) pt +t 5 > (4t< - 8) pt; 

and if it is less than 1.32, we have by [6, Lemma 
5.12.1] 

tlp(2, 2, 2tt - 1.153) +3(0.55)pt+ lA7pt+t 5 
> (4tt£ - 8) pt. 

□ 

Lemma 7.8. The bound tri(v) > 2 holds if v is a 
vertex of an aggregate face. 

Proof. The exceptional region enters into the pre- 
ceding two proofs in a purely formal way. Pentagons 
enter through the bounds 

£5, s 5 , 1.47 pt, 5.66pt 

and angles 1.153, 1.32. Hexagons enter through the 
bounds 

and so forth. These bounds hold for the aggregate 
faces. Hence the proofs hold for aggregates as well. 

□ 

7.2 Bounds when tri(v) G {3,4}. 

In this subsection D is a contravening decomposi- 
tion star with associated graph G(D). Let V be a 
separated set of vertices. For every vertex v in V, 
we assume that none of its faces is an aggregate in 
the sense of Remarks 4.2 and 4.3. We assume that 
there are three or four triangles containing v, for 
every v eV. 

To prove the Inequality 4 in the definition of ad- 
missible weight assignments, we will rely on the fol- 
lowing reductions. Define an equivalence relation 
on exceptional faces by F ~ F' if there is a se- 
quence Fo — F,...,F r = F' of exceptional faces 
such that consecutive faces share a vertex of type 
(3, 0, 2). Let T be an equivalence class of faces. 

Lemma 7.9. Let V be a separated set of vertices. 
For every equivalence class of exceptional faces T , 
let V(F) be the subset ofV whose vertices lie in the 
union of faces of T . Suppose that for every equiv- 
alence class T, the Inequality 4 (in the definition 
of admissible weight assignments) holds for V(T). 
Then the Inequality holds for V . 



Proof. By construction, each vertex in V lies in 
some F, for an exceptional face. Moreover, the sep- 
arating property of V insures that the triangles and 
quadrilaterals in the inequality are associated with 
a well-defined T . Thus, the inequality for V is a 
sum of the inequalities for each V(F). □ 

Lemma 7.10. Let v be a vertex at which there are 
p triangles, q quadrilaterals, and r other faces. Sup- 
pose that for some p' < p and q' < q, we have 

T LP (p', q', a) > {p'd(3) + q'd{A) + o(p)) pt 

for some upper bound a on the angle occupied by p' 
triangles and q' quadrilaterals at v. Suppose further 
that the Inequality 4 (in the definition of admissible 
weight assignments) holds for V' = V \ {v}. Then 
the inequality holds for V. 

Proof. Let Fi, . . . , F m , m = p' + q', be faces corre- 
sponding to the triangles and quadrilaterals in the 
lemma. The hypotheses of the lemma imply that 

m 

J2(r Fi (D) - d(len(Fi))) > o(p). 
i 

Clearly, the Inequality for V is the sum of this in- 
equality, the inequality for V, and d(n) > 0. □ 

7.3 tri(v) = 3 

To continue further, we need to recall from [6] how 
tr(D), for exceptional standard regions, can be 
written as a sum of terms in a way that reflects 
the internal geometry of D over R. We have pre- 
pared a summary of the relevant material from [6] 
in Appendix 16. We recommend that the reader re- 
view Appendix 16 before reading the next lemma. 
We use terminology from that appendix (penalties, 
loops, erasure, flat quarters, D(n, k), 5i oop (n, k), a). 

Lemma 7.11. Property 4 °f admissibility holds 
when tri(v) = 3 for all v G V. 

Proof. When there are three triangles at a vertex of 
degree five along with an exceptional face, the rele- 
vant constant from the definition of admissibility is 
a(3) = 1.4. Erase all upright quarters except those 
in loops. The result follows from [6, Lemma 5.12.1]. 
and the subsequent remark if the internal angle at 
the vertex of an exceptional region is < 1.32. 

Assume the angles on the exceptional regions are 
> 1.32. If there are three triangles, a quad, and 
an exceptional face at the vertex, then we have 
tlp(3, 1,2tt- 1.32) > lApt+U. 
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The final case is three triangles and two excep- 
tional faces at the vertex. There can be no hep- 
tagons and at most one hexagon 2t e = t 5 + t 7 > 
(47r£ — 8) pt. In these calculations, the flat quar- 
ters are scored by a. By Appendix 16, each flat 
quarter contributes -0(3, 1) to t n . The constant 
a(tri(v)) pt estimates what the flat quarter and sur- 
rounding quasi-regular tetrahedra squander in ex- 
cess of 0(3, 1). Most of the calculations now follow 
directly from the interval calculations of Appendix 
17.20. A few additional comments are needed if 
there is a loop masking a flat quarter. If the con- 
stant 5i oop (n,k) exceeds r(lA)pt, where r is the 
greatest possible number of masked flat quarters 
within the loop, then the result follows. This gives 
every case except (n,k) — (4,1), (5,0). If (n,k) — 
(4, 1), we erase and take the penalty ttf = 2£y + £r 
whenever the internal angle is at least 1.32. This 
penalty appears in Appendix 17.20. 

If (n,k) = (5,0), there are various possibilities. 
If there is one masked flat quarter, use Si oop (5, 0) > 
1.4 pt. If there arc two masked flat quarters, erase 
with penalty ttf = 2£y + £r/2 at each flat quarter. 
If the internal angle is at most 1.32, use 

3.07 pf- .0(3,1) > lApt+ir F . 

If the internal angle is at least 1.32, use the interval 
calculations Appendix 17.20 with penalties. □ 

7.4 tri(v) = 4 

If there are four triangles and one exceptional re- 
gion at a vertex, we justify the constant a(tri(A)) = 

1.5 pt, in the definition of admissible weight assign- 
ment. These calculations are similar to those of Sec- 
tion 7.3 and are based on the interval calculations 
Appendix 17.21. Details are left to the reader. It 
is simpler because there is only one exceptional re- 
gion, and it is not necessary to separate the case 
according to the size of the internal angle (at least 
1.32 and at most 1.32). 



8 The Aggregate Cases 

8.1 The non-polygonal standard re- 
gions 

Lemma 8.1. Consider a separated set of vertices 
V on an aggregated face F as in Remark 4-2. The 
Inequality 4 holds (in the definition of admissible 



weight assignments): 

(w(F) - d(len(F))) > a(tri(v)). 

Proof. We may assume that tri(v) £ {3,4}. 

First consider the aggregate of Remark 4.2 of 
a triangle and eight-sided region, with pentagonal 
hull F. There is no other exceptional region in a 
contravening decomposition star with this aggre- 
gate: 

t 8 + t 8 > (4ttC - 8) pt. 

A separated set of vertices V on F has cardinality 
at most 2. This gives the desired bound 

t 8 > h + 2(1.5) pt. 

Next, consider the aggregate of a hexagonal hull 
with an enclosed vertex. Again, there is no other 
exceptional face. If there are at most k < 2 vertices 
in a separated set, then the result follows from 

t 8 > t 6 + k{1.5)pt. 

There are at most three vertices in V on a hexagon, 
by the non- adjacency conditions defining V. A ver- 
tex v can be removed from V if it docs not lie on 
a flat quarter (Lemma 7.10 and Appendix 17.21). 
However, if there is an enclosed vertex w, it is im- 
possible for there to be three nonadjacent vertices, 
none lying along a flat quarter: 

£(2, 2, 2, V8, V8, V8, 2t , 2t , 2) > 2t . 

(The definition of £ can be reviewed at [2, Sec- 
tion 1].) 

Finally consider the aggregate of a pentagonal 
hull with an enclosed vertex. There are at most 
k < 2 vertices in a separated set in F. There is no 
other exceptional region: 

t 7 + h > (4ttC - 8) pt. 

The result follows from 

t 7 > t 5 + 2(1.5) pt. 

□ 

Lemma 8.2. Consider a separated set of vertices 
V on an aggregate face of a contravening plane 
graph as in Remark 4.3. The Inequality 4 holds in 
the definition of admissible weight assignments. 
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Proof. There is at most one exceptional face in the 
plane graph: 

t s + t 5 > (4vrC - 8) pt. 

Assume first that aggregate face is an octagon (Fig- 
ure 10). At each of the vertices of the face that lies 
on a triangular standard region in the aggregate, we 
can remove the vertex from V using Lemma 7.10 
and the estimate 

t lp (4,0,27t- 2(0.8638)) > 1.5 pt. 

This leaves at most one vertex in V, and it lies on a 
vertex of F which is "not aggregated," so that there 
are five standard regions of the underlying decom- 
position star at that vertex, and one of those regions 
is pentagonal. The value o(4) = 1.5 pt can be esti- 
mated at this vertex in the same way it is done for 
a non-aggregated case. (We use the calculations of 
Appendix 17.21.) 

Now consider the case of an aggregate face that is 
a hexagon (Figure 10). The argument is the same: 
we reduce to V containing a single vertex, and argue 
that this vertex can be treated by the inequalities 
of Appendix 17.21. □ 

9 Four-cycles 

It has now been verified that all the properties of 
tameness except one have been verified. This sec- 
tion verifies that contravening plane graphs satisfy 
this final property, and in this way we complete the 
proof that all contravening plane graphs are tame. 

9.1 Possible four-circuits 



Proposition 9.1. Let D be a contravening plane 
graph. The interior of a 4-circuit surrounds is one 
of the aggregates of faces shown in Property 3 of 
tameness. 

Proof. If there are no enclosed vertices, then the 
only possibilities are for it to be a single quadrilat- 
eral face or a pair of adjacent triangles. 

Assume there is one enclosed vertex v. If v is con- 
nected to 3 or 4 vertices of the quadrilateral, then 
that possibility is listed as part of the conclusion. 

If v is connected to 2 opposite vertices in the 
4-cycle, then the vertex v has type (0, 2) and the 
bounds of [5, Proposition 5.2] show that the graph 
cannot be contravening. 

If v is connected to 2 adjacent vertices in the 4- 
cycle, then we appeal to Lemma 9.3 to conclude 
that the graph is not contravening. 

If v is connected to or 1 vertices, then we appeal 
to [5, Lemma 2.2], which shows that this creates an 
impossible geometric situation. This completes the 
proof. □ 

9.2 A non-contravening Four-circuit 

This subsection rules out the existence of a par- 
ticular four-circuit on a contravening plane graph. 
The interior of the circuit consists of two faces: a 
triangle and a pentagon. The circuit and its en- 
closed vertex are show in Figure 12 with vertices 
marked p\ , . . . ,p§. The vertex p\ is the enclosed 
vertex, the triangle is (pi,P2,Pb) and the pentagon 
is (pi, . . . ,p 5 ). 

2 5 



Every 4-circuit divides a plane graph into two ag- 
gregates of faces that we may call the interior and 
exterior. We call vertices of the faces in the aggre- 
gate that do not lie on the 4-cycle enclosed vertices. 
Thus, every vertex lies in the 4-cycle, is enclosed 
over the interior, or is enclosed over the exterior. 

Proposition [3, 4.2] asserts that either the interior 
or the exterior has at most 1 enclosed vertex. (The 
proof does not require the hypothesis that the graph 
be contravening. It expresses a geometric constraint 
holding for all decomposition stars.) When choos- 
ing which aggregate is to be called the interior, we 
may make our choice so that the interior contains at 
most 1 enclosed vertex. With this choice, we have 
the following proposition. 



Figure 12: A non-contravening 4-circuit 

Suppose that D is a decomposition star whose as- 
sociate graph contains such a 4-circuit. Recall that 
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D determines a set U (D) of vertices in Euclidean 
3-space of distance at most 2i from the center (the 
origin) of D, and that each vertex pi can be real- 
ized geometrically as a point on the unit sphere at 
the origin, obtained as the radial projection of some 
Vi € U(D). 

In the rest of this subsection, we assume that 
the reader is familiar with techniques in [6]. (In 
particular, the reader will need to be familiar with 
geometric considerations and the details about how 
the function tr(D) is estimated for a pentagonal 
standard region R. This is reviewed in Appendix 
16.) 

Lemma 9.2. One of the edges (vi,v 3 ), (vi,Vi) has 
length less than 2^/2. Both of the them have lengths 
less than 3.02. Also, \vi\ > 2.3. 

Proof. This is a standard exercise in geometric 
considerations ([2, Section 1] and [5, Section 2]). 
(The reader should review these sections for the 
framework of the following argument.) We de- 
form the figure using pivots to a configuration 
V2, ■ ■ ■ , i>5 at height 2, and \vi — Vj\ — 2to, = 
(2,3), (3,4), (4,5), (5,2). We scale vi until \ Vl \ = 
2to- We can also take the distance from v\ to v^ 
and to i>2 to be 2. If we have \v\ — > 2V2, then 
we stretch the edge \v\ — Vi\ until \v± — v%\ = 2\[2. 
The resulting configuration is rigid. Pick coordi- 
nates to find that \v\ — va\ < 2\f2. If we have 
l^i — ^3 1 > 2*o, follow a similar procedure to re- 
duce to the rigid configuration \v\ — v-&\ — 2fo, to 
find that \v\ — v±\ < 3.02. The estimate \vi\ > 2.3 
is similar. □ 

There arc restrictive bounds on the dihedral 
angles of the simplices (0, v\, Vi, Vj) along the 
edge (0, Vi). The quasi- regular tetrahedron has 
a dihedral angle of at most 1.875 ([5, 10.1.2]). 
The dihedral angles of the simplices (0, v\, i>2, V3), 
(0, vi, 1)5, V4) adjacent to it are at most 1.63. The di- 
hedral angle of the remaining simplex (0, v±, i> 3 , V4) 
is at most 1.51. (These constants come from 1.624 
and 1.507 in Appendix 17.8.) This leads to lower 
bounds as well. The quasi-regular tetrahedron has a 
dihedral angle that is at least 2n — 2(1.63) — 1.51 > 
1.51. The dihedral angles adjacent to the quasi- 
regular tetrahedron is at least 2n — 1.63 — 1.51 — 
1.875 > 1.26. The remaining dihedral angle is at 
least 2tt - 1.875 - 2(1.63) > 1.14. 

Lemma 9.3. A decomposition star that contains 
this configuration is not contravening: t(D) > 
(4ttC - 8) pt. 



Proof. We will show that the sum of tr(D), as R 
ranges over the pentagonal standard region and tri- 
angular regions, is at least 11.16 pi. Let P denote 
the aggregate formed by these standard regions. 

First we show how the lemma follows from this 
bound. There arc no other exceptional faces, be- 
cause 11.16 pt+ts > (47r( — 8) pt. Every vertex not 
on P has type (5, 0), by [5, Proposition 5.2]. In par- 
ticular, there are no quadrilateral regions. There 
are at most 4 triangles at every vertex of P, be- 
cause the tlp(5,0,27t- 1.32) > 6.02 pt. (The 1.32 
comes from that fact that P has no flat quarters be- 
cause v is too short to be enclosed over one [2, 1.3] 
[6, 3.11.4].) There are at least 3 triangles at every 
vertex of P, otherwise we contradict a property of 
tameness (Property 2 and Property 3). 

The only triangulation with these properties is 
obtained by removing one edge from the icosahe- 
dron (Exercise). This implies that there are two 
opposite corners of P each having four quasi- regular 
tetrahedra. Since the diagonals of P have lengths 
greater than 2\/2, the results of Section 7.4 show 
that these eight quasi-regular tetrahedra squander 
at least 2(1.5) pi. There are two additional ver- 
tices of type (5,0) whose tetrahedra are distinct 
from these eight quasi-regular tetrahedra. They 
give an additional 2(0.55) pt Now (11.16 + 2(1.5) + 
2(0.55)) pt > (4ttC - 8) pt. The result follows. 

Return to the bound 11.16 pi. (A review of Ap- 
pendix 16 is recommended before continuing. The 
definitions of r, f and tq are reviewed in Appendix 
15.) Let dih denote the dihedral angle of a sim- 
plex along a given edge. Let Sij be the simplex 
(0, Vl ,Vi,Vj), for = (2, 3), (3, 4), (4, 5), (2, 5). 

We have dih(Sij) = 2tt. Suppose one of the 

edges (^1,^3) or (1)1,1)4) has length > 2\/2. Say 
{vi,v 3 ). 

We have (by Appendix 17.1) 

t(S 25 ) -0.2529 dih(5 25 ) > -0.3442, 

r (5 23 ) -0.25 29 dih(5 23 ) > -0.1787, 

t(S* 45 ) -0.25 29 dih(5 45 ) > -0.2137, 

to^) -0.2529 dih(5 34 ) > -0.1371. 

We have by [6, A10] that 

r(D) >£(4)^(^-)-5£r 

> 2tt(0.2529) - 0.3442 

-0.1787 - 0.2137 - 0.1371 - 5£r 

> 11.16 pf, 

where t x = t, f, r as appropriate. 

Now suppose (vi,v 3 ) and (i>i,i> 4 ) have length < 
2\/2. If there is an upright diagonal that is not 
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enclosed over either flat quarter, the penalty is at 
most 3£r + 2£y. Otherwise, the penalty is smaller: 
4£r + We have 

r(D) >E W r(Sn)-(Mr + ^ v ) 

> 2tt(0.2529) - 0.3442 

-2(0.2137) - 0.1371 - (3£ r + 2£ y ) 

> 11.16 pi. 

□ 

The proof that contravening plane graphs are 
tame is complete. 



10 Linear Programs 

We have completed the first half of the proof of the 
Kepler Conjecture by proving that every contraven- 
ing plane graph is tame. 

The second half of the proof of the Kepler Con- 
jecture consists in showing that tame graphs are 
not contravening, except for the isomorphism class 
of graphs isomorphic to G/ cc and Ghcp associated 
with the face-centered cubic and hexagonal close 
packings. 

The paper [5] treats all contravening tame graphs 
in which every face has length three or four, except 
for three cases Gf cc , G pen t, and Ghcp- The two 
cases Gf cc and Ghcp arc treated in [2, 3.13], and 
the case G pm t is treated in [7]. 

We complete the proof of the Kepler Conjecture 
by analyzing all contravening plane graphs that con- 
tain at least one exceptional face. The primary tool 
that will be used is linear programming. The linear 
programs are obtained as relaxations of the origi- 
nal nonlinear optimization problem of maximizing 
<j{D) over all decomposition stars whose associated 
graph is a given tame graph G. The upper bounds 
obtained through relaxation are upper bounds to 
the nonlinear problem. 

To eliminate a tame graph, we must show that it 
is not contravening. By definition, this means we 
must show that a{D) < 8pt. When a single linear 
program does not yield an upper bound under 8 pt, 
we branch into a sequence of linear programs that 
collectively imply the upper bound of 8 pt. This will 
call for a sequence of increasingly complex linear 
programs. 



10.1 Relaxation 

(NLP) Let / : P — > R be a function on a nonempty 
set P. Consider the nonlinear maximization prob- 
lem 

max/O). 

pep 

(LP): Consider a linear programming problem 
max c • x 

such that Ax < b, where A is a matrix, b, c are vec- 
tors of real constants and a; is a vector of variables 
x = (xi, . . . ,x n ). We write the linear programming 
problem as 

max(c ■ x : Ax < b). 
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An interpretation I of a linear programming 
problem (LP) is a set |J|, together with an assign- 
ment I of functions x\ : | J| — ► R to variables 
Sj. We say the constraints < 6 of the linear 
program are satisfied under the interpretation I if 
for all p e |J|, we have 

Aai^p) < b. 

The interpretation / is said to be a relaxation of 
the nonlinear program (NLP), if the following three 
conditions hold. 

1. P=\I\. 

2. The constraints are satisfied under the inter- 
pretation. 

3. f(p) <c-x T {p), for all pg |7|. 

Lemma 10.1. Let (LP) be a linear program with 
relaxation L to (NLP). Then (LP) has a feasible 
solution. Moreover, if (LP) is bounded above by 
a constant M, then M is an upper bound on the 
function f : \I\ — > R. 

Proof. A feasible solution is Xi = x\ (p) , for any 
p£ |/|. The rest is clear. □ 

Remark 10.2. In general, it is to be expected that 
the interpretations Ax 1 < b will be nonlinear in- 
equalities on the domain P. In our situation, sat- 
isfaction of the inequality will be proved by inter- 
val arithmetic. Thus, the construction of an upper 
bound to (NLP) breaks into two tasks: to solve the 
linear programs and to prove the nonlinear inequal- 
ities required by the interpretation. 

There are many nonlinear inequalities that enter 
into our interpretation, which have been proved by 
interval arithmetic on computer. These inequalities 
are contained in Appendix 17. 

Remark 10.3. There is a second method of estab- 
lishing the satisfaction of inequalities under an in- 
terpretation. Suppose that we wish to show that 
the inequality e • x < b' is satisfied under the in- 
terpretation L. Suppose that we have already es- 
tablished that a system of inequalities Ax < b is 
satisfied under the interpretation L. We solve the 
linear programming problem max(e ■ x : Ax < b). 
If this maximum is at most b', then the inequality 
e • x < b' is satisfied under the interpretation I. We 
will refer to e • x < b' as an LP-derived inequality 
(with respect to the system Ax < b). 



10.2 The Linear Programs 

Let G be a tame plane graph containing an excep- 
tional face. Let Xg be the space of all decomposi- 
tion stars whose associated plane graph is isomor- 
phic to G. 

Theorem 10.4. For every tame plane graph G 
containing at least one exceptional face, there ex- 
ists a finite sequence of linear programs with the 
following properties. 

1. Every linear program has an admissible solu- 
tion and its solution is strictly less than 8 pt. 

2. For every linear program in this sequence, there 
is an interpretation I of the linear program that 
is a relaxation of the nonlinear optimization 
problem 

a : \I\ — > R, 
where \I\ is a subset of Xq- 

3. The union of the subsets \L\, as we run over 
the sequence of linear programs, is Xg- 

The proof is constructive. For every tame plane 
graph G a sequence of linear programs is generated 
by computer and solved. The optimal solutions are 
all bounded above by 8pt. It will be clear from 
construction of the sequence that the union of the 
sets 1 7" | exhausts Xg- We estimate that nearly 10 5 
linear programs are involved in the construction. 

Remark 10.5. The paper [9, Section 3.1.1] shows 
how the linear programs that arise in connection 
with the Kepler Conjecture can be formulated in 
such a way that they always have a feasible solu- 
tion and so that the optimal solution is bounded. 
We assume that all our linear programs have been 
constructed in this way. 

Corollary 10.6. If a tame graph G contains an ex- 
ceptional face, then it is not contravening. 

Proof. This follows immediately from the theorem 
and Lemma 10.1. □ 



11 LP Formulation 

11.1 Basic Linear Programs 

To describe the basic linear program, we need the 
following indexing sets. 
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Let VERTEX be the set of all vertices in G. Let 
FACE be the set of all faces in G. (Recall that by 
construction each face F of the graph carries an 
orientation.) Let ANGLE be the set of all angles in 
G, defined as the set of pairs (v,F), where the ver- 
tex v lies in the face F. Let DIRECTED be the set 
of directed edges. It consists of all ordered pairs 
(v,s(v,F)), where s(v,F) denotes the successor of 
the vertex v in the oriented face F. Let TRIANGLES 
be the subset of FACE consisting of those faces of 
length 3. Let UNDIRECTED be the set of undirected 
edges. It consists of all unordered pairs {v, s(v, F)}, 
for v € F. 

We introduce variables indexed by these sets. 
Following AMPL notation, we write for instance 
y{VERTEX} to declare a collection of variables y[v] 
indexed by vertices v in VERTEX. With this in mind, 
we declare the variables 

a{ ANGLE}, y{VERTEX}, e{UNDIRECTED}, 
<t{FACE}, t{FACE}, sol{FACE}. 

We obtain an interpretation / on the compact 
space Xq. First, we define an interpretation at the 
level of indexing sets. A decomposition star deter- 
mines the set U(D) of vertices of height at most 
2to from the origin of D. Each decomposition star 
D e Xq determines a (metric) graph with geodesic 
edges on the surface of the unit sphere, which is 
isomorphic to G as a (combinatorial) plane graph. 
There is a map from the vertices of G to U (D) given 
by v i— > v 1 , if the radial projection of v 1 to the unit 
sphere at the origin corresponds to v under this iso- 
morphism. Similarly, each face F of G corresponds 
to a set F 1 of standard regions. Each edge e of 
G corresponds to a geodesic edge e 1 on the unit 
sphere. 

Now we give an interpretation I to the linear pro- 
gramming variables at a decomposition star D. As 
usual, we add a superscript / to a variable to indi- 
cate its interpretation. Let a[v,FY be the sum of 
the internal angles at v 1 of the metric graph in the 
standard regions F 1 . Let y[vY be the length \v I \ 
of the vertex v 1 e U(D) corresponding to v. Let 
e[v, w] 1 be the length |u 7 — w I \ of the edge between 
v 1 and w 1 G U(D). Let 

a[FY = a F (D), 
sol[F] 7 = soliF 1 ), 
r[FY =r F (D). 

The objective function for the optimization prob- 
lems is 

max : ^ uf- 
fgFACE 



Its interpretation under I is the score a(D). 

We can write a number of linear inequalities that 
will be satisfied under our interpretation. For ex- 
ample, we have the bounds 

< y[v] < 2t , v € VERTEX 

< e[v, w] < 2t (v, w) e EDGE 

< a[v, F] < 2tt (v, F) e ANGLE 

< sol[F] < 4tt F e FACE 

There are other linear relations that are suggested 
directly by the definitions or the geometry. Here, v 
belongs to VERTEX. 

t[F] = sol[F](pt- a[F], 

sol[F] = Y,veF <*b>, F] - {len(F) - 2)tt 

There are long lists of additional inequalities that 
come from interval arithmetic verifications. For ex- 
ample, the interval verifications of [2, Propositions 
4.1-4.3] supply us with many linear inequalities re- 
lating the variables 

a[F], T [F], a[v,F], 
sol[F], y[v], e[v,w] 

whenever F 1 is a single standard region having 
three sides. Similarly, [5, Section 4] gives inequali- 
ties for cr[F] and related variables, when the length 
of F is four. A complete list of inequalities that are 
used for triangular and quadrilateral faces is found 
in [5, Section 9]. 

For exceptional faces, we have an admissible 
weight function w(F). According to definitions 
w(F) = r[F]/pt, so that the inequalities for the 
weight function can be expressed in terms of the 
linear program variables. 

When the exceptional face is not an aggregate, 
then it also satisfies the inequalities of [6, Theorem 
4.4]. 

11.2 Face Refinements 

The linear program formulated in the previous sec- 
tion does not carry enough information in general 
to obtain the desired bounds. The paper [6] shows 
that exceptional standard regions have a great deal 
of internal structure. Until now in this paper, we 
have mostly ignored this internal structure. We will 
now make a detailed examination of this internal 
structure. 

A refinement F of a face F of a plane graph G is 
a set F of faces such that 
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1. The intersection of the vertex set of G with 
that of F is the set F. 

2. FU {F°p} is a plane graph. 

We use refinements of faces to describe the internal 
structure of faces. 

We introduce indexing sets FACE-F, VERTEX-F, 
ANGLE-F, EDGE-F, the sets of faces, vertices, angles, 
and edges in F, respectively, analogous to those in- 
troduced for G. 

We create variables tt[F], and indexed variables 

sol{FACE-F}, sc{FACE-F} tsc{FACE-F}, 
a{ ANGLE-F}, y{ VERTEX-F}, e{EDGE-F}. 

(Variables with names and u e[v, w]" were al- 

ready created for some v, w g VERTEX-F n VERTEX. 
In these cases, we use the variables already created.) 

Each vertex v in the refinement will be inter- 
preted either as a vertex v 1 g 11(D), or as the end- 
point of an upright diagonal lying over the standard 
region F 1 . We will interpret the faces of the refine- 
ment in terms of the geometry of the decomposition 
star D variously as flat quarters, upright quarters, 
anchored simplices, and the other constructs of [6]. 
This interpretation depends on the context, and will 
be described in greater detail below. 

Once the interpretation of faces is fixed, the inter- 
pretations are as before for the variable names in- 
troduced already: y, e, a, sol. The lower and upper 
bounds for a and sol are as before. The lower and 
upper bounds for y[v] are 2 and 2i if v 1 g U(D), 
but if (0, v 1 ) is an upright diagonal, then the bounds 
are [2to,2\/2]. The lower and upper bounds for e 
will depend on the context. 

11.3 Variables related to score 

The variables sc are a stand-in for the score a on 
a face. We do not call them a because the sum of 
these variables will not in general equal the variable 
cr[F], when F is a refinement of F: 

sc[F']^a[F]}. 

F'eF 

We will use have a weaker relation: 

a[F}< J2 sc[F'] +7 r[F]. 
F'eF 

The variable 7r[F] is called the penalty associated 
with the refinement F. (Penalties are discussed at 



length in Appendix 16.) The interpretations of sc 
and 7r[F] are rather involved, and will be discussed 
on a case-by-case basis below. The interpretation 
of rsc follows from the identity: 

rsc[F'] = sol[F'](pt - sc[F'], VF' g F. 

The interpretation of variables that follows might 
appear to be hodge-podge at first. However, they 
are obtained in a systematic way. We analyze the 
proofs and approximations in [6], and define scfF] 7 
as the best penalty-free scoring approximation that 
is consistent with the given face refinement, here 
are the details. 

If the subregion is a flat quarter, the interpreta- 
tion of sc[F] is the function a, defined in Appendix 
16. If the subregion is an upright quarter Q, the in- 
terpretation of sc[F] is the function a(Q) from [2]. If 
the subregion is an anchored simplex that is not an 
upright quarter, sc[F] is interpreted as the analytic 
Voronoi function vor if the simplex has type C or 
C , and as voro otherwise. (The types A, B, C and 
C are defined in [6, Sections 2.5-2.10].) Whether or 
not the simplex has type C, the inequality sc[F] < 
is satisfied. In fact, if vor scoring is used, we note 
that there arc no quoins, and <j)(\, t ) < 0. 

If the subregion is triangular, if no vertex repre- 
sents an upright diagonal, and if the subregion is 
not a quarter, then sc[F] is interpreted as vor or 
voro depending on whether the simplex has type A. 
In either case, the inequality sc[F] < vor is satis- 
fied. 

In most other cases, the interpretation of sc[F] is 
vor . However, if R is a heptagon or octagon, and 
F has > 4 sides, then sc[F] is interpreted as vor 
except on simplices of type A, where it becomes the 
analytic Voronoi function. 

If R is a pentagon or hexagon, and F is a quadri- 
lateral that is not adjacent to a flat quarter, then 
the interpretation of sc[F] is the actual score of the 
subregion over the subregion. In this case, the score 
<jr has a well-defined meaning for the quadrilateral, 
because it is not possible for an upright quarter in 
the Q-system to straddle the quadrilateral region 
and an adjacent region. Consequently, any erasing 
that is done can be associated with the subregion 
without ambiguity. By the results of [4], we have 
sc[F] < 0. We also have sc[F] < vor . 

One other bound that we have not explicitly men- 
tioned is the bound (Jr{D) < s n . For heptagons 
and octagons that are not aggregates, this is a bet- 
ter bound than the one used in the definition of 
tameness (Property 6). In heptagons and octagons 
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that are not aggregates, if we have a subregion 
with four or more sides, then sc[F] < Z(n, k) and 
tsc[F] > D(n,k). (See [6, 5.5.1], [6, 5.5.2].) 

The variables are subject to a number of compati- 
bility relations that are evident from the underlying 
definitions and geometry. 

sol[F'] = £„ eF , a[v, F'} - (len[F'\ - 2)tt, VF' 
E s<at[v,F'] = a[v,F], Vw 

Assume that a face F\ e F has been interpreted 
subregion R = F[ of a standard region. As- 
sume that each vertex of F\ is interpreted as a 
vertex in U(D) or as the endpoint of an upright 
diagonal over F 1 . One common interpretation of 
sc is vor 0> F(^(-D)), the truncated Voronoi function. 
When this is the interpretation, we introduce fur- 
ther variables: 

quo[w, s(v, Fx)] Vv e F\, 
qao[s{v,F\),v] Vv e Fx, 
Adih[w,Fi] ' VveFx, 

We interpret the variables as follows. If w = s(v), 
and the triangle (0,t/,w J ) has circumradius rj at 
most to, then 

quo^w] 7 = quo( J R(|i; / |/2,77,i )), 
quolw^} 1 = quo(R(\w I \/2,n,t )). 

If the circumradius is greater than to, we take 

quo^jw] 7 = quo[w,w] 7 = 0. 

The variable Adih has the following interpretation: 

I otherwise. 

Under these interpretations, the following identity 
is satisfied: 

sc[Ji] - sol[Fi]0o + £„ eFl Adih K F i] 
-4:S ct^quo[v,w). 

This relation and the constants that appear in it are 
based on [5, Appendix A. 4] . The final sum runs over 
all pairs (v, w), where v = s(w, F\) or w — s(v, Fi). 

For this to be useful, we need good inequalities 
governing the individual variables. Such inequal- 
ities for Adih^jf] and quo[w,w] are found in [5, 
A. 4]. To make of these inequalities, it is necessary 
to have lower and upper bounds on a[v, F] and y[v]. 
We obtain such bounds as LP-dcrivcd inequalities 
in the sense of Remark 10.3. 



11.4 Triangle and Quad Branching 

The paper [5, Appendix A. 7] divides triangular 
faces into two cases: 

+ y[v 2 ] + y[v 3 ] < 6.25. 
y[vi] + y[v2] + y[v3] > 6.25. 

whenever sufficiently good bounds are not obtained 
as a single linear program. That paper also divides 
quadrilateral faces into four cases: two flat quar- 
ters, two flat quarters with diagonal running in the 
other direction, four upright quarters forming an 
octahedron, no quarters. In general, if there are r\ 
triangles and T2 quadrilaterals, we obtain as many 
as 2 ri+2r2 cases by breaking the various triangles 
and quadrilaterals into subcases. 

We break triangular faces and quadrilaterals into 
subcases as needed in the linear programs that fol- 
low without further comment. 



12 Elimination of Aggregates 

The proof of the following theorem occupies the en- 
tire section. It eliminates all the pathological cases 
that we have had to carry along until now. The 
terminology in the next lemma comes from [6] (an- 
chors, S3 , S3 , and so forth). This terminology is 
reviewed in Appendix 16. 

Theorem 12.1. Let D be a contravening decompo- 
sition star, and let G be its tame graph. Every face 
of G corresponds to exactly one standard region of 
D. No standard region of D has any enclosed ver- 
tices from U{D). (That is, enclosed vertices have 
height at least 2t . The decomposition star D does 
not contain any S^ or configurations. If an up- 
right diagonal with at least 5 anchors appears, the 
5 anchors are the five corners of U(D) lying over a 
pentagonal standard region. 

12.1 A pentagonal hull with n = 8 

The next few sections treat the nonpolygonal stan- 
dard regions illustrated in Remark 4.2. In this sub- 
section, there is an aggregate of the octagonal region 
and a triangle has a pentagonal hull. Let P denote 
this aggregate. 

Lemma 12.2. Let G be a contravening plane graph 
with this aggregate. Some vertex on the pentagonal 
face has type not equal to (3, 0, 1). 
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Proof. If every vertex on the pentagonal face has 
type (3,0,1), then at the vertex of the pentagon 
meeting the aggregated triangle, the four triangles 
together with the octagon give 

*8 + ]T tlp(4, 0, 2tt - 2(1.153)) > (4ttC - 8) pt, 

(4) 

so that the graph does not contravene. □ 

For a general contravening plane graph with this 
aggregate, we have bounds 

o-f(D) <pt+s 8 , 
T F (D) > t 8 . 

We add the inequalities rsc[F] > t$ and sc[F] < 
pt + sg to the exceptional face. There is no other 
exceptional face, because t$ +t 5 > (4n( — 8) pt. We 
run the linear programs for all tame graphs with 
the property asserted by Lemma 12.2. Every upper 
bound is less than 8 pt, so that there are no contra- 
vening decomposition stars with this configuration. 

12.2 n — 8, hexagonal hull 

We treat the two cases from Remark 4.2 that have 
a hexagonal hull (Figure 8). One can be described 
as a hexagonal region with an enclosed vertex that 
has height at most 2t and distance at least 2to 
from each corner over the hexagon. The other is 
described hexagonal region with an enclosed 
vertex of height at most 2t<j, but this time with 
distance less than 2i from one of the corners over 
the hexagon. 

The argument for the case n = 8 with hexago- 
nal hull is similar to the argument of Section 12.1. 
Add the inequalities t(R) > t 8 and a(R) < s$ for 
each hexagonal region. Run the linear programs for 
all tame graphs, and check that these additional in- 
equalities yield linear programming bounds under 
8pt. 

12.3 n = 7, pentagonal hull 

We treat the two cases illustrated in Figure 8 that 
have a pentagonal hull. One can be described as a 
pentagon with an enclosed vertex that has height 
at most 2to and distance at least 2to from each cor- 
ner of the pentagon. The other is described as a 
pentagon with an enclosed vertex of height at most 
2to, but this time with distance less than 2to from 
one of the corners of the pentagon. 



In discussing various maps, we let Vi be the 
corners of the regions, and we set yt = \vi\ and 
■tjij = \vi — Vj\. The subscript F is dropped, when 
there is no great danger of ambiguity. 

Add the inequalities t[F] > tr, a[F] < sr for the 
pentagonal face. There is no other exceptional re- 
gion, because t 5 + tr > (47r£ — 8) pt. With these 
changes, of all the tame plane graphs with a pen- 
tagonal face, all but one of the linear programs give 
a bound under 8 pt. 

The plane graph G15 that remains is easy to de- 
scribe. It is the plane graph with 11 vertices, ob- 
tained by removing from an icosahedron a vertex 
and all the edges that meet at that vertex. 

We treat the case G15. Familiarity with Ap- 
pendix 16 will be assumed in the following pas- 
sage. Let V12 be the vertex enclosed over the pen- 
tagon. We let Vi,... ,t>5 be the five corners of 
U(D) over the pentagon. Break the pentagon into 
five simplices along (0,^12): Si — (0, V12, Vi, Vj+i). 
We have LP-dcrived bounds (in the sense of Re- 
mark 10.3) y[vi] < 2.168, and a[v u F] < 2.89, for 
2 = 1,2, 3, 4, 5. In particular, the pentagonal region 
is convex. 

Another LP-derived inequality is 

a[F] > -0.2345. 

Since —0.4339 is less than this the lower bound, the 
configuration does not occur. Similarly, since 
—0.25 is less than the lower bound, S4 does not 
occur ([6, 3.7 and 3.8]. 

Suppose that there is a loop in context (n, k) = 
(4,2). Again the score is less than the LP-lower 
bound, showing that this does not occur: 

a R (D) < s 7 + z loop (A, 2) < -0.2345. 

The constants come from Table [6, Table 5.11] and 
[6, Theorem 4.4]. 

We have the LP-derived inequality 

t[F] < 0.644. 

If we branch and bound on the triangular faces, this 
LP-derived inequality can be improved to 

t[F] < 0.6079. 

If there is a loop other than (4,2) and (4,1), the 
linear program becomes infeasiblc: 

t[F] < 0.644 < t 7 + S loop (n, k) < t[F}. 
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We conclude that all loops have context (n, k) = 
(4,1). 

Case 1. The vertex v\ 2 has distance at least 2to 
from the five corners ofU(D) over the pentagon. 

The interval calculations relevant in Case 1 ap- 
pear in Appendix 17.2. 

The penalty to switch the pentagon to a pure voro 
score is at most 5£r (see [6, Aio]). There cannot be 
two flat quarters because then 

|«i 2 | >£(S(2,2,2,2fo,2\/2,2\/2),2*o,2*o,2*o) > 2t . 

(Case 1-a) Suppose there is one flat quarter, 
\vi — v 4 \ < 2\[2. There is a lower bound of 1.2 on 
the dihedral angles of the simplices (0, t>i 2 , v%, v%+i). 
This is obtained as follows. The proof relies on the 
convexity of the quadrilateral region. We leave it 
to the reader to verify that the following pivots can 
be made to preserve convexity. Disregard all ver- 
tices except Vi,V2,V3,V4, V\2- We give the argument 
that dih(0, V\2, Vi, Vi) > 1.2. The others are similar. 
Disregard the length \v\ — Vi\. Wc show that 

sd := dih(0,Ui2,Vi,U2) + dih(0, v 12 , v 2 , v 3 ) 
+ dih(0, v 12 , v 3 , Vi) <2ir- 1.2. 

Lift V12 so I f 12 ] = 2*o- Maximize sd by taking \vi — 
v 2\ = \ v 2 — v 3\ = \ v 3 — v i\ = 2f . Fixing w 3 and 
Vi, pivot i>i around (0,^12) toward U4, dragging t> 2 
toward v\2 until |w 2 — U12I = 2*o. Similarly, we 
obtain \v 3 — v\ 2 \ = 2t . We now have sd < 3(1.63) < 
2ir — 1.2, by an interval calculation. 

Return to the original figure and move 
v\2 without increasing |ui 2 | until each simplex 
(0,v\2,Vi,Vi + i) has an edge (v\2,vj) of length 2* . 
Interval calculations show that the four simplices 
around V12 squander 

2tt(0.2529) - 3(0.1376) - 0.12 > (4ttC - 8) pt+ 5fr- 

(Case 1-b) Assume there are no flat quarters. An 
LP-derived bound on the perimeter \ v i ~ v i+i\ 1S 
1.407. We have arc(2, 2, x)" = 2x/(16 - x 2 ) 3 / 2 > 0. 
The arclength of the perimeter is therefore at most 

2 arc(2, 2, 2t ) + 2 arc(2, 2, 2) + arc(2, 2, 2.387) < 2tt. 

There is a well-defined interior of the spherical pen- 
tagon, a component of area < 2tt. If we deform 
by decreasing the perimeter, the component of area 
< 2n does not get swapped with the other compo- 
nent. 

Disregard all vertices but V\, . . . ,v 5 , w 12 . If a ver- 
tex Vi satisfies \vi — t>i 2 | > 2to, deform Vi as in [6, 



Section 4.9] until — Vi\ = \vi — = 2, or 

\vi — Vi2\ = 2to- If at any time, four of the edges 
realize the bound \vi — = 2, we have reached 
an impossible situation, because it leads to the con- 
tradiction (see Appendix 17.2) 

(5) 

2tt = ^ dih < 1-51 + 4(1.16). 

(This inequality relies on the observation, which we 
leave to the reader, that in any such assembly, piv- 
ots can by applied to bring \v\2 — Vi\ = 2t for at 
least one edge of each of the five simplices.) 

The vertex u 12 may be moved without increasing 
|i>i 2 | so that eventually by these deformations (and 
reindcxing if necessary) we have \vu — Vi\ = 2t , 
2 = 1,3, 4. (If we have % = 1,2, 3, the two dihedral 
angles along (0,w 2 ) satisfy < 2(1.51) < n, so the 
deformations can continue.) 

There are two cases. In both cases \vi— v\ 2 \ = 2t , 
fori =1,3,4. 

(i) \vi2 - v 2 \ = \vi2 - v 5 \ = 2t , 

(ii) \vi2 - v 2 \ = 2i , \v 4 - v 5 \ = \v 5 - vi\ = 2, 

Case (i) follows from interval calculations 

T o > 2tt(0.2529) - 5(0.1453) > 0.644 + 7£r- 

In case (ii), we have again 

2tt(0.2529) - 5(0.1453). 

In this interval calculation we have assumed that 
1^12 — ^5 1 < 3.488. Otherwise, setting S = 
(vi2,Vi,V5,vi), we have 

A(S) < A(3.488 2 , 4, 4, 8, (2* ) 2 , (2* ) 2 ) < 0, 

and the simplex does not exist. (\vi — Vi\ > 2\/2 
because there are no flat quarters.) This completes 
Case 1. 

Case 2. The vertex v\2 has distance at most 2* 
from the vertex v\ and distance at least 2* from the 
others. 

Let (0,i>i3) be the upright diagonal of a loop 
(4, 1). The vertices of the loop are not V2, V3, Vi, v$ 
with V\2 enclosed over (0,^2,^5,^13) by [6, Lemma 
3.6]. The vertices of the loop are not v 2 , t> 3 , t>4 , v 5 
with v\2 enclosed over (0, v\, V2, v§) because this 
would lead to a contradiction 

Vi2 > £{S{2, 2, 2, 2*o, 2*o, 3.2), 2* , 2* , 2) > 2* , 
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or 

J/12 > £(S(2, 2, 2, 2t , 2t , 3.2), 2, 2i , 2) > 2t . 

We get a contradiction for the same reasons unless 
is an edge of some upright quarter of every 
loop of type (4, 1). 

We consider two cases. (2- a) There is a flat quar- 
ter along an edge other than (vi,vi2). That is, the 
central vertex is V2, V3, V4, or D5. (Recall that the 
central vertex of a flat quarter is the vertex other 
than the origin that is not an cndpoint of the di- 
agonal.) (2-b) Every flat quarter has central vertex 
vi- 

Case 2-a. We erase all upright quarters including 
those in loops, taking penalties as required. There 
cannot be two flat quarters by geometric consider- 
ations 

£(S(2, 2, 2, 2\/2, 2V2, 2i ), 2t , 2t , 2) > 2t Q 
£(S(2, 2, 2, 2V2, 2V2, 2t ), 2, 2t , 2t ) > 2i 

The penalty is at most 7£r- We show that the 
region (with upright quarters erased) squanders 
> 7£r + 0.644. We assume that the central ver- 
tex is v 2 (case 2-a-i) or w 3 (case 2-a-ii). In case 
2-a-i, we have three types of simplices around V12, 
characterized by the bounds on their edge lengths. 
Let (0,1)12,1)1,1)5) have type A, (0, V12, ^5,^4) and 
(0,^12,^4,^3) have type B, and let (0, V12, V3, v\) 
have type C. In case 2-a-ii there are also three 
types. Let (0, v\ 2 , vi, 1*2) and (0, Vi2,v\,vs) have 
type A, (0,v\2,V5,V4) type B, and (0,^12,^2,^4) 
type D. (There is no relation here between these 
types and the types of simplices A, B, C defined 
in Appendix 16.2 and 16.3.) Upper bounds on the 
dihedral angles along the edge (0, V12) are given in 
Appendix 17.2. These upper bounds come as a re- 
sult of a pivot argument similar to that establishing 
the bound 1.2 in Case 1-a. 

These upper bounds imply the following lower 
bounds. In case 2-a-i, 

dih > 1.33 (A), 
dih > 1.21 (B), 
dih > 1.63 (C), 

and in case 2-a-ii, 

dih > 1.37 (A), 
dih > 1.25 (B), 
dih > 1.51 (£>), 

In every case the dihedral angle is at least 1.21. 
In case 2-a-i, the interval calculations of Appendix 



17.2 lead to a lower bound on what is squandered by 
the four simplices around (0,^12). Again, we move 
V12 without decreasing the score until each simplex 
(0,^12, Vi,Vi+i) has an edge satisfying \v\ 2 — Vj\ < 
2t . Interval calculations give 

E (4) tq > 2tt(0.2529) - 0.2391 - 2(0.1376) - 0.266 

> 0.808. 

In case 2-a-ii, we have 

£ (4) t > 2tt(0.2529) - 2(0.2391) - 0.1376 - 0.12 

> 0.853. 

So we squander more than 7£r + 0.644, as claimed. 

Case 2-b. We now assume that there are no 
flat quarters with central vertex v 2 ,. ■■ , V5. We 
claim that v\ 2 is not enclosed over (0 7 vi,v 2 ,v 3 ) 
or (0, Vi, V5, V4). In fact, if v\ 2 is enclosed over 
(0,^1,^2,^3), then we reach the contradiction 

7r < dih(0,Ui2,t>i,U2) + dih(0, v 12 , v 2 , v 3 ) 
< 1.63 + 1.51. 

We claim that V12 is not enclosed over 
(0, V5,vi, V2)- Let Si = (0,ui2,ui,U2), and S2 = 
(Q,v\2,vi, V5). We have the linear programming 
bound 

2/4 (Si) + y4(S 2 ) = \vi - v 2 \ + \vi - v 5 \ < 4.804. 
An interval calculation gives (Appendix 17.3) 
£ (2 ) dih(^) < E (2 ) (dMS l ) + 0.5(0.4804/2 - y 4 (S 

< IT. 

So V12 is not enclosed over (0, vi,v 2 ,v 5 ). 

Erase all upright quarters, taking penalties as re- 
quired. Replace all flat quarters with voro-scoring 
taking penalties as required. (Any flat quarter has 
v\ as its central vertex.) We move V12 keeping \v\ 2 \ 
fixed and not decreasing \vi 2 — vi\. The only ef- 
fect this has on the score comes through the quoins 
along (0,vi,vi 2 )- Stretching \vn — Vi\ shrinks the 
quoins and increases the score. (The sign of the 
derivative of the quoin with respect to the top edge 
is computed in the proof of [6, Lemma 4.9].) 

If we stretch \v\ 2 — vi\ to length 2to, we are done 
by case 1 and case 2-a. (If deformations produce a 
flat quarter, use case 2-a, otherwise use case 1.) By 
the claims, we can eventually arrange (reindexing if 
necessary) so that 

(i) \v 12 - v 3 \ = \v 12 - v 4 \ = 2t , or 

(ii) \vi2 -v 3 \ = \vi2 -v 5 \ = 2t . 
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We combine this with the deformations of [6, Sec- 
tion 4] so that in case (i) we may also assume that 
if \v 5 — v\ 2 \ > 2to, then \v4 — v$\ = \v 5 — v\\ — 2 and 
that if \v 2 — vi 2 \ > 2t , then \vi — v 2 \ — \v 2 — v 3 \ = 2. 
In case (ii) we may also assume that if \v4 — v\ 2 \ > 
2to, then \v 3 — V4\ — \v4 — v 5 \ — 2 and that if 
\v 2 — vi 2 \ > 2t , then \vi — v 2 \ = \v 2 — W3I = 2. 

Break the pentagon into subregions by cutting 
along the edges (vi 2 ,Vi) that satisfy \vi 2 — Vi\ < 2i - 
So for example in case (i), we cut along (^12,^3), 
(vi 2l v 4 ), (v 12 ,vi), and possibly along (v 12 ,v 2 ) and 
(vi 2 ,v 5 ). This breaks the pentagon into triangular 
and quadrilateral regions. 

In case (ii), if |v 4 — v i2 \ > 2t , then the argument 
used in Case 1 to show that |v 4 — u 12 | < 3.488 applies 
here as well. In case (i) or (ii), if \v i2 — v 2 \ > 2i , 
then for similar reasons, we may assume 

A(|«i2 - « 2 | 2 ,4,4,8, (2i ) 2 , \v 12 - Vl \ 2 ) > 0. 

This justifies the hypotheses for the interval calcu- 
lations in Appendix 17.2. We conclude that 

^r > 2tt(0.2529)-3(0.1453)-2(0.2391) > 0.6749. 

Recall that the LP-upper bound on what is squan- 
dered is 0.6079. If the penalty is less than 0.067 = 
0.6749 - 0.6079, we are done. 

We have ruled out the existence of all loops except 
(4, 1). Note that a flat quarter with central vertex 

gives penalty at most 0.02 by [6, Lemma 3.11.3]. 
If there is at most one such a flat quarter and at 
most one loop, we are done: 

3£r + 0.02 < 0.067. 

Assume there are two loops of context (n, k) = 
(4, 1). They both lie along the edge (v\, v\ 2 ), which 
precludes any unmasked flat quarters. If one of 
the upright diagonals has height > 2.696, then the 
penalty is at most 3£r + 3£y < 0.067. Assume both 
heights are at most 2.696. The total internal angle 
of the exceptional face at v\ is at least four times 
the dihedral angle of one of the flat quarters along 
(0, vi), or 4(0.74) by an interval calculation. This 
is contrary to the linear programming upper bound 
(0.289) on the dihedral angle along (0,t>i). This 
completes Case 2. This shows that heptagons with 
pentagonal hulls do not occur. 

12.4 S3, S+ 

We use the bounds on the score and on what is 
squandered from [6, Sections 3.7 and 3.8]: S 4 (a < 



-0.25, t > -0.4), S3 (er < -0.4339, r > 0.5606). 
If we have the configuration , there is only one 
exceptional region (0.5606 + t 5 > (47r( — 8)pt). The 
configuration with six standard regions around a 
vertex from Section 4.1 does not occur because a 
calculation with tlp (see Section 6.1) show that the 
five quasi-regular tetrahedra in the configuration 
squander > 6 pt, giving 6 pi +0.5606 > (4:Tr( — 8) pt. 

We add the inequalities t > 0.5606 and a < 
—0.4339 to the exceptional region. All linear 
programming bounds drop under 8 pt when these 
changes are made. 

The configuration requires more work. We 
begin with a lemma. 

Lemma 12.3. Let a be the dihedral angle along the 
large gap in an S4 configuration. Let vi and v 2 be 
the anchors of U(D) lying along the large gap. If 
\vi\ + \v 2 \ < 4.6, then a > 1.78 and the score of the 
four quarters is at most —0.31547. 



Proof. The bound a > 1.78 is an interval calcu- 
lation (Appendix 17.4). The upper bound on the 
score is a linear programming calculation involving 
the inequality a > 1.78 and the known inequalities 
on the score of an upright quarter. □ 



Add the inequalities a < —0.25, r > 0.4 at the 
exceptional regions. The configuration does not 
appear in a pentagon. Run the linear programs for 
all tame plane graphs with an exceptional region 
that is not a pentagon. If this linear program fails 
to produce a bound of 8pt, we use the lemma to 
branch into two cases: either y[v\] + y[v 2 ] > 4.6, or 
a[R] < —0.31547. In every case the bound drops 
below 8 pt. 



12.5 Type (n, k) = (5, 1) 

We return briefly to the case of six standard re- 
gions around a vertex discussed in Section 4.1. In 
the plane graph they are aggregated into an octahe- 
dron. We take each of the remaining cases with an 
octahedron, and replace the octahedron with a pen- 
tagon and six triangles around a new vertex. There 
are eight ways of doing this. All eight ways in each 
of the cases gives an LP bound under 8pt. This 
completes this case. 
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12.6 Four anchors 



13.1 Pent and Hex Branching 



In the context (n,k) — (4,3), the standard region 
R must have at least seven sides n(R) > 7. Then 

t{D) > £ 7 + ^ oop (4,3) 
> (4ttC - 8) pt. 

Thus, we may assume that this context does not 
occur. 



12.7 Five Anchors 

Now turn to the decomposition stars with an up- 
right diagonal with five anchors. Five quarters 
around a common upright diagonal in a pentagonal 
region can certainly occur. We claim that any other 
upright diagonal with five anchors leads to a decom- 
position star that does not contravene. In fact, the 
only other possible context is (n, k) — (5, 1) (see [6, 
5.11]). 

If this appears in an octagon, we have 

t(D) > 5 loop (5, l)+t s > (4ttC - 8) pt. 

If this appears in a heptagon, we have 

t(D) > Si oop (5, l) + t 7 + 0.55 pt > (4ttC - 8) pt, 

because there must be a vertex that is not a corner 
of the heptagon. It cannot appear on a pentagon. 
If it appears on a hexagonal region R, we have 

t(D) > Si oop (5,l)+t 6 

and <j(D) < zi oop (5, 1) + sq. Add these constraints 
to the linear program of the tame graphs with a 
hexagonal face. The LP-bound on a(D) with these 
changes is < 8pt. 



13 Branching and Bound 

In Section 11.4 we mentioned that when the linear 
programs do not give sufficiently good bounds, we 
branch into various subcases depending on condi- 
tions triangles and quadrilaterals in the graph. 

This section describes additional branching for 
exceptional regions. We assume the results from 
the Section 12 that eliminate the most unpleasant 
types of configurations. 



In Section 11.2 we discussed face refinements in ab- 
stract terms. We now describe these refinements in 
detail. 

The possibilities are listed in the diagram only up 
to symmetry by the dihedral group action on the 
polygon. We do not prove the completeness of the 
list, but its completeness can be seen by inspection, 
in view of the comments that follow here and in 
Section 13.4. 




Figure 13: Pentagonal Face Refinements 




Figure 14: Hexagonal Face Refinements. The only 
figures with a penalty are the first two on the top 
row and those on the bottom row. The first two 
on the top row have penalties 2(0.008) and 0.008. 
Those on the bottom row have penalties 3£r, 3£r, 
£r + 2£y, and + 2£ v . 

The conventions for generating the possibilities 
are different for the pentagons and hexagons than 
for the heptagons and octagons. We describe the 
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pentagons and hexagons first. We erase S3. If there 
is one loop we leave the loop in the figure. If there 
are two loops (so that both necessarily have context 
(n, k) = (4,1)), we erase one and keep the other. 
We draw all edges from an enclosed upright diago- 
nal to its anchors, and all edges of length [2to, 2y/2] 
that are not masked by upright quarters. Since the 
only remaining upright quarters belong to loops, the 
four simplices around a loop are anchored simplices 
and the edge opposite the diagonal has length at 
most 3.2. 

13.2 Branching on pentagonal faces 

The first three inequalities of Appendix 17.12 and 
the first inequality of Appendix 17.29 have been de- 
signed with pentagons in mind. The additional in- 
equalities in Appendix 17.30, 17.31, 17.32, and the 
last two inequalities of 17.12 have been designed for 
subregions in hexagonal regions. 

13.3 Hept and Oct Branching 

When the figure is a heptagon or octagon, we pro- 
ceed differently. We erase all Sjj" configurations and 
all loops (either context (n, k) = (4, 1) or (4, 2)) and 
draw only the flat quarters. An undrawn diagonal of 
the polygon has length at least 2to. Thus, in these 
cases much less internal structure is represented. 

In the cases where Sjj~ or loops have been erased, 
a number indicating a penalty accompanies the di- 
agram. These penalties are derived in Section 13.4. 

Here are some special arguments that are used 
for heptagons and octagons. 

Suppose that the face refinement gives two sub- 
regions: a flat quarter and one with parameters 
(n, k) = (n — 1, 1), with n = 6 or 7. We have the 
inequality: 

a R {D) < a + Z(n - 1, 1) + £r + 2£ y . 

The penalty term £r+2£y comes from a possible an- 
chored simplex that masks the flat quarter. Let v be 
the central vertex of the flat quarter. Let (i>i, V2) be 
its diagonal. Masked flat quarters satisfy restrictive 
edge constraints. Thus, we have one of the following 
three possibilities: 

1. y[v] > 2.2, 

2. y[v 1 ,v 2 }>2.7, 

3. (j R {D) < & + Z(n- 1,1). 




Figure 15: Hept Face Refinements. The pair (a, b) 
represents the penalty et£r + b£v- 

We proceed similarly, if the face refinement has 
three subregions: two flat quarters R\, R2 and one 
(R3) with parameters (n,k) = (n — 2,2). We have 
the cases: 

1. The height of a central vertex is at least 2.2. 

2. The diagonal of a flat quarter is at least 2.7. 
3. 

a R {D) <a 1 +a 2 + Z{n- 2,2) 
t r (D) >T 1 +f 2 + D(n-2,2). 

With heptagons, it is helpful on occasion to use 
an upper bound on the penalty of 3£r = 0.04683. 
This bound holds if neither flat quarter is masked 
by a loop. For this, it suffices to show that first two 
of the given three cases do not hold. 

13.3.1 Octagon Branching 



If there is a loop of context (11, k) — (4,2), we 
have, by Lemma 16.4, the upper bound 

Sri 1 Zloop 

(4,2) 
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Figure 16: Oct Face Refinements. The penalties 
reading from left to right are 6£r and 6£r on the 
top row, 4£ r + 2£y, 4£ r + 2£y and 4£ r + 2£y on the 
middle row, and 2£r + 4£y , 2£r + 4£y, and on the 
bottom row. 

on the score of exceptional standard regions R, with 
n = n(R). 

If there is no loop of context (n, k) = (4, 2), and 
if there are two flat quarters in the face refinement, 
we have the upper bound 

(Ti + (T 2 + D(n -2,2) + 2(£r + 2£y). 
13.4 Penalties 

The reader should review Appendix 16 before read- 
ing this section. 

Erasing an upright quarter of compression type 
gives a penalty of at most £r and one of Voronoi 
type gives at most £y. We take the worst possible 
penalty. It is at most n£r in an n-gon. If there is 
a masked flat quarter, the penalty is at most 2£y 
from the two upright quarters along the flat quarter. 
We note in this connection that both edges of a 
polygon along a flat quarter lie on upright quarters, 
or neither does. 

If an upright diagonal appears enclosed over a 
flat quarter, the flat quarter is part of a loop with 
context (n, k) = (4, 1), for a penalty at most 2£p + 
£y. This is smaller than the bound on the penalty 
obtained from a loop with context (n,k) = (4,1), 
when the upright diagonal is not enclosed over the 



flat quarter: 

£r + 2£y. 

So we calculate the worst-case penalties under the 
assumption that the upright diagonals are not en- 
closed over flat quarters. 

A loop of context (n, k) = (4, 1) gives £r + 2£y or 
3£r- A loop of context (n, k) = (4,2) gives 2£r or 
2£y. 

If we erase S^", there is a penalty of 0.008 (or if 
it masks a flat quarter.) This is dominated by the 
penalty 3£r of context (n, k) = (4, 1). 

Suppose we have an octagonal standard region. 
We claim that a loop does not occur in context 
(n, k) = (4,2). If there are at most three vertices 
that are not corners of the octagon, then there are 
at most 12 quasi-regular tetrahedra, and the score 
is at most 

s s + 12 pt< 8pt. 

Assume there are more than three vertices that are 
not corners over the octagon. We squander 

ts + 5 loop (A, 2) + 4r LP (5, 0) > (4^C - 8) pt. 

As a consequence, context (n, k) = (4, 2) does not 
occur. 

So there are at most 2 upright diagonals and at 
most 6 quarters, and the penalty is at most 6£r- 
Let / be the number of flat quarters This leads to 

f6Cr, / = o,i, 

/ = 2, 

TTp = < 

|2£r + 4^, / = 3, 
U / = 4. 

The is justified by a parity argument. Each up- 
right quarter occurs in a pair at each masked flat 
quarter. But there is an odd number of quarters 
along the upright diagonal, so no penalty at all can 
occur. 

Suppose we have a heptagonal standard region. 
Three loops are a geometric impossibility. Assume 
there are at most two upright diagonals. If there is 
no context (n, k) — (4, 2), then we have the follow- 
ing bounds on the penalty 

'6£r, / = 0, 

4£ r + 26/, / = 1, 

I 3£r, / = 2, 

& + 2£ v , / = 3. 
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If an upright diagonal has context (n,k) = (4,2), 
then 

f5£ r , / = 0,1, 

*f = < 3£r + 2£ y , / = 2, 
[Cr + ^y, / = 3. 

This gives the bounds used in the diagrams of cases. 



13.5 Branching on Upright Diago- 
nals 

We divide the upright simplices into two domains 
depending on the height of the upright diagonal, 
using yi = 2.696 as the break point. We break the 
upright quarters into cases: 

1. The upright diagonal has height at most 2.696. 

2. The upright diagonal (0, v) has height at least 
2.696, and some anchor w along the flat quarter 
satisfies \w\ > 2.45 and \v - w\ > 2.45. (There 
is a separate case here for each anchor w.) 

3. The upright diagonal (0, v) has height at least 
2.696, and every anchor w along the Hat quar- 
ter satisfies |io| < 2.45 and \v — w\ < 2.45. 

Many inequalities have been specially designed to 
hold on these smaller domains. They are included 
into the linear programming problems as appropri- 
ate. 

When all the upright quarters can be erased, then 
the case for upright quarters follows from some 
other case without the upright quarters. An up- 
right quarter can be erased in the following situa- 
tions. If the upright quarter is compression type 
and > 2.696, then 

i>r < vor 

([6, A 10 ]). (If there are masked flat quarters, they 
become scored by a.) If upright quarters is of 
Voronoi type and the anchors w satisfy \w\ < 2.45 
and \v — io | < 2.45, then the quarter can be erased: 

v < vor 

In general, we only have the weaker inequality ([6, 
An]) 

v < vor Hv, = 0.003521. 



13.6 Branching on Upright Quarters 

In a pentagon or hexagon, consider an upright di- 
agonal with three upright quarters, that is, context 
(n,k) = (4,1). If the upright diagonal is at most 
2.696, and if an upright quarter shares both faces 
along the upright diagonal with other upright quar- 
ters, then we may assume that the upright quarter 
has compression type. For otherwise, there is a face 
of circumradius at least y/2, and hence two upright 
quarters of Voronoi type. The inequality 

octavor < octavor —0.008, (5) 

if yi £ [2to, 2.696], and 77126 > \/2 shows that the 
upright quarters can be erased without penalty be- 
cause £r < 2(0.008). If erased, the case is treated 
as part of a different case. 

This allows the inequalities for i>\> in Appendix 

17.6 and 17.7 to be used. Furthermore, it can often 
be concluded that all three upright quarters have 
compression type. For this, we use the Inequalities 
17.27 and 17.23, which can often be used to show 
that if the anchored simplex has a face of circum- 
radius at least \/~2, then the linear programming 
bound on the score is less than 8 pt. 

13.7 Branching on Flat Quarters 

There are a few other interval-based inequalities 
that are used in particular cases. Before an in- 
equality that uses compression scoring is used on a 
flat quarter, it is necessary to verify that the quar- 
ter satisfies y\ < 2.2, y 4 < 2.7,77234,77456 < \/2, to 
insure that the flat quarter has compression type. 
The circumradius is not a linear-programming vari- 
able, so its upper bound must be deduced from 
edge-length information. 

Information about the internal structure of an ex- 
ceptional face gives improvements to the constants 
lApt and 1.5 pt of Property 4 in the definition of 
admissible weight assignments. (The bounds re- 
main fixed at lApt and 1.5 pt, but these arguments 
allow us to specify more precisely what simplices 
contribute to these bounds.) These constants con- 
tribute to the bound on the t(D) through the ad- 
missible weight assignment. Assume that at the 
vertex v there are four tetrahedra and an excep- 
tional face, and that the exceptional face has a flat 
quarter with central vertex v. The calculations of 
Section 7.4 show that the four quarters and excep- 
tional squander at least 1.5 pt. If there is no flat 
quarter (masked or unmasked) whose first edge lies 
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along (0,v), then the four quasi-regular tetrahcdra 
at v squander at least 1.5 pt We can make similar 
improvements when tri(v) — 3. 

If we have LP-derived inequalities on a flat quar- 
ter whose corners Vi have height at most 2.14, and 
if the diagonal has length less than 2.77, then the 
circumradius of the face containing the origin and 
diagonal is at most 77(2.14,2.14,2.77) < y/2. This 
allows us to combine the cases defining a into three 
cases. 

1. The simplex has compression type. 

2. The diagonal has length < 2.7 and 77456 > \/2. 

3. The diagonal has length > 2.7. 

In the last two cases, the scoring is by vor . We use 
the specially tailored Inequalities 17.15 and 17.16. 

13.8 Branching on Type A simplices 

Suppose that the vertices of a triangular subregion 
come by projection from a simplex with vertices at 
the origin and in U(D). Suppose that the simplex 
has two edges of length in [2i ,2\/2], say j/ 5 , j/ 6 <G 
[2£o,2\/2], 7/4 <E [2, 2to]. These may be simplices 
of type A. Suppose that one edge is at most 2.77. 
Such simplices fall into the following domains. 

1. The simplex has type A, so that 7/5,7/5 g 
[2to,2.77], and the scoring is vor, the analytic 
Voronoi function. 

2. An edge (say t/ 6 ) of the face shared with the 
flat quarter has length t/ 6 > 2.77. 

3. The edges have lengths 7/5,7/g g [2i ,2.77], and 
7/456 > \/2. 

In the last two cases, the scoring is by the truncated 
Voronoi function voro . We use the specially tailored 
Inequalities 17.17. 

If there is an enclosed vertex of height at most 
y/2, then we can use Inequalities 17.18. 

13.9 Branching on Quadrilateral 
subregions 

The inequalities of 17.27 hold for a quadrilateral 
subregion, if certain conditions are satisfied. One 
of the conditions is 7/4 € [2\/2, 3.0], where 7/4 is a 
diagonal of the subregion. Since this diagonal is 
not one of the linear programming variables, these 



bounds cannot be verified directly from the linear 
program. Instead we use the interval calculation 
17.13 which relates the desired bound 7/4 < 3 to the 
linear programming variables a[v, F], t/ 2 , 7/3, 7/5, and 

ye- 

13.10 Conclusion 

By combinations of branching along the lines set 
forth in the preceding sections, a sequence of linear 
programs is obtained that establishes that t(D) is 
less than 8 pt. For details of particular cases, the in- 
terested reader can consult the log files in [8] , which 
record which branches are followed for any given 
tame graph. (For most tame graphs, a single linear 
program suffices.) 

This completes the proof of the Kepler conjec- 
ture. 
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14 Appendix. Colored Spaces 

It is claimed at various places that the optimiza- 
tion problem is that of a continuous function on a 
compact topological space. That is, the function a 
on the space X of decomposition stars is continu- 
ous. The purpose of this appendix is to describe 
the sense in which this is true. 

We begin with an example that should make 
things clear. Suppose that we have a discontinu- 
ous piecewise linear function on the unit interval 
[—1,1], as in Figure 17. It is continuous, except at 
a; = 0. 




Figure 17: A piecewise linear function 

We break the interval in two at x = 0, forming 
two compact intervals [—1,0] and [0,1]. We have 
a continuous functions /_ : [—1,0] — > M and / + : 
[0,1], such that 



1 f+(x) otherwise. 

We have replaced the discontinuous function by a 
pair of continuous functions on smaller intervals, at 
the expense of duplicating the point of discontinuity 
x = 0. We view this pair of functions as a single 
function F on the compact topological space with 
two components 

[-1,0] x {-} and [0,1] x {+}. 

where F(x, a) = f a (x), and a € { — , +}. 

This is the approach that we follow in general 
with the Kepler conjecture. The function a is de- 
fined by a series of case statements, and the function 
does not extend continuously across the boundary 
of the cases. However, in the degenerate cases that 
land precisely between two or more cases, we form 



multiple copies of decomposition star for each case, 
and place each case into a separate compact domain 
on which the function a is continuous. 

This can be formalized as a colored space. A col- 
ored space is a topological space X together with 
an equivalence relation on X with the property that 
no point x is equivalent to any other point in the 
same connected component as x. We refer to the 
connected components as colors, and call the points 
of X colored points. We call the set of equivalence 
classes of X the underlying uncolored space of X. 
Two colored points are equal as uncolored points if 
they are equivalent under the equivalence relation. 

In our example, there are two colors "— " and 
"+." The equivalence class of (x, a) is the set of 
pairs (x,b) with the same first coordinate. Thus, if 
x =/= 0, the equivalence class contains one element 
(x, sign(x)), and in the boundary case x — there 
are two equivalent elements (0,—) and (0,+). 

In our treatment of the Kepler Conjecture, there 
are various cases: whether an edge has length less 
than or greater than 2to, less than or greater than 
whether a face has circumradius less than or 
greater than V2, and so forth. By duplicating the 
degenerate cases (say an edge of exact length 2to), 
creating a separate connected component for each 
case, and expressing the optimization problem on a 
colored space, we obtain a continuous function a on 
a compact domain X. 

The colorings have in general been suppressed 
from the notation of the paper. To obtain consis- 
tent results, a statement about x € [2, 2t ] should be 
interpreted as having an implicit condition saying 
that x has the coloring induced from the coloring 
on the component containing [2, 2to\. A later state- 
ment about y € [2to, V8) deals with y of a different 
color, and no relation between x and y of different 
colors is assumed at the endpoint 2tQ. 

In general, to the greatest extent possible, we ex- 
press each inequality as a strict inequality on a com- 
pact domain. There are, however, a few inequalities 
that are not strict, such as the bound of 1 pt on the 
score of a quasi-regular tetrahedron or the bound 
of on the score of a quad cluster. (These partic- 
ular sharp bounds appear in the proof of the local 
optimality of the face-centered cubic and hexagonal 
close packings.) 
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15 Appendix. 
Functions 



Constants and 



This appendix gives a summary of various constants 
and functions that appear in the proof. 

The density of a regular tetrahedron and Rogers's 
bound on the density of a packing is 

S tet = v / 8arctan(v / 2/5). 

The density of a regular octahedron is 

S oct = — « 0.72. 

The unit of measure in calculations of <r(D) is a 
point: 

pt = -tt/3 + V2S te t ~ 0.05537. 
A contravening decomposition star D satisfies 
a(D) >8pt 

or equivalcntly, 

t(D) < (4ttC - 8) pt, 

where 

C = l/(2arctan(V2/5)). 

Objects are truncated using the value to = 1.255. 
Centers of the packing at distance at most 2i from 
the origin form a set U(D). 

Internal angles of triangular standard regions are 
at least 0.8638, and those of other standard regions 
are at least 1.153. The internal angles of standard 
regions are generally denoted with variables a. The 
dihedral angles of simplices S are denoted dih(S'). 
If the edge of a simplex has index i, then we write 
dih^S*) for the dihedral angle along the edge with 
index i. 

The circumradius of a triangle is denoted 77. If 
there is a simplex with indexed edges, then r/ijk de- 
notes the circumradius of the face of the simplex 
whose edges have indices i, j, and k. 

In general, the edges of simplices arc indexed 
1 , . . . , 6 so that the first three edges meet at a ver- 
tex, and so that opposite edges have indices that 
are congruent mod 3. We write j/i, . . . , y& for the 
edge lengths and xi = yf as their squares. In gen- 
eral, the vertex where the first three edges meet is 
fixed as a distinguished vertex of the simplex. We 
write sol(S') for the solid angle of the simplex at 
that vertex. 



The function a on the topological space of de- 
composition stars is the function optimized in this 
paper. It is closely related to 

t(D) = sol(D)(pt-a(D). 

The function a is expressed as a sum of functions 
ctr, indexed by standard regions. These sums are 
in turn expressed as finer sums. 

The formulas for <jr and their upper bounds con- 
tain functions T, fi, v, quo, Adih, vor, vor . We do 
not repeat these definitions (they are found in [2]), 
because they appear in this paper in a purely formal 
way. (We make use of interval arithmetic inequali- 
ties that these functions satisfy, but we do not make 
use of anything deeper about them.) The function a 
is defined in [6, Section 3.11]. The function vo^ is a 
function on anchored simplices that are not upright 
quarters. It is vor or voro depending on whether 
the upright simplex has type C. Also, we set 



vot x = 



vor type A, C, C 
vor otherwise 



Corresponding to these functions related to a are 
functions related to r: 

t x = sol (pt — VOT x 

TR = SOI (pt - <T R 

f = sol (pt — u 

T = Sol (pt — G 

r = sol (pt — voro 
Tr = sol (pt — T 

T v = SOl (pt — V 

T> = SOl (pt - H 

Ty = sol (pt — vor 

t( ,t) — sol (pt — vor( , t). 

In the last case, the truncation is at t = \/2, 1.385, 
to, and so forth. In each case, sol is the function 
that measures the solid angle of the object chosen 
from the domain of the function. Each function on 
the left-hand side is referred to as a squander, and 
each function that are subtracted on the right-hand 
side is a score. When the function r* is related to 
the scoring function in this way, we say that it is 
adapted. 

The function voro is used in many contexts, and 
we sometimes take liberties with the notation for 
this function. For example, we write voro, r(U(D)) 
to indicate that the function depends only on the set 
U(D) and the combinatorial information encoded 
by the region (or subregion) R. 

The polynomial A of six variables has the prop- 
erty that A(xi, . . . , xq) < 0, then the simplex with 
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squared edge lengths Xi does not exist. In the case 
of equality A = 0, the simplex degenerates to a 
planar quadrilateral. 

16 Appendix. A Review of [6] 

16.1 Standard regions 

This paper relies heavily on results about the inter- 
nal structure of standard regions that is developed 
in [6]. We make a brief review of these structures. 
We restrict our attention to exceptional standard 
regions. A quarter is a simplex, whose vertices are 
vertices v € A in the packing, such that five edges 
lengths are in [2, 2t ] and the sixth edge length lies 
in the interval [2to, V&]- The sixth edge of the quar- 
ter is called its diagonal. 

We fix our attention on the decomposition star D 
centered at a given vertex in the packing, which we 
take to be located at the origin. The decomposition 
star determines the set U(D) of all vertices of A at 
distance at most 2i from the origin. 

From the decomposition star the set of quarters 
with a vertex at the origin can be also determined. 
There are two types. If the diagonal of the quarter 
has an endpoint at the origin, it is said to be up- 
right. Otherwise, it is said to be flat. A flat quarter 
has its vertices at the origin and three points of the 
set U(D). An upright quarter has one of its ver- 
tices v outside U(D). The edge (0, v) is an upright 
diagonal. 

The decomposition star determines a distin- 
guished set of quarters with vertex at the origin. 
These distinguished quarters are said to belong to 
the Q- system. If one quarter along a diagonal lies in 
the Q-system, then every quarter along that diago- 
nal lies in the Q-system. Quarters in the Q-svste 111 
do not overlap (in the sense that they have disjoint 
interiors) . The cones over the quarters at the origin 
in the Q-system do not overlap. 

Each quarter lies entirely in the cone over a single 
standard region. We fix our attention on a single 
standard region and the quarters that are associated 
with it in this way. We write Ur(D) for the subset 
of U (D) consisting of points whose projection to the 
unit sphere lands in the (closed) region R. 

There is a continuous function a defined on the 
topological space of all decomposition stars. It can 
be expressed as a sum of terms 

R 



indexed by the set of standard regions. We will de- 
scribe the structural properties of or for a given 
standard region R. There is a closely related func- 
tion 

t r = sol(R)Cpt- or. 

The theory develops in parallel for these two func- 
tions. 

16.2 Flat quarters 

Flat quarters appear in several ways. We color (in 
the sense of Appendix 14) each quarter with identi- 
fying information that marks how that flat quarter 
arises. We will define functions a and a a dj on the 
disjoint union of these different types of flat quar- 
ters. 

If a flat quarter is in the Q-system, then we set 

&(S) = o- adj (S) = fj,(S). 

Where fi(S) is the function defined in [2, For- 
mula 3.8]. 

Another type of flat quarter are those appearing 
in isolated pairs (see [2, Figure 1.10]). A special 
subset of such flat quarters are called simplices of 
type B. They have the following properties. 

1. The simplex is a flat quarter. 

2. The diagonal has length at least 2.77. 

3. The edges running between the endpoints of 
the diagonal and the origin have length at most 
2.23. 

For simplices of type B, set 

v(S) = Vadj{S) =vor(S, 1.385). 

For simplices that are part of an isolated pair, but 
that are not of type B, we set 

&(S) = o adj (S)=vor (S). 

There are flat quarters that enclose an upright 
diagonal along an upright diagonal in the Q-system. 
We call these masked flat quarters of the first type. 
In this case, there are two upright quarters Q\ and 
Q 2 over the flat quarter, as well as a third simplex 
S. We set 

fradj = v{Qi) + v{Q 2 ) + vor x (S), 

where vor x = vor if S has type C, and vor^ = vor 
otherwise. We note that 77456 (S) > \/2. Also, set 

&(S) = min(vor (S),0). 
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We have 

0~adj <&{S). 

[6, Lemma 3.11.1] 

Let (0, v, v\, v 2 ) be a flat quarter that is not any 
of those discussed above. Let v be the central ver- 
tex, and let (vi, v 2 ) be its diagonal. (Recall that the 
central vertex is defined as the vertex on a flat quar- 
ter that is not the origin and is not an endpoint of 
the diagonal.) These remaining flat quarters have 
the properties: 

1. \vi - v 2 \ > 2.6. 

2. \vi - v 2 \ > 2.7 or \v\ > 2.2. 

We call these quarters masked flat quarters of the 
second type. Set 

_ f 0.0063 M>2.2, 
6 [0.0114 M<2.2. 

Vadj(S) =vor (5)-e. 
a(S) =vor (S). 

For all types of flat quarters S, set 

t(S) = sol(S)(pt-a(S). 

We have 

O'adj < 

6~adj < Z(3, 1), 
Tadj >D(3,1), 

where Z{3, 1) = 0.00005 and £>(3, 1) = 0.06585. 
(See [6, 5.5].) 

16.3 Type A 

A simplex S of type A is a simplex with a vertex 
at the origin and at three other points v\ , v 2 , V3 in 
U R {D) such that 

1. The circumradius of the triangle with vertices 
vi,v 2 , V3 is less than 

2. Two of the edges (vi,vj) have edge length in 
the interval [2i ,2.77]. 

3. The third edge has length in the interval 
[2,2t ]. 



16.4 Loops 

Let (0,v) be an upright diagonal in the cone over 
a standard region R. We assume throughout this 
discussion that there is an upright quarter in the 
Q-system with diagonal (0, v). A vertex w (distinct 
from and v) of distance at most 2to from both 
and v is called an anchor (of v). By definition, we 
have w e U(D). 

Lemma 16.1. The projection of v to the unit 
sphere does not lie in the subregions associated with 
any flat quarter in the Q -system, any simplex of 
type A or any simplex of type B. The projection 
e of an edge (v,w) does not meet the interior of a 
subregion associated with any flat quarter in the Q- 
system, a simplex of type A, or a simplex of type 
B. 

The edges (v, w) from v to its anchors w emanate 
like spokes from v. We order the anchors of v ac- 
cording to this cyclic order into a cycle and write 
s v (w) for the successor of the anchor w in this or- 
dering. If the angle swept out along (0, v) in passing 
from one spoke to the next is greater than it, we say 
that the angle around (0, v) between w and s v (w) is 
concave. Otherwise, we say it the angle is convex. 

We say that a simplex (0, v, w, s v (w)) along (0, v) 
is an anchored simplex, if the angle is convex and if 
\w — s v (w)\ < 3.2. We say that (0, v, w, s v (w)) is a 
large gap, if the angle is convex and if \w — s v (w)\ > 
3.2. 

Two other anchored simplices S(yi, . . . , ye) that 
arise are those of type C or C. Write the anchored 
simplex S with edge yi representing the length of 
the upright diagonal. We say the anchored simplex 
has type C if j/4 € [2io, VS] and 

1. y 4 < 2.77, or 

2. Both face circumradii of S along the fourth 
edge are at most \/2. 

We say that the anchored simplex has type C if 

1. y 4 G [2t ,V8], 

2. y 2 ,y e e [2.45, 2t ], 

3. The fourth edge is a diagonal of a flat quarter 
that shares a face with S. 

If the anchored simplices around (0,v) are not 
contained in a half-space bounded by a plane 
through (0, v) and if there are no large gaps, then 
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the anchored simplices around v are said to form a 
loop around v. By construction, the anchored sim- 
plices around v do not overlap. Anchored simplices 
(even those around different vertices) in a contra- 
vening decomposition star do not overlap. An- 
chored simplices do not overlap flat quarters in the 
Q-system, nor do they overlap simplices of types A 
and B. 

A special simplex is a simplex (0, v\, V2, V3) with 
vertices in Ur(D) such that for some ordering of the 
indices, |«i-« 2 | < 2£ , \v2-v 3 \ < 2t , and U3I G 
[\/8, 3.2]. A special simplex (0,vi,V2,Vs) is said 
to lie along an anchored simplex (0,v,w,s v (w)) if 
the simplices can be indexed so that {vi,vs} = 
{w, s v (w)}, with \vi — Vs\ > \/8. Since anchored 
simplices do not overlap, each special simplex is re- 
lated to at most one upright diagonal (0, v) in this 
way. A special simplex attached to an upright di- 
agonal in this way does not overlap any of the pre- 
ceding objects (other special simplices so attached, 
other anchored simplices, simplices of types A and 
B, and flat quarters in the Q-system). 

16.5 Bounds on a R and r R . 

Let R be a standard region. More generally, we 
can partition R into subregions by considering three 
types of bounding edges: 

1 . a bounding edge of a standard region, 

2. the projection (to the unit sphere) of a face 
(along the diagonal) of an upright quarter in 
the Q-system, 

3. the projection of a face formed by and Vi,i>2 € 
U(D), where \vi — V2I < and additionally if 
it bounds a flat quarter and that flat quarter is 
not in the Q-system then the face also bounds 
a simplex of type A. 

Let R be one of these subregions. We have 

R 

If R is the projection of a simplex Q in the Q- 
system, then or(D) = o(Q), for some function a 
depending only on the simplex. Otherwise or(D) 
has the form 

o-r(D) = *r(V d ), 

where Vb is the V^-cell defined in [2, Section 2]. The 
V-cell can be truncated at y/2, which leads to the 



upper bound 

<tr{Vd) <vor fl (Vb,V^) <0. 
Each term indexed by quarters satisfies 
a(Q) < 0. 

In particular, if we have any subset Q of the set of 
quarters over a standard region R, then 

o-r(D) <]>>(Q). 
Q 

There is a second partition of subregions that 
is used more frequently. Let R be a subregion of 
the preceding paragraph. We break any anchored 
simplices that might occur from it (provided that 
the subregion does not come from a special simplex 
that lies along the anchored simplex). If this in- 
volves breaking the subregion into smaller pieces, 
we switch from the function aR to the upper bound 
vor (together with the volumes Sp(v) of [6, 2.11] 
. We also break off any flat quarters that are not 
masked (we switch to the upper bound a a dj on the 
flat quarter and to vor on the rest of the subre- 
gion). 

According to a construction that extends over 
several sections of [6] with a summary in [6, Sec- 
tion 3.10], it is often possible to eliminate the up- 
right quarters around an upright diagonal by re- 
placing the function or with an upper bound vor 
over those pieces. This process is called erasing the 
upright quarters. At other times, we get a weaker 
statement, by replacing the functions or for regions 
R at an upright diagonal (0, v) with 

a R( D ) ^ E vor , R (U(D)) + n(v, D). 

R R 

The constant 7r(w, D) is called a penalty. 

There are other upper bounds for or of inter- 
est. If we do not make approximations anywhere 
that produces penalties, we have a rather lengthy 
expression for or: 



o-r(D) = 





(6) 


+ E^vor(5) 


(7) 


+ E B vor(5) 


(8) 


+ Ei S olatcd, notB VOT 0(<3) 


(9) 


+ Eioo P :„Esv^ vor^S*) 


(10) 


■non-loop:v + Esh->u VOT a: 


(S))(ii) 


+ Etypcl ^adj 


(12) 


+ J2vov ,r>(Ur(D)). 


(13) 
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The sum indexed by Q runs over flat quarters in 
the Q-system. The sum indexed by A runs over sim- 
plices of type A over R. The sum indexed by B runs 
over simplices of type B over R. The next sum runs 
over isolated flat quarters that do not have type B. 
The next sum runs over the upright diagonals (0, v) 
of loops. (We restrict the sum to those loops that 
have four or five anchored simplices around it.) Its 
inner sum runs over the anchored simplices around 
(0, v), as well as the special simplices that lie along 
these anchored simplices. 

The next sum runs over upright diagonals (0, v) 
of upright quarters in the Q-system that are not 
part of a loop. The next sum runs over masked flat 
quarters of type 1 for which there enclosed upright 
diagonal (0, v) bounds exactly two quarters. The 
final sum is parameterized by subregions R 1 that 
have not otherwise been included, so that the sub- 
regions in the sum give a partition of R. 

Each summand in this formula for or defines a 
geometric object whose cone cuts out a sub-region 
R! on the unit sphere. We write cri,(D) for the 
summand. The summand is a single term except 
for those indexed by upright diagonals. The sum- 
mands indexed by upright diagonals include the full 
inner sum over anchored simplices and so forth. We 
attach constants (n, k) to each such subregion. 

To each of these subregions, we associate a natu- 
ral number n(R). For polygons, it is the number of 
edges. The edge is to be counted with multiplicity 
two if it does not bound the standard region. We 
connect an internal dot that is not connected to any 
other vertices to another vertex by an imaginary 
edge of length greater than 2to- Such an imaginary 
edge is counted with multiplicity two. The num- 
bers n(R) are indicated in Figure 18. This constant 
is defined for the subregions indexing the sums of 
Identity 13. 




Figure 18: Parameters (n, k) for non-polygonal re- 
gions 



A second parameter is k(R) the number of sides 
of length at least 2to, counted with multiplicity 2 if 
it is internal. (Count the imaginary edge if it exists, 
with multiplicity two.) 

Define constants 

U = 0.1317 t 5 = 0.27113 t 6 = 0.41056, 

t 7 = 0.54999 t s = 0.6045, 

s 5 = -0.05704 s 6 = -0.11408 s 7 = -0.17112, 
s 8 = -0.22816. 

Define constants 

D{n,k) = t n+k - 0.06585A;, 

for < k < n and tg > n + k > 4. It is readily 
verified that 

D(n 1 ,k 1 )+D(n 2 ,k 2 ) > D(m + n 2 -2, ki + k 2 -2). 
Define constants e — 0.00005 and 

Z(n, k) = s n+k - ke, 

for (n, k) ^ (3, 1) and Z(3, 1) = e. It is readily 
verified that 

Z(m,ki) + Z(n 2 , k 2 ) < Z(m +n 2 -2,k 1 + k 2 -2). 

Theorem 16.2. For subregions R' associated with 
summands as above, we have 

a+ <Z(n(R'),k(R')), 
sol{R')(pt - 0-+, > D(n(R'), k(R')). 

16.6 Bounds for Standard Regions 

Every standard region of a contravening decompo- 
sition star is a polygon with between 3 and 8 sides 
or one of the cases of Figure 18. In these figures, 
we draw a dot to indicate the projection of a vertex 
in U(D) that does not lie on the boundary. These 
enclosed dots do not lie in the cone over a simplex 
of type A, nor in the cone over a quarter. We draw 
an edge between vertices if the corresponding ver- 
tices in U(D) have distance at most 2to from each 
other. (This figure is a copy of Figure 8.) 

Corollary 16.3. Define constants as in Section 
16.5. For every standard region R such that n — 
n(R) > 5: 

O-R < S n , 

tr > t n . 

The corollary is an immediate consequence of the 
additivity properties of a^/, the sub and superad- 
ditivity of Z and D, and finally the geometric in- 
terpretation of sub and superadditivity as gluing 
subregions together to form a larger subregion. See 
[6, Section 4.5]. 
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Lemma 16.4. Let R be a subregion cut out by a 
loop group. Assume that the subregion is not that 
defined by a flat quarter. Let n be the number of 
anchors in the loop. Let k be the number of an- 
chored simplices, with an edge ( opposite the upright 
diagonal) of length at least 2i - Let s be the number 
of special simplices. Ln a contravening plane graph 
(n,k) is one of (4,0), (4,1), (4,2), (4,3), (5,0), 
(5,1). Furthermore, 

so\(R)Cpt - o-ft(D) > D(n + s,k - s) + Si oop (n, k), 

where Si oop is defined by the following constants. 
Similarly, 

&r(D) < Z(n + s,k - s) + zi 00 p(n, k). 



4. If the vertex has three anchors, the loop con- 
sists of three anchored simplices that cut out 
the same region on the unit sphere as the 
masked flat quarter. 

5. If the vertex has four anchors, the loop consists 
of three upright quarters and one other that is 
not an upright quarter. 

If the upright diagonal does not lie in the cone 
over the masked flat quarter, then the following 
facts are known. Call these masked flat quarters 
of the second type. 

1. The flat quarter has a diagonal of length at 
least 2.6. 



(n, k) 5i 00 p 

(4. 0) 0.0045 

(4. 1) 0.00272 

(4.2) 0.12034 

(4.3) 0.29426 

(5.0) 0.09537 

(5. 1) 0.24939 

(/I, Zloop 

(4,2) -0.08582 
(see [6, Section 5.11].) 

16.7 Remaining subregions 

A masked flat quarter is a flat quarter that is not in 
the Q-system and that overlaps some anchored sim- 
plex around an upright diagonal (0, v ). The vertices 
of a masked flat quarter consist of three consecutive 
anchors of v. 

There are two cases to consider. Either the up- 
right diagonal (0, v) lies in the cone over the masked 
flat quarter or it does not. If it lies in the cone over 
the masked flat diagonal, then the following facts 
are known. (Recall that we are assuming through- 
out this appendix that the standard region is not a 
quadrilateral.) As above, we call these masked flat 
quarters of the first type. 

1. The circumradius of the top face of the flat 
quarter is at least \[2. 

2. The anchored simplices form a loop. 

3. The vertex v has three or four anchors. 



2. If the flat quarter has a diagonal of length less 
than 2.7, then the edge opposite the diagonal 
has length at least 2.2. 

16.8 Penalties 

Set 

s ffr = 0.01561 2t < y < 2.57 

~\& = 0.00935 2.57 <y<V8. 

and 

n () = [O 2t < y < 2.57, 

\tv = 0.003521 2.57 < y < y/8. 

If R is one of the subregions of Identity 13, then 
we have a two functions voyo.r(U(D)) and a^(D). 
For example, R may be the projection of the set 
of anchored simplices around an upright diagonal 
(0, v), together with the set of special simplices that 
lie along the anchored simplices. In this situation, 
we write ir(v, D) for the difference: 

o-+(D)=voY ,R(U(D))+n(v,D). 

If (0, v) is associated with a loop with n anchored 
simplices, then we have 

D) < n7rr(H). 

If there is a loop around v with n anchored sim- 
plices and with k masked flat quarters of the second 
type, then 

tt(w, D) < (n - 2fc)7r r (M) + 2kw v (\v\). 
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We recall that for masked flat quarters Q of the 
second type, with projection to the subregion R, 
we have 

a{Q) = vor (Q). 

If there is a loop around v with a masked flat 
quarter of the first type, then the context is (n, k) — 
(3, 0) or (4, 1). In the case of context (4, 1), we have 

7r(«! D) < tt v (\v\). 

The constants ir(v, D) are called penalties. When 
we approximate <7r(D) by means of an upper bound 
on w(v, D), the upright diagonal is said to be erased. 

16.9 Anchored Simplices 

Often, ir(v, D) < 0. When this is the case, we gen- 
erally eliminate (that is, erase) the upright diagonal 
and all the anchored simplices around it. (This has 
been done implicitly already.) When we do this, we 
use the function vor on the region that becomes 
"exposed" in the process. It is possible to eliminate 
the upright diagonals in all cases, except 

1. loops, 

2. a configuration S3 , 

3. a configuration Sg", 

4. a configuration S4". 

16.9.1 The configuration S^ 

The configuration around an upright diagonal 
(0, v) consists of two upright quarters, one other 
anchored simplex, and one large gap. The configu- 
ration never masks a flat quarter. If the other 
anchored simplex is not an upright quarter, then 
the penalty satisfies 

tt(v,D) < 0.00222. 

If all three anchored simplices are upright quarters, 
then the penalty satisfies 

n(v,D) < 0.008. 



(0,v). If the configuration masks a flat quarter, 
then the penalty satisfies 

tt(v,D) < 7n-(|v|) + 27iy(M). 
If it does not mask a flat quarter, then 
n(v,D) < 37r r (M). 

There is at most one configuration in a contra- 
vening decomposition star. Each upright quarter Q 
is assigned a score a(Q). The three quarters satisfy 

£ (3) a(Q) < -0.4339, 
£ (3) r(Q) > 0.5606. 

16.9.3 The configuration Sj 

The configuration around an upright diagonal 
(0, v) consists of four upright quarters, and a large 
gap. The configuration masks at most k < 2 flat 
quarters. The penalty satisfies 

n(v, D) < 2kir v (\v\) + (4 - 2k)w T {\v\). 

There is at most one configuration in a contra- 
vening decomposition star. It does not appear in 
a contravening decomposition star with Sg". Each 
upright quarter Q is assigned a score cr(Q). The 
four quarters satisfy 

E ( 4)<K<2) <-o.25, 

£ (4) r(Q) >0.4. 

16.10 Compatibility Notes 

We conclude this section with a few comments 
about differences in notation between this article 
and the others. 

In [6], simplices of types A, B, C, and C are 
called Sa, Sb, Sc, and So, respectively. The 
functions <Jr(D) and tr(D) are denoted a(R) and 
t(R), respectively, in other articles. The trunca- 
tions vor 0j i{(-D) and t ^(D) are denoted vor (i?) 
and To(i?), respectively, in other articles. 



16.9.2 The configuration S 3 

The configuration S3 around an upright diagonal 
(0, v) consists of three upright quarters that all lie in 
a common half-space bounded by a plane through 
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17 Appendix. Calculations 



The proof of the Kepler Conjecture relies on a large 
number of inequalities that have been established 
by interval arithmetic with a computer. Generally, 
speaking there is no geometry of interest associated 
with these inequalities. The methods that are used 
in verifying these conjecture are described in [3] and 
[9]. (The paper describing the algorithms goes be- 
yond what is used in these verifications, because it 
also proposes some new unimplemented techniques 
that should simplify the verification process.) 

Interval calculations are arranged according to 
the section in which they appear. Each inequality 
is accompanied by one or more reference numbers. 
These identification numbers are needed to find fur- 
ther details about the calculation in [3]. Most of 
the verifications were completed by Samuel Fergu- 
son. 1 The verifications that he made are marked 
with a dagger (f). The computer-based proofs are 
carried out by interval arithmetic as in our earlier 
papers. A few inequalities are not standard in the 
sense that they have substantially higher dimension 
that the others. Linear programming methods are 
used to break these problems into six-dimensional 
pieces according to the suggestion of [5, Appendix 
2]- 

Edge lengths whose bounds are not specified are 
assumed to be between 2 and 2< - The first edge of 
an upright quarter is its diagonal. The fourth edge 
of a flat quarter is its diagonal. 

a is the function of [6, Section 3.11] and f = 
sol£pi — a. The function definitions are those ap- 
pearing in Appendix 15. The other functions are 
described in Appendix 15. 

We introduce some notation for the heights and 
edge lengths of a polygon. The heights will gener- 
ally be between 2 and 2£o, the edge lengths between 
consecutive corners will generally be 2, 2io, or 2\/2. 
We represent the edge lengths by a vector 

(ai, 61,02,62, ■ • ■ ,o„,6„), 

if the corners of an n-gon, ordered cyclically have 
heights Oj and if the edge length between corner i 
and i + 1 is 6j. Figure 19 illustrates the case n = 5. 
We say two vectors are equivalent if they are related 
by a different cyclic ordering on the corners of the 
polygon, that is, by the action of the dihedral group. 



1 I express my sincere thanks to Samuel Ferguson for his 
fundamental contributions to this project. 
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Figure 19: The labels on an n-gon. 
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17.1 Group 1 



t- 0.2529 dih> -0.3442, if y x e [2.3, 2t ], and dih > 1.51. (572068135) 

t - 0.2529dih > -0.1787, if y 1 € [2.3, 2i ], y 6 e [2-^/2,3.02], 1.26 < dih < 1.63. (723700608) 

f- 0.2529 dih 2 > -0.2137, if y 2 G [2.3, 2t ], y 4 G [2i ,2\/2], (560470084) 

t - 0.2529 dih > -0.1371, if y 1 € [2.3, 2i ], 2/5,2/6 G [2i ,3.02], 1.14 < dih < 1.51. (535502975) 



17.2 Group 2f 

A. dih < 1.63, if y 6 > 2f , 2/2,2/3 G [2,2.168]. (821707685) 

B. dih < 1.51, if 2/5 = 2t , 2/e > 2t , J/2,2/3 G [2,2.168]. (115383627) 

C. dih < 1.93, if 1/6 > 2t , 2/4 = 2\/2, 2/2,2/3 G [2,2.168]. (576221766) 

D. dih<1.77, if2/ 5 = 2i , 2/6 > 2^o, 2/4 = 2\/2, 2/2,2/3 G [2,2.168]. (122081309) 
t - 0.2529 dih > -0.2391, if y 6 > 2t , dih > 1.2, 2/2,2/3 G [2,2.168]. (644534985) 
t - 0.2529 dih > -0.1376, if y 5 = 2t , y e > 2t , dih > 1.2, and y 2 ,y 3 G [2,2.168]. (467530297) 
T O - 0.2529dih> -0.266, if y 6 > 2t , y 4 e [2i ,2\/2], dih > 1.2, y 2 ,y 3 G [2,2.168]. (603910880) 
t - 0.2529 dih > -0.12, if y 5 = 2t , 2/6 > 2i , y 4 G [2£ ,2\/2], dih > 1.2, 2/2,2/3 G [2,2.168]. (135427691) 
dih < 1.16, if 2/5 = 2to, ye > 2t 0) 2/4 = 2, 2/2, 2/3 G [2, 2.168]. (60314528) 
t - 0.2529 dih > -0.1453, if 2/2,2/3 e [2, 2.168], 2/5 G [2t , 3.488], y 6 = 2t . (312132053) 



17.3 Group 3 

dih 2 > 0.74, if 2/1 G [2t , 2.696], 2/2,2/3 G [2,2.168]. (751442360) 

t - 0.2529 dih > -0.2391, if A(j/f, 4, 4, 8, 2t§, y%) > 0, 2/2,2/3 G [2,2.168], 2/5 G [2t , 3.488]. (893059266) 

dih +0.5(2.402- 2/4) < tt/ 2 , if 2/5 > 2< o, 2/2,2/3 G [2,2.168]. (690646028) 



17.4 Group 4 

dih > 1.78, if 2/4 = 3.2, 2/1 G [2i ,2\/2], y 2 + 2/3 < 4.6. (161665083) 

17.5 Group 5 

The following inequalities hold for flat quarters. In these inequalities the fourth edge is the diagonal. 
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— clin 2 


i o Q£„. o 1 a i k„, i n 7noo„, n 1 7/., \ a ni oq 
+U.001/2 — U.loj/i — U.loj/3 + U. /Ozzt/5 — U.l ly 4 > — O.Olzo, 


ain 2 


— U.I0J/2 + U.bolyi + U.3I7/3 — 0. 007/5 + 0.41,37/4 + U.0zo7/g > z.oooo3, 


— ailii 


t f\ *71 A ~. A O0 1 „ . ft noi i AAA-, n Oil,. A A A 1 ~ . ^ AO/lOO 

+U.714j/i — O.zzIt/2 — 0.z21t/3 + 0.9zt/4 — O.ZZI7/5 — 0.z21t/6 > 0.3482, 


dihi 


A 01 £ i A OA*7A i A OA*7A 

— O.oloj/i + U.o97zj/2 + (J.o97zj/3— 




0.715t/ 4 + 0.3972t/5 + 0.3972t/ 6 > 2.37095, 


— sol 


— 0.187t/i — 0.187j/2 — 0.187t/ 3 + 0.1185t/ 4 + 0.479t/ 5 + 0.479t/ 6 > 0.437235f 


sol 


iA/IOO 1 A /I O O iA/IOO A O O /I AOO/1 O O /I /1 4. 

+(J.488j/i + 0.4881/2 + 0.4887/3 — 0.334t/5 — 0.334j/6 > 2.244y, 


— <7 


-0.145yi - O.O8I1/2 - O.O8I7/3 - 0.133t/ 5 - 0.133t/ 6 > -1.17401, 


—a 


-0.12yi - O.O8I7/2 - O.O8I7/3 ~ 0.113y 5 ~ 0.113y 6 + 0.029t/ 4 > -0.94903, 


a 


+0.153t/ 4 + 0.153i/5 + 0.1537/6 < 1.05382, 


a 


+0.419351 sol +0.19yi + 0.19t/ 2 + 0.19t/ 3 < 1-449, 


a 


+0.419351 sol < -0.01465 + 0.0436t/ 5 + 0.0436t/ 6 + 0.079431 dih, 


a 


< 0.0114, 


T 


> 1.019 pi, 



(867513567) 



17.6 Group 6 

In a quadrilateral cluster, with a given edge 7/4 as the diagonal, the other diagonal will be denoted y' 4 . 

The following relations for upright quarters (scored by v) hold. (We use the inequalities of [5, Section 4] 
for upright quarters in quad clusters, which are scored by a different function.) In these inequalities the 
upright diagonal is the first edge. We include in this group, the inequalities [6, A 2 ], [6, A 3 ] for upright 
quarters. 

2/i > 2i , 

2/1 < 2V2, 

dihi -0.636j/i + 0.462t/ 2 + 0.462t/ 3 - 0.82t/ 4 + 0.462t/ 5 + 0.462t/ 6 > 1.82419, 

- dihi +0.55i/i - 0.214i/2 - 0.214j/3 + 1.24t/ 4 - 0.214y 5 - 0.214t/ 6 > 0.75281, 
dih 2 +0.4yi - 0.15?/2 + 0.097/3 + 0.631t/ 4 - 0.57t/ 5 + 0.23t/ 6 > 2.5481, 

- dih 2 -0.454yi + 0.34t/ 2 + 0.154t/ 3 - 0.346t/ 4 + 0.805t/ 5 > -0.3429, 
dih 3 +0.47/1 - 0.157/3 + 0.09t/ 2 + 0.631t/ 4 - 0.57t/ 6 + 0.23t/ 5 > 2.5481, 

- dih 3 -0.4547/1 + 0.34t/ 3 + 0.154t/ 2 - 0.346t/ 4 + 0.805t/ 6 > -0.3429, 
sol +0.0657/2 + 0.065t/3 + O.O6I7/4 - 0.115t/ 5 - 0.115t/ 6 > 0.2618, 

- sol -0.2937/1 - 0.037/2 - 0.037/3 + 0.12t/ 4 + 0.325t/ 5 + 0.325t/ 6 > 0.2514, 
-v -0.05 387/2 - 0.05 38t/3 - 0.083t/ 4 - 0.0538t/ 5 - 0.0538t/ 6 > -0.5995, 
f < 0, (see [2, 3.13.3] and [2, 3.13.4]) 

t v -0.5945 pt>0. 

(498839271) 



V 


-4.10113 dihi 


< 


-4.3223, 


V 


-0.80449 dihi 


< 


-0.9871, 


V 


-0.70186 dihi 


< 


-0.8756, 


V 


-0.24573 dihi 


< 


-0.3404, 


V 


-0.00154 dihi 


< 


-0.0024, 


V 


+0.07611 dihi 


< 0.1196. 



([6, A 2 ]) 
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t v +4.16523 dihi > 4.42873, 
t v +0.78701 dihi > 1.01104, 
t„ +0.77627 dihi > 0.99937, 
t D +0.21916 dihi > 0.34877, 
t v +0.05107 dihi > 0.11434, 
t v -0.07106 dihi > -0.07749. 

([6, A 3 ]) 

17.7 Group 7 

The following additional inequalities are known to hold if the upright diagonal has height at most 2.696. 
i>t denotes the restriction of v to a simplex of compression type. 



yi 






< 2.696, 


(14) 


dihi 


-0.49J/1 


+ 0.44y 2 


+ 0.441/3 - 0.822/4 + 0.44y 5 + 0.Uy e > 2.0421, 


(15) 


-dihi 


+0.495yi - 


- 0.214y 2 


- 0.2142/3 + 1.052/4 - 0.2142/5 - 0.214y 6 > 0.2282, 


(16) 


dili2 


+0.38yi 


- 0.15y 2 


+ 0.092/3 + 0.542/4 - 0.572/5 + 0.24y 6 > 2.3398, 


(17) 


— dih 2 


-0.375j/i H 


- 0.33j/2 4 


- O.II2/3 - 0.362/4 + 0.722/5 + 0.0342/ 6 > -0.36135, 


(18) 


sol 


+0.42yi 4 


-0.165y 2 


■f 0.1652/3 - O.O62/4 - 0.1352/5 - 0.1352/ 6 > 1.479, 


(19) 


— sol 


-0.265yi - 


- 0.06y 2 - 


- O.O62/3 + 0.1242/4 + 0.2962/5 + 0.2962/ 6 > 0.0997, 


(20) 


— V 


+0.112?/!- 


0.1427/2 - 


- 0.1422/3 - O.I62/4 - 0.0742/5 - 0.074y e > -0.9029, 


(21) 


V 






+0.07611 dihi < 0.11, 


(22) 


v v 


-0.015j/i - 0.16(y 2 + y 3 + y 4 ) - 0.0738(y 5 + Ve) > -1.29285, 


(23) 


T v 






-0.07106 dihi > -0.06429, 


(24) 








> 0.0414. 


(25) 










(319046543) 



Remark 17.1. In connection with the Inequality (16), we occasionally use the stronger constant 0.2345 
instead of 0.2282. To justify this constant, we have checked using interval arithmetic that the bound 
0.2345 holds if 2/1 < 2.68 or 2/4 < 2.475. Further interval calculations show that the anchored simplices 
can be erased if they share an upright diagonal with such a quarter. 



17.8 Group 8f 

We give lower and upper bounds on dihedral angles. The domains that we list are not disjoint. In general 
we consider an edge as belonging to the most restrictive domain that the information of the following 
charts permit us to conclude that it lies in. 

The following chart summarizes the bounds. The dihedral angle is computed along the first edge. The 
chart is divided into three sections. In the first, there is no upright diagonal. In the second, the upright 
diagonal is the first edge. In the third, the upright diagonal is the third edge. The bounds in the second 
section have been established in [6, A 8 ]. 

In the first group 2/1, 2/2 , 2/3 € [2, 2to]. In the third row, the dihedral bound 1.624 holds for 2/6 on the larger 
interval [2i ,3.02]. In the seventh row, the dihedral bound 1.507 holds for 2/5,2/6 on the larger interval 



45 



[2i ,3.02]. 



2/5 

2, 2i ] 
2, 2* ] 
2, 2* ] 
2, 2to] 
2, 2t ] 
2, 2* ] 
2*o, 2 
2*0,2 
"2* ,2\/2 
2*o, 2 V2 



ye 

[2, 2* ] 
[2, 2* ] 
[2*o, 2 V2] 
[2*o, 2 V2] 
[2*o, 2 V2] 
[2*o, 2 V2] 
[2*0,2^2] 
[2*0,2^2] 
[2*o,2V2] 
[2*0,2^] 



2/4 

[2to, 2 v^] 

> 2V2 
[2, 2* ] 
[2* ,2\/2] 

> 2*o 

> 2\/2 
[2,2* ] 
[2* ,2\/2] 

> 2*o 

> 2^2 



dih mi „ 

1.153 

1.32 

0.633 

1.033 

1.033 

1.259 

0.817 

1.07 

1.07 

1.23 



dih max 

2.28 

2tt 

1.624 
1.929 
2tt 
2tt 

1.507 
1.761 

2tt 
2tt 



In the second group y\ e [2f ,2\/2], 2/2,2/3 G [2,2* ]. 



2/5 


2/6 2/4 


dih min 


dih max 


[2,2* ] 


[2, 2* ] [2, 2* ] 


0.956 


2.184 


[2, 2* ] 


[2,2* ] [2* ,2\/2] 


1.23 


7T 


[2, 2*o] 


[2, 2*o] > 2*o 


1.23 


7T 


[2, 2* ] 


[2, 2* ] > 2V2 


1.416 


7T 


2*o], 2/3 e [2*o, 2 s/2]. 






2/5 


2/6 2/4 


dih min 


dih max 


[2,2* ] 


[2, 2* ] [2, 2* ] 


0.633 


1.624 


[2,2* ] 


[2, 2* ] > 2*o 


1.033 


2tt 


[2,2* ] 


[2* ,2\/2] [2,2* ] 





1.381 


[2, 2*o] 


[2*0,2^] >2* 


0.777 


2tt 



(853728973) 



17.9 Group 9| 



(The verifications in this section that do not involve the score were verified by S. Ferguson.) 

Depending on the lengths of the edges 2/5, 2/6, 2/4, there are additional inequalities that hold. Generally 
these inequalities are for a fragment of a subregion. If V2, Vi, W3 are consecutive corners, we form a simplex 
S(yi,U2, 2/3,2/4,2/5,2/6) = (0, Vi,V2,Vs). There is no score intrinsically associated with this simplex, unless 
the subregion is triangular. But we can state various useful bounds on the dihedral angles along the first 
edge (0, vi). A triple (a, 6, c) preceding the inequality gives bounds on the edges (1/5, 2/6,2/4), respectively. 

In this group of inequalities, 2/1,2/2,2/3 <= [2,2*o]. Set vor^ — vor if the simplex has type A, and set 
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vot^ = vor , otherwise. 



(2/5,2/6,2/4) 

([2,2io],[2,2io],> 2V2) 

dih -0.372i/i + 0.465i/ 2 + 0.465i/ 3 + 0.465y 5 + 0.465i/ 6 > 4.885, 
([2,2i ],[2to,2x/2],[2i ,2^]) 

0.291yi - 0.393i/ 2 - 0.586i/ 3 + 0.79y 4 - 0.321y 5 - 0.397y 6 - dih < -2.47277, 
([2,2* ],[2i ,2x/2],> 2t ) 

0.291j/i - 0.393i/ 2 - 0.586i/ 3 - 0.321y 5 - 0.397i/ 6 - dih < -4.45567, 
([2,2* ],[2to,2x/2],> 2y/2) 

0.291j/i - 0.393y 2 - 0.586j/ 3 - 0.321y 5 - 0.397i/ 6 - dih < -4.71107, 

dih -0.214yi + 0.4y 2 + 0.58y 3 + 0.155y 5 + 0.395t/ 6 > 4.52345, 
([2t ,2V2l,[2t ,2V2],[2,2to]) 

r x >£>(3,2), see [6, Ai 6 ] 

vor* < Z(3,2), see [6, Ai 6 ] 

- sol -0.492yi - 0.492y 2 - 0.492i/ 3 + 0.43i/ 4 + 
0.038y 5 + 0.038y 6 < -2.71884, 

- vor^ -0.058yi - 0.105y 2 - 0.105y 3 - 0.115y 4 - 
0.062i/ 5 - 0.062y 6 > -1.02014, 

vor +0.419351 sol < 0.3085, 

0.115J/1 - 0.452y 2 - 0.452i/ 3 + 0.613y 4 - 0.15y 5 - 0.15y 6 - dih < -2.177, 
([2t ,2V2],[2t ,2V2], [2i ,2v^]) 

0.115t/i - 0.4521/2 - 0.452i/ 3 + 0.618i/ 4 - 0.15y 5 - 0.15y 6 - dih < -2.17382, 

vor < -0.121, see [6, Ai 8 ] 

r > D(3, 3) = 0.21301, see [6, A i7 ] 
([2t ,2v^|,[2to,2v^],>2to) 

0.115yi - 0.4521/2 - 0.452y 3 - 0.15y 5 - 0.15y 6 - dih < -3.725, 
([2t ,2V2],[2t ,2V2],> 2V2) 

0.115i/i - 0.4521/2 - 0.452y 3 - 0.15y 5 - 0.15i/ 6 - dih < -3.927. 



(529738375) 



17.10 Group 10 

In the next group 1/1 G [2t n ,2\/2], 1/2,2/3 G [2, 2i ] ■ Set vot^ = vor, if the simplex has type C or C. Set 
vor x = vor , otherwise. Other than this, retain the notation and conventions from Appendix 17.9. 

(2/5, 2/6, 2/4) 

([2,2i ],[2,2i ],[2to,2x/2]) 

vor x < 0, (cf. Appendix 17.23) 

0.471/1 - 0.522i/ 2 - 0.522i/ 3 + 0.812y 4 - 0.522i/ 5 - 0.522i/ 6 - dih < -2.82998, 
([2,2io],[2,2t<>],> 2t ) 

0.471/1 - 0.522i/2 - 0.522i/ 3 - 0.522i/5 - 0.522i/ 6 - dih < -4.8681, 

([2,2to],[2,2to],>2>/2) 

0.47i/i - 0.522y 2 - 0.522i/ 3 - 0.522i/ 5 - 0.522i/ 6 - dih < -5.1623. 



(456320257) 
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17.11 Group 11 

In the next group 2/1,2/2 G [2,2i ], J/3 G [2/j , 2-\/2]. Otherwise, retain the notation and conventions of 
Appendix 17.9. 

(2/5,2/6,2/4) 

([2,2t ],[2,2io],>2io) 

-0.42/3 + 0.152/1 - 0.092/2 - 0.631t/ 6 - 0.23y 5 - dih < -3.9788, 
([2,2t ],[2t ,2v^|,[2,2to]) 

0.2892/1 - O.U82/2 - 1.362/3 + O.6882/4 - O.U82/5 - l-36y 6 - dih < -6.3282, 
([2,2t ],[2*o,2\/2],> 2t ) 

0.289j/i - O.U82/2 - 0.723t/ 3 - O.U8y 5 - 0.723y 6 - dih < -4.85746. 



17.12 Group 12 



(664959245) 



v < -0.055 and r v > 0.092, 
if 2/1 G [2.696, 2^2], y 2 ,y 6 G [2.45, 2i ]. 



2. 



(7 < -0.039 and f > 0.094, 



if 2/2 e [2.45, 2t ], 2/4 G [2t ,2V2]. 



(704795925) 



(332919646) 



vor < -0.197 and t v > 0.239, 

if 2/1 G [2.696, 2\/2], 2/2 G [2.45, 2t ], j/ 6 G [2.45, 2t ], and 2/4 G [2i ,2\/2]. (These simplices have type 
S' c in the sense of [6, 2.10], and this justifies the use of the function vor.) (335795137) 



4. 



vor < -0.089 and r > 0.154, 
if 2/1 G [2.45, 2t ], 2/5,2/6 G [2t ,2s/2\. 



(605071818f) 



vor < -0.089 and r > 0.154, 
if 2/1 G [2.45, 2t ], 2/4,2/5 G [2t ,2s/2\. 



(642806938f) 



17.13 Group 13 

octavor < octavoro —0.017, if 2/1 G [2io, 2.696], and 77126 > V%- 

dih > 1.678, if 2/4 > 3.0, and 2/2 + 2/3 + 2/5 + 2/6 < 8.77, 



(104506452) 



(601083647) 



17.14 Group 14 

r < 0.3 1 38 - 0.1572/5, if 2/5 G [2, 2.138], 2/4 € [2t , 2.6], (543730647) 
r < -0.06, if 2/2 G [2.121,2.145], 2/4 G [2t Q ,2y/2], y 5 G [2.22,2.238]. (163030624) 

r < 10- 6 + 1.4 -O.I2/1 -0.15(2/2 + 2/3 + 2/5 + 2/6), if 2/4 G [2,2V2], 2/1,2/2,2/3 G [2,2.2], 2/5,2/6 G [2,2.35], 

(181462710) 
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17.15 Group 15 

vor < vor if y 4 G [2.7,2\/2], 2/1,2/2,2/3 € [2,2.14]. 
vor < -0.064, if 2/4 € [2*o, 2.72], 77455 > V2- 

vor < 1.0612 - 0.08(2/1 + 2/2 + 2/s) - 0.142(y B + y 6 ), if 2/4 G [2.7, 2\/2]. 
vor < -0.0713, if 2/4 G [2.59,2.64], 2/5 G [2.47, 2* ], 2/ 6 G [2.1,3.51]. 
vor < -0.06, if 2/1,2/2,2/3 G [2,2.13], 2/4 G [2.7,2.74], 77456 > V$. 
vor < -0.058, if y A G [2* , 2.747], 7/455 > \/2. 
vor < -0.0498, if 2/4 G [2t , 2.77], 77 456 > y/2. 



(463544803) 
(399326202) 
(569240360) 
(252231882) 
(472436131) 
(913534858) 
(850226792) 



17.16 Group 16 

We assume that 2/1,2/2,2/3 < 2.14 in the following inequalities. 

-r - 0.145yi - 0.08(2/2 + 2/3) - 0.133(2/5 + ye) > -1-146, if y 4 G [2* , 2\/2]. (594246986) 

-r-0.1452/i-0.081(2/2+2/3)-0.16(2/5+2/6) > -1.255, if 2/4 G [2* ,2\/2], and 2/5,2/6 G [2,2.3]. (381970727) 

-r - 0.032/1 - 0.03(2/2 + 2/3) - 0.094(2/5 + 2/ 6 ) > -0.5361. if y 4 G [2t , 2V2], j/ 5 + 2/6 > 4.3. (951798877) 

-r- 0.032/1 -0.03(2/2 + 2/3) -0.16(2/5 + 2/6) > -0.82-10- 6 . if 2/4 G [2t ,2\/2], 2/5 + 2/6 < 4.3. (923397705) 

r < -0.053 if 2/4 G [2t , 2V2], 2/5 > 2.35. (312481617) 

r < -0.041, if 2/4 G [2t ,2\/2], 2/5 G [2.25, 2rj ]. (292760488) 

r + 0.419351 sol < 0.079431 dih +0.0436 (y 5 + ye) - 0.0294, if 7/456 < V?, 2/4 G [2t 0) 2\/2]. (155008179) 

r < 1.1457 -0.1(2/1+2/2 +2/3) -0.172/5-0.112/6, if 2/1, 2/2, 2/3 G [2,2.13], 2/5 G [2,2.1], 2/6 G [2.27,2.43], 

2/4 G [2io,2.67]. (819450002) 

1.692/4 + 2/5 + 2/6 > 9.0659, if 2/4 G [2t , 2.7] and r/ 456 > V2. (495568072) 

1.692/4 + 2/5 + 2/6 > 9.044, if 2/4 G [2/ , 2.77] and 77455 > V2. (838887715) 

2/5 + 2/6 > 4.4, if t/456 > \/2, and 2/4 < 2.72. (794413343) 



17.17 Group 17 

In these inequalities 2/5,2/6 € [2t ,2-\/2] and 2/4 G [2,2rj ]. We distinguish three cases: 

1. The simplex has type A, so that 2/5,2/6 G [2/ ,2.77], and the scoring is vor, the analytic Voronoi 
function. 

2. 2/6 G [2.77, 2 V2]. 

3. The edges have lengths 2/5,2/6 G [2to,2.77], and 7/456 > V2- In the last two cases, the scoring is by 
the truncated Voronoi function voro. 

-vor-0.0582/i - O.O82/2 - O.O82/3 - O.I62/4 - 0.21(y 6 + ye) > -1.7531, if 2/1,2/2,2/3 < 2.14 and 2/5,2/6 G 
[2*0,2.77]. (378020227) 

- vor -0.0582/1 -O.I2/2 -O.I2/3 - 0.1652/4 - 0.1152/6 - 0.122/5 > -1.38875, if 7/456 > \/2, j/ 6 G [2.77, 2V2], 
and 2/1, 2/2, 2/3 < 2.14. (256893386) 

2/4 + 2/5 + 2/6 > 7.206, if 2/5, 2/6 G [2t , 2.77] and 77456 > V2. (749955642) 

-vor -0.0582/1 - 0.052/2 - 0.05y 3 - 0.16y 4 - 0.13y 6 - 0.13y 8 > -1-24547, if 2/1,2/2,2/3 < 2.14, 2/5,2/6 € 
[2*o, 2.77], and 77456 > V2. (653849975) 

vor < -0.077, if y A G [2* ,2\/2], 2/5 G [2.77,2^/2]. (480930831) 
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vor +0.419351 sol < 0.289, if y 4 ,y 5 e [2t ,2.77]. 

vor < 1.798 - 0.1 (yi + y 2 + ya) - 0A9y 4 - 0.17 fa + y 6 ), if y 6 , y 6 e [2.7, 2.77]. 



(271703736) 
(900212351) 



17.18 Group 18 

The next inequality is used as an estimate when there is an upright diagonal enclosed over a simplex of 
the third type in Appendix 17.17. 

vor < -0.078/2, y 1 e [2t , 2\/2], y 4 S [2t , 2.6961] for simpliccs of type C. (455329491) 

If there is no enclosed upright diagonal, then we can sometimes use the following instead. 

vor(S, \/2) < -0.078, if y 5 , y 6 € [2t , 2.6961], 77456 > a/2. (857241493) 

17.19 Group 19 

The interval calculations here show that the set of separated vertices (4) can be generalized to include 
opposite vertices of a quadrilateral unless the edge between those vertices forms a flat quarter. Consider 
a vertex of type (3,1,1) with a(3) = lApt. By the arguments in the text, we may assume that the 
dihedral angles of the exceptional regions at those vertices are at least 1.32 (see [6, 3.11.4]). Also, the 
three quasi- regular tetrahedra at the vertex squander at least 1.5 pt by a linear programming bound, 
if the angle of the quad cluster is at least 1.55. Thus, we assume that the dihedral angles at opposite 
vertices of the quad cluster are at most 1.55. A linear program also gives r + 0.316 dih > 0.3864 for a 
quasi-regular tetrahedron. 

If we give bounds of the form t x + 0.316 dih > b, for the part of the quad cluster around a vertex, where 
t x is the appropriate squander function, then we obtain 

T x > -0.316(2tt - 1.32) + b + 3(0.3864) 

for a lower bound on what is squandered. If the two opposite vertices give at least 2(1.4) pi+0.1317, then 
the inclusion of two opposite vertices in the separated set of vertices is justified. (Recall that t 4 = 0.1317.) 
The following inequalities give the desired result. 

t p + 0.316 dih > 0.5765, if dih < 1.55, y 4 € [2i ,2\/2]. (912536613) 
t + 0.316 dih > 0.5765 if dih < 1.55, y 4 > 2V2. (640248153) 
t v + 0.316 dih 2 > 0.2778, if yi e [2i ,2\/2], (594902677) 



17.20 Group 20f 

If the circumradius of a quasi-regular tetrahedron is > 1.41, then by [3, Section 9.17], r > 1.8 pt, and 
many of the inequalities hold (without further interval arithmetic calculations). 

In Sections 17.20 and 17.21, let Si, ... , S5 be 5 simplices arranged around a common edge (0,v), with 
\v\ e [2,2i ]. Let yi(Sj) be the edges, with yi(Sj) = \v\ for all j, y 3 (Sj) = y 2 (S j+1 ), and y 5 (Sj) = 
j/e(Sj+i). where the subscripts j are extended modulo 5. In Sections 17.20 and 17.21, ^dih(Sj) < 2ir. 
Set TTp = 2t; v + £ r if a = vor in the cases (y 4 > 2.6, y\ > 2.2) and (y 4 > 2.7). Set tt f — 0, otherwise. 

t(Si) +t(S 2 ) +t(S 4 ) > 1.4 pt, if y 4 (S 3 ),y 4 (S 5 ) > 2^2. (551665569) 

t(Si) +t(S 2 ) +t(S 3 ) > I Apt, iiy 4 {S 4 ),y 4 {S b ) > 2s/2. (824762926) 

r(S 1 ) + r(S 2 ) + (f(S 3 )-7r F ) + r(S 4 ) > 1.4p*+£>(3, 1), if y 4 (S 3 ) € [2t , 2s/2], y 4 (S 5 ) > 2t , dih(S 5 ) > 1.32, 

(675785884) 

r(S 1 ) + T(S 2 )+T(S 3 ) + (f(S 4 )-7r F ) > lApt+D(3, 1), if y 4 (S 4 ) G [2t , 2^], y 4 (S 5 ) > 2t , dih(S 5 ) > 1.32. 

(193592217) 
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17.21 Group 21f 

As in Section 17.20, the quasi-regular tetrahedra are generally compression scored. Define np as in 
Section 17.20. The constraint X)(5) dih(Sj) = 27r is assumed. 

t(Si) + t(S 2 ) + t(S 3 ) + t(S 4 ) > 1.5 pt, if i/ 4 (S B ) > 2V2. (325738864) 
t(Si)+t(S 2 )+t(S 3 )+t(S 4 ) + (t(S 5 )-^f) > 1.5pt+D(3,l), if t/ 4 (S 5 ) G [2*o,2\/2]. (314974315) 

17.22 Group 22 

If there are four quasi-regular tetrahedra {Si, . . . , S4} at the central vertex v of the flat quarter Q, and 
if there are only five standard regions at v, then 

&(Q)+J2°(Si) < 0-114. 

(4) 

(867359387) 



17.23 Group 23 



Let vor x = vor for simplices of types C and C. Let vor x = vor, otherwise. 

'0, yi G [2i ,2\/2] 

vor,; < I -0.05, yi € [2i , 2.696], 

-0.119, yi € [2i , 2.696], m2 6 > A 

if 2/4 G [2*o, 2\/2]. The bound of is established in Appendix 17.10. (Even when simplex has type C, the 
bound —0.05 is based on the upper bound vor .) (365179082) 

17.24 Group 24 

' I -0.043, yi e [2to,2.696], v ; 

for quad regions i? constructed from an anchored simplex S and adjacent special simplex S'. Assume 
that 2/4 (S) = 2/4 (S') € [2\/2, 3.2], and that the other edges have lengths in [2,2*o]- The bound is found 
in [2, Lemma 3.13]. The bound —0.043 is obtained from deformations, reducing the inequality to the 
following interval calculations. 

vor < -0.043/2, if y 6 = 2t , y 1 e [2t , 2.696]. (368244553+) 

voro(S)+vor (S(2, 2 / 2 ,2/3,2y4,2,2)) < -0.043, if Vl e [2t , 2.696], 2/4 G [2^2,3.2], y' 4 > 2i . (820900672+) 

vor (S) + vor (S(2t , 2/2,2/3,^4, 2, 2)) < -0.043, if 2/1 G [2t 0) 2.696], 2/4 G [2^2,3.2], y' 4 > 2i . 

(961078136f) 

Under certain conditions, Inequality 26 can be improved. 

vor (S') < -0.033, if dih(S') < 1.8, 2/1 G [2,2.12], 2/4 G [2^/2,3.2], 
vor (S) < -0.058, if dih(S) < 2.5, 2/1 G [2* , 2.696], 2/4 G [2^,3.2], 
voro(S) < -0.073, if 77126 > y/2, yi e [2t , 2.696], 2/4 G [2\/2,3.2].+ 

(The last of these was verified by S. Ferguson.) (424186517) 



51 



These combine to give 

,™ / si/ \ [—0.091, or 

J y ' [-0.106, 

for the combination of special simplex and anchored simplex under the stated conditions. 
17.25 Group 25 (pentagons) 

There are a few inequalities that arise for pentagonal regions. 

Proposition 17.2. If the pentagonal region has no flat quarters and no upright quarters, the subregion 
F is a pentagon. It satisfies 

vor < -0.128, 
r > 0.36925. 

Proof. The proof is by deformations and interval calculations. If a deformation produces a new flat 
quarter, then the result follows from [6, Ai 3 ] and Inequality 18.2. So we may assume that all diagonals 
remain at least 2\/2. If all diagonals remain at least 3.2, the result follows from the tec-bound on the 
pentagon [6, Section 5.5]. Thus, we assume that some diagonal is at most 3.2. We deform the cluster 
into the form 

(ai, 2, 012,2,03,2,04,2,05, 2), \vi\ = a,e {2,2t }- 

Assume that \v\ — vz\ < 3.2. If max(oi,a 3 ) = 2i , the result follows from [6, A13] and Section 18.2, 
Equations 31 and 32. 

Assume oi = a 3 = 2. There is a diagonal of the quadrilateral of length at most 3.23 because 

A(3.23 2 ,4,4,3.23 2 ,4,3.2 2 ) < 0. 

The result now follows from the following interval arithmetic calculations. 
(These inequalities are closely related to [6, A 2 i].) 

voro(S(2,2,2,y4,2,2))+voro(5(2,2,2,y 4 ,2,2)) + voro(5(2,2,2,y4,y 4 ,2)) < -0.128, if y 4 e [2^2,3.2], 
y' A e [2\/2,3.23]. (587781327+) 

r (5(2, 2,2,y 4 ,2,2)) + t (S(2, 2, 2, j/ 4 , 2, 2)) + t (S(2, 2, 2, 2/4, 2/ 4 , 2)) > 0.36925, if y 4 e [2^2,3.2], y' A e 
[2V2,3.23]. (807067544f) 

t (2,2,2/ 3 ,2/4,2,2) < ro(2to,2,y 3 ,y4,2,2), if y 4 E [2^2,3.06]. (986970370+) 

vor (2, 2, 2/3, 2/4, 2, 2) > vor (2t , 2, y 3 , y 4 , 2, 2), if y A e [2 A 3.06]. (677910379+) 

vor < -0.128, if yi = y 2 = 2/4 = 2, y 5 ,y 6 e [3.06,3.23]. (276168273) 

r > 0.36925, if 2/1=2/2= 2/3 = 2/4 = 2,y 5 ,y 6 E [3.06, 3.23]. (411203982) 

t > 0.31, if 2/1 = 2/2 = 2/4 = 2, 2/3 = 2i , 2/5,2/6 € [3.06,3.23]. (860823724) 

vor < -0.137 - (2/5 - 2\/2)0.14, if 2/1=2/4 = 2, 2/5 € [2\/2, 3.23], 2/6 € [3.06, 3.23]. (353116955) 

t > 0.31 + (2/5 - 2^2)0.14, if 2/1=2/4 = 2, 2/5 € [2V2, 3.23], 2/6 € [3.105, 3.23]. (943315982) 

t > 0.31 + (2/5 - 2^2)0.14 + (2/6 - 3.105)0.19, if 2/1 = 2/4 = 2, 2/5 € [2^2,3.23], y 6 e [3.06,3.105]. 

(941799628) 

vor < 0.009 + (2/5 - 2V2)0.14, if 2/1=2/4 = 2, 2/5 € [2V2, 3.23], 2/2 = 2/6 = 2. (674284283) 
t > 0.05925 + (2/5 - 2\/2)0.14, if 2/1=2/4 = 2, 2/5 € [2\/2, 3.23], 2/2=2/6 = 2. (775220784) 
t > 0.05925, if 2/3 = 2i , 2/5 = 2/6 = 2, 2/4 e [2\/2, 3.23]. (286076305) 
t > -(2/4 - 3.105)0.19, if 2/1 = 2to, 2/2 = 2/3 = 2/5=2/6 = 2, 2/4 G [3.06,3.105]. (589319960) 

□ 
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17.26 Group 26f 

Let Q be a quadrilateral region with parameters 

(ai, 2t , a 2 , 2, a 3 , 2, a 4 , 2i ), a, G {2, 2i }- 
Assume that |v 2 - v 4 \ G [2^2,3.2], |t>i - u 3 | G [3.2,3.46]. Note that 

A(4,4,8,2tg,2tg,3.46 2 ) < 0. 

vor (Q) < -0.168. (302085207) 
t (Q) > 0.352. (411491283) 

17.27 Group 27 

Consider a pentagonal region. If the pentagonal region has one flat quarter and no upright quarters, 
there is a quadrilateral region F. It satisfies 

vor < -0.075, 
r > 0.176. 

Break the cluster into two simplices 5* = S(yi, . . . , ye), S' = S(y[, j/ 2 , J/3, J/4, J/5, Vq), by drawing a diagonal 
of length 2/4. Assume that the edge y' 5 G [2to, 2y/2]. Let y' 4 be the length of the diagonal that crosses 3/4. 

vor < 2.1327 - 0.1j/i - 0.15j/ 2 - 0.08j/ 3 - 0.15j/ 5 
-0.15y 6 - 0.1j/i - 0.17y 5 - 0.16y 6 , 

if dih(S) < 1.9, dih(S') < 2.0, y x G [2,2.2], j/ 4 > 2V2, 
vor < 2.02644 - 0.1j/i - 0.14(y 2 + j/ 3 ) - 0.15(j/ 5 + y 6 ) - O.lj/i - 0.12(j/ 5 + 

if j/i G [2,2.08], j/ 4 <3. 
vor +0.419351 sol < 0.4542 + 0.0238(j/ 5 + j/ 6 + 2/ 6 ), 

if J/4,J/ 4 >2V2. 

The inequalities above are verified in smaller pieces: 

vor < 1.01-0.1j/i-0.05j/2-0.05y 3 -0.15j/5-0.15j/ 6 , if dih < 1.9, y 4 > 2\/2, and yi < 2.2. (131574415) 

vor < 1.1227 - 0.1j/i - 0.1y 2 - 0.03j/ 3 - 0.1 7y 5 - 0.16j/ 6 , if dih < 2, j/ 4 > 2\/2, J/5 G [2i ,2\/2], and 
2/2 + J/3 < 4.67. (929773933) 

vor < 1.0159 -0.1j/i-0.08(j/ 2 + j/ 3 ) + 0.04j/ 4 -0.15(j/ 5 + j/ 6 ), if J/4 G [2V2,3], j/i G [2,2.08]. (223261160) 

vor < 1.01054-0.1j/i-0.06(j/ 2 +j/ 3 )-0.04j/4-0.12(j/ 5 +j/ 6 ), if j/ 4 G [2>/2,3], y 5 G [2t ,2\/2]. (135018647) 

Let Q be the quadrilateral subcluster expressed as a union of two simplices S and S', as above 

voro(Q) + 0.419351 sol(Q) < 0.4542 + 0.0238(j/ 5 + y 6 + y' 6 ), 

if y' 5 G [2t , 2\/2] and both diagonals have length > 2\/2. (559676877) 

17.28 Group 28f 

If v is a vertex of an exceptional cluster and there are exactly 4 quasi-regular tetrahedra along (0,v), 
then there are 5 vertices V\, . . . ,v 5 adjacent to v. We have 

J2(\ v ~ v i\ + N) > 20 - 42 - 

(5) 

(615073260) 
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If v\ and V5 are the two vertices on an exceptional cluster, and if the additional hypothesis \v\ — v$\ > 2\/2 
holds, then 

J2(\v-Vi\ + \ Vi \) > 20.76. 

(5) 

(844430737) 



17.29 Group 29 

vor < -0.136 and t > 0.224, 

for a combination of anchored simplex S and special simplex S', with yi(S) € [2.696, 2\/2], 2/2(5), 2/6(5) e 
[2.45, 2io], 2/4(5*) G [2\/2, 3.2], and with cross-diagonal at least 2t - This inequality can be verified by 
proving the following inequalities in lower dimension. In the first four 2/1 G [2.696, 2\/2], 2/2, 2/6 £ [2.45, 2io], 
2/4 G [2\/2, 3.2], and 2/4 > 2i (the cross-diagonal). 



vor (S(j/i, 


... ,2/ 6 ))+vor (5(2, y 2 ,y 3 ,2/4,2, 2)) < -0.136 


(967376139+) 


vor (S , (2/i, 


... ,2/ 6 ))+vor (5(2t , 2^2,2/3,2/4, 2, 2)) < -0.136 


(666869244+) 


7-0(5(2/1, • • 


. , y 6 )) + r (5(2, 2/2, 2/3, 2/4, 2, 2)) > 0.224 


(268066802+,) 


7-0(5(2/1, • • 


• ,2/e)) +T (5(2to,2/2,2/3,2/4,2,2)) > 0.224 


(508108214f) 


vor (5(?/i, 


... ,2/ 6 )) < -0.125, if 2/1 G [2.696,2^2], 2/2,2/6 € [2.45, 2t ], 2/5 = 2t . 


(322505397f) 


vor (5(2/i, 


... ,2/ 6 )) < 0.011, if 2/1 G [2.696, 2\/2], 2/5 = 2t . 


(736616321f) 


to(5(2/i, • • 


. ,2/ 6 )) >0.17, if 2/1 G [2.696, 2V2], 2,2,2/6 € [2.45, 2t ], 2/5 = 2t . 


(689417023+) 


7-o(5(2/i, ■ • 


. ,2/ 6 )) > 0.054, if 2/1 G [2.696, 2V2], 2/5 = 2t . 


(748466752f) 


17.30 


Group 30 






vor < -0.24 and r > 0.346, 





for an anchored simplex S and simplex S' with edge parameters (3, 2) in a hexagonal cluster, with 
2/2(5) = 2/2(5'), 2/3(5) = 2/3(5'), 2/4(5) = 2/4(5'), 2/i(5) e [2.696, 2V2], y 4 (S) G [2^2,3.2], y 2 (S),y 6 (S) G 
[2.45, 2t ], and 

max(2/ 5 (5'),2/ 6 (5')) G [2i ,2\/2], min(2/ 5 (5'), 2/ 6 (5')) G [2,2t ]. 

This breaks into separate interval calculations for 5 and 5'. 
This inequality results from the following four inequalities: 

vor (5) < -0.126 and t (5) > 0.16 (369386367+) 

vor (S") < —0.114 and r (5') > 0.186 (There arc two cases for each, depending on which of 2/5,2/6 is 
longer.) (724943459+) 

17.31 Group 31 

vor < -0.149 and t > 0.281, 
for a quadrilateral subcluster with both diagonals > 2\/2, and parameters 

(01,61,02,62,03,63,04,64), 
with o 4 G [2.45, 2t ], and 64 G [2t , 2\/2]. (836331201f) 
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17.32 Group 32 



vor < -0.254 and t > 0.42625, 

for a combination of anchored simplex S and quadrilateral cluster Q. It is assumed that yi(S) G 
[2.696, 2\/2], y2(S),ye(S) e [2.45, 2t ]- The adjacent quadrilateral subcluster is assumed to have both 
diagonals > 2y/2, and parameters 

(01,61,02,62,03,63,04,64), 

with 64 G [2\/2, 3.2]. The verification of this inequality reduces to separate inequalities for the anchored 
simplex and quadrilateral subcluster. For the anchored simplex we use the bounds voro(S") < —0.126, 
To(S') > 0.16 that have already been established above. We then show that the quad cluster satisfies 

vor < -0.128 and t > 0.26625. (327474205+) 

For this, use deformations to reduce either to the case where the diagonal is 2\/2, or to the case where 
61 = 6 2 = 63 = 2, o 2 , a 3 e {2, 2t }. When the diagonal is 2\/2, the flat quarter can be scored by [6, A i3 ]: 
(voro < 0.009, t > 0.05925). (There are two cases depending on which direction the diagonal of length 
runs.) 
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18 Appendix Hexagonal In- 
equalities 

There are a number of inequalities that have been 
particularly designed for standard regions that are 
hexagons. This appendix describes those inequali- 
ties. They are generally inequalities involving more 
than six variables, and because of current techno- 
logical limitations on interval arithmetic, we were 
not able to prove these inequalities directly with 
interval arithmetic. 

Instead we give various lemmas that deduce the in- 
equalities from inequalities in a smaller number of 
variables (small enough to prove by interval arith- 
metic.) 

18.1 Statement of results 

There are a number of inequalities that hold in spe- 
cial situations when there is a hexagonal region. Af- 
ter stating all of them, we will turn to the proofs. 

1. If there are no flat quarters and no upright 
quarters (so that there is a single subregion F), 
then 

vor < -0.212 (27) 
r > 0.54525. (28) 



2. If there is one flat quarter and no upright quar- 
ters, there is a pentagonal subregion F. It sat- 
isfies 

voi'o < -0.221 
r > 0.486. 

3. If there are two flat quarters and no upright 
quarters, there is a quadrilateral subregion F. 
It satisfies 

vor < -0.168, 
r > 0.352. 

These are twice the constants appearing in 11; 

4. If there is an edge of length between 2to and 
2^/2 running between two opposite corners of 
the hexagonal cluster, and if there are no flat or 
upright quarters on one side, leaving a quadri- 
lateral region F, then F satisfies 

vor < -0.075, 
r > 0.176. 



5. If the hexagonal cluster has an upright diago- 
nal with context (4, 2), and if there are no flat 
quarters (Figure 20, then the hexagonal cluster 
R satisfies 

a R < -0.297, 
t r > 0.504. 




Figure 20: A hexagonal cluster with context (4,2). 



6. If the hexagonal cluster has an upright diago- 
nal with context (4, 2), and if there is one un- 
masked flat quarter (Figure 21, let {F} be the 
set of four subregions around the upright diag- 
onal. (That is, take all subregions except for 
the flat quarter.) In the following inequality 
and Inequality 7, let er^ be defined as cjr on 
quarters, and vor x on other anchored simplices. 
is the adapted squander function. 

E(4)4 < -0.253, 
£ (4) r+ > 0.4686. 




Figure 21: A hexagonal cluster with context (4,2). 



7. If the hexagonal cluster has an upright diag- 
onal with context (4,2), and if there are two 
unmasked flat quarters (Figure 22, let {F} be 
the set of four subregions around the upright 
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diagonal. (That is, take all subregions except 
for the flat quarters.) 

E ( 4)4 < -o.2, 




Figure 22: A hexagonal cluster with context (4,2). 

8. If the hexagonal cluster has an upright diag- 
onal in context (4, 1), and if there are no flat 
quarters, let {F} be the set of four subregions 
around the upright diagonal. Assume that the 
edge opposite the upright diagonal on the an- 
chored simplex has length at least 2y/2. (See 
Figure 23.) 

vor Q , R (D) + £ (3) <j{Q) < -0.2187 
tom(D) + E( 3 )t(Q) > 0.518. 




Figure 23: A hexagonal cluster with context (4, 1). 

9. In this same context, let F be the pentagonal 
subregion along the upright diagonal. It satis- 
fies 

vor < -0.137, (29) 
t > 0.31. (30) 

10. If the hexagonal cluster has an upright diagonal 
in context (4, 1), and if there is one unmasked 



fiat quarter, let {F} be the set of four sub- 
regions around the upright diagonal. Assume 
that the edge opposite the upright diagonal on 
the anchored simplex has length at least 2^/2. 
(There are five subregions, shown in Figure 24.) 

vor , R {D) + £ (3) a(Q) < -0.1657, 
^(0)tE (3) t(Q) > 0.384. 




Figure 24: A hexagonal cluster with context (4, 1). 

11. In this same context, let F be the quadrilateral 
subregion in Figure 24. It satisfies 

vor < -0.084, 
r > 0.176. 

18.2 Proof of inequalities 

Proposition 18.1. Inequalities 1 -11 are valid. 

We prove the inequalities in reverse order 11- 1. 
The bounds vor < 0.009 and r > 0.05925 from 
[6, A13] for what a flat quarters with diagonal v° 
will be used repeatedly. Some of the proofs will 
make use of tec-bounds, which are described in [6, 
Section 5.2]. 

Proof. {Inequality 10 and Inequality 11.) Break the 
quadrilateral cluster into two simplices S and S' 
along the long edge of the anchored simplex S. The 
anchored simplex S satisfies t(S) > 0, a(S) < 0. 
The other simplex satisfies tq(S') > 0.176 and 
voro(5") < —0.084 by an interval calculation 18.16. 
This gives Inequality 11. For Inequality 10, we 
combine these bounds with the linear programming 
bound on the four anchored simplices around the 
upright diagonal. From the inequalities [6, A2] - 
[6, A7], [6, A22], [6, A24], we find that they score 
< -0.0817 and squander > 0.208. Adding these to 
the bounds from Inequality 11, we obtain Inequality 
10. □ 
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Proof. (Inequality (8) and (9).) The pentagon is a 
union of an anchored simplex and a quadrilateral 
region. LP-bounds similar to those in the previous 
paragraph and based on the inequalities of [6] show 
that the loop scores at most —0.0817 and squanders 
at least 0.208. If we show that the quadrilateral 
satisfies 

vor < -0.137, (31) 
T > 0.31, (32) 

then Inequalities (8) and (9) follow. If by defor- 
mations a diagonal of the quadrilateral drops to 
2\/2, then the result follows from Inequalities 18.15 
and [6, Ai 3 ]. By this we may now assume that the 
quadrilateral has the form 

(01,2,02,2,03,2,04,64), a 2 ,a 3 e {2,2* }. 

If the diagonals drop under 3.2 and max(a 2 ,o 3 ) = 
2io, again the result follows from Inequalities 18.15 
and [6, A13]. If the diagonals drop under 3.2 and 
02 = 03 = 2, then the result follows from [6, A19]. 
So finally we attain by deformations 64 = 2\/2 with 
both diagonals greater than 3.2. But this does not 
exist, because 

A(4,4,4,3.2 2 ,4,8,3.2 2 ) < 0. 

□ 

Proof. (Inequality 5, Inequality 6, and Inequality 
7.) Inequalities 7 are derived in [6, Section 5.11]. 
Inequalities 5, 6 are LP-bounds based on [6, A 2 ], 
[6, A 7 ], [6, A 22 ]. 

Proof of 4- Deform as in [6]. If at any point a 
diagonal of the quadrilateral drops to 2\/2, then 
the result follows from [6, A13] and Inequality 11: 

vor < 0.009 - 0.084 = -0.075, 
T > + 0.176 = 0.176, 

Continue deformations until the quadrilateral has 
the form 

(01, 2, o 2 , 2, 03, 2, o 4 , 64), 02, o 3 e {2, 2* }. 

There is necessarily a diagonal of length < 3.2, be- 
cause 

A(4,4,3.2 2 ,8,4,3.2 2 ) < 0. 

Suppose the diagonal between vertices v 2 and w 4 
has length at most 3.2. If a 2 = 2* or a 3 = 2f , 
the result follows from [6, A13] and Inequality 11. 
Take a 2 = a 3 = 2. Inequality 4 now follows from 
[6, A19]. □ 



Proof. (Inequality 3). We prove that the quadrilat- 
eral satisfies 

vor < -0.168 
t > 0.352. 

There are two types of quadrilaterals. In (a), there 
are two flat quarters whose central vertices are op- 
posite corners of the hexagon. In (b), the flat quar- 
ters share a vertex. We consider case (a) first. 

Case (a). We deform the quadrilateral as in [6]. If 
at any point there is a diagonal of length at most 
3.2, the result follows from Inequality 10 and In- 
equality 11. Otherwise, the deformations give us a 
quadrilateral 

(oi, 2, o 2 , 2i , 03, 2, 04, 2), Oj G {2, 2t Q }. 

The tec approximation now gives the result (see [6, 
5.3)]. 

Case (b). Label the vertices of the quadrilateral 
V\,... ,1)4, where (^1,^2) an d (vi,tti) are the di- 
agonals of the flat quarter. Again, we deform the 
quadrilateral. If at any point of the deformation, 
we find that \v\ — t> 3 | < 3.2, the result follows 
from Inequalities 10, 11. If during the deformation 
\v2 — ^4 1 < 2\/2, the result follows from [6, A13] 
and the interval calculations 18.14. If the diagonal 
(^2,^4) has length at least 3.2 throughout the de- 
formation, we eventually obtain a quadrilateral of 
the form 

(oi, 2t , o 2 , 2, o 3 , 2, o 4 , 2*o ), a, e {2, 2* }. 
But this does not exist: 

A(4,4,3.2 2 ,(2* ) 2 ,(2*o) 2 ,3.2 2 ) < 0. 

We may assume that |i>2 — W4I G [2\/2, 3.2]. The 
result now follows from interval calculations 18.14. 

□ 

Proof. (Inequality 2). This case requires more ef- 
fort. We show that 

vor < -0.221 
t > 0.486 

Label the corners (vi,...,vs) cyclically with 
(vi, V5) the diagonal of the flat quarter in the hexag- 
onal cluster. We use the deformation theory of [6]. 
The proof appears in steps (1), . . . , (6). 

(1) If during the deformations, \vi — v±\ < 3.2 or 
\ v 2 — W5I < 3.2, the result follows from Inequalities 
18.2 and 11. We may assume this does not occur. 
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(2) If an edge {vi,v 3 ), {v 2 ,v 4 ), or (v 3 ,v 5 ) drops to 
2^/2, continue with deformations that do not fur- 
ther decrease this diagonal. If \v\— v 3 \ = \v 3 — v 5 \ = 
2\[2, then the result follows from [6, A13] and in- 
terval calculations 18.7. 

If we have \vi — v 3 | = 2\/2, deform the figure to the 
form 

(ai, 2, a 2 , 2, a 3 , 2, a 4 , 2, a 5 , 2t ), a 2 , a 4 , a 5 e {2, 2i }- 

Once it is in this form, break the flat quarter 
(0,V\,V2, V3) from the cluster and deform w 3 until 
a 3 e {2,2t }- The result follows from an interval 
calculation 18.8. 

We handle a boundary case of the preceding cal- 
culation separately. After breaking the flat quarter 
off, we have the cluster 

(01, 2V2, a 3 , 2, a 4 , 2, a 5 , 2t ), a 3 , a 4 , a 5 E {2, 2t }. 

If |«i — v 4 \ = 3.2, we break the quadrilateral cluster 
into two pieces along this diagonal and use inter- 
val calculations 18.9 to conclude the result. This 
completes the analysis of the case \vi — v 3 \ = 2\p2. 

(3) If \vi — v 4 \ < 3.2, then deform until the cluster 
has the form 

(ai, 2, 012,2,03,2,04,2,05,2*0), 01,03,05 e {2,2i }- 

Then cut along the special simplex to produce a 
quadrilateral. Disregarding cases already treated 
by the interval calculations 18.9, we can deform it 
to 

(01, 2, a 2 ,2\/2, 04,2,05,2*0), a { E {2,2i }, 

with diagonals at least 3.2. The result now follows 
from the interval calculations 18.10. 

In summary of (1), (2), (3), we find that by dis- 
regarding cases already considered, we may deform 
the cluster into the form 

(01, 2, o 2 , 2, a 3 , 2, a 4 , 2, a 5 , 2t ), a* e {2, 2t }, 

|«i - v 3 | > 2V2, |vs - v 6 \ > 2^2, \v 2 - v 4 \ > 3.2. 

(4) Assume \v\ — u 3 |,|w 3 — v 5 \ < 3.2. If 
max(ai, a 3 , 05) = 2*o, we invoke interval calcula- 
tions 18.8 and [6, A13] to prove the inequalities. 
So we may assume ai = 03 = 05 = 2. The result 
now follows from interval calculations 18.11. This 
completes the case \v\ — v$\, \v 3 — v$\ < 3.2. 

(5) Assume |t>i — v 3 \, \v 3 — v$\ > 3.2. We deform to 
(01, 2, 02, 2, 03, 2, 04, 2, 05, 2i ), o» € {2, 2t }- 



If 02 = 2*0 and ai = 03 = 2, then the simplex does 
not exist by Section [6, 5.6]. Similarly, a 4 = 2*o, 
«5 = 03 = 2 does not exist. The tec bound gives 
the result except when a 2 = a 4 = 2. The condition 
\v2 — v 4\ > 3.2 forces a 3 = 2. These remaining cases 
are treated with the interval calculations 18.12. 

(6) Assume \vi — Vs\ < 3.2 and \v 3 — v 5 \ > 3.2. 
This case follows from deformations, interval calcu- 
lations 18.8, and 18.13. This completes the proof of 
Inequalities 2. □ 

Proof. (Inequality 1). Label the corners of the 
hexagon V\ , . . . ,v§. The proof to this inequality 
is similar to the other cases. We deform the cluster 
by the method of IV until it breaks into pieces that 
are small enough to be estimated by interval cal- 
culations. If a diagonal between opposite corners 
has length at most 3.2, then the hexagon breaks 
into two quadrilaterals and the result follows from 
Inequality 18.2. 

If a flat quarter is formed during the course of de- 
formation, then the result follows from Inequality 
2 and [6, A13]. Deform until the hexagon has the 
form 

(ai, 2, o 2 , 2, . . . , 06, 2), a t e {2, 2* }- 

We may also assume that the hexagon is convex (see 
[6, Section 4.11]). 

If there arc no special simplices, we consider the 
tec-bound. The tec-bound implies Inequality 1, ex- 
cept when a, = 2, for all i. But if this occurs, 
the perimeter of the convex spherical polygon is 
6 arc(2, 2,2) = 2-k. Thus, there is a pair of antipodal 
points on the hexagon. The hexagon degenerates to 
a lune with vertices at the antipodal points. This 
means that some of the angles of the hexagon are 
7r. One of the tecs has the form C(2, 1.6, n), in the 
notation of [6, Lemma 4.11]. With this extra bit of 
information, the tec bound implies Inequality 1. 

If there is one special simplex, say \vs — Vi\ <G 
[2-^/2,3.2], we remove it. The score of the special 
simplex is [6, Inequalities Ai 3 ] 

vor < 0, t > 0.05925, if max(|ui|, |u 5 |) = 2t , 
vor < 0.0461, r > 0, if M = |u 5 | = 2, 

The resulting pentagon can be deformed. If by de- 
formations, we obtain \v2 — fs| = 3.2 or \vi — v 4 \ — 
3.2, the result follows from Inequalities 18.2 and the 
following two interval calculations 

vor < -0.212 - 0.0461 + 0.137, if y 4 = 3.2, y 5 e 
[2\/5,3.2]. (725257062f) 
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t > 0.54525-0-0.31. if y 4 = 3.2, y 5 e [2\/2,3.2]. 

(977272202f) 

If |ug — ui| = 2\/2, we use Inequality 2 and [6, Ai 3 ] 
unless \v\\ — |t>5| = 2. If |t>i| = \v 5 \ = 2, we use 
the interval calculations 18.3. If a second special 
simplex forms during the deformations, the result 
follows from the interval calculations 18.4. 

The final case of Inequality 1 to consider is that of 
two special simplices. We divide this into two cases, 
(a) The central vertices of the specials are v 2 and 
vq. (b) The central vertices are opposite v\ and v A . 
In case (a), the result follows by deformations and 
interval calculations 18.5. In case (b), the result fol- 
lows by deformations and interval calculations 18.6. 
This completes the proof of Inequalities 1 and the 
proof of the Proposition. □ 

18.3 Groupt 

Let R be a pentagonal region with parameters 

(2, 2, 02, 2, 03, 2, 04, 2, 2, 2y/2), a 2 , o 3) a A e {2, 2t }. 

Since this is a pentagonal region, we may discard 
any edge-length combinations that produce A < 0. 
Assume \v\— v$\ > 3.2, |v 3 — v§\ > 3.2. Under these 
conditions the following inequalities hold. 

vor +0.0461 < -0.212 (583626763) 

t > 0.54525 (390951718) 

18.4 Groupf 

Let R be a pentagonal region with parameters 

(01, 2, a 2 , 2, 03, 2, a 4 , 2, a 5 , 65), Oj € {2, 2t }- 

Since this is a pentagonal region, we may discard 
any edge-length combinations that produce A < 0. 
Assume |t>i— w 3 |, |v 3 — vs\, b 5 = \v 5 —vi\ E [2\/2, 3.2]. 
Under these conditions the following inequalities 
hold. 

vor +0.0461 < -0.212 (621852152) 
T > 0.54525 (207203174) 

18.5 Groupf 

Let R be a hexagonal region with parameters 

(oi, 2, a 2 , 2, a 3 , 2, 04, 2, a 5 , 2, a 6 , 2), Oj G {2, 2t }. 

Since this is a hexagonal region, we may discard 
any edge-length combinations that produce A < 0. 



Assume \vi — vs\,\vi — vs\ £ [2\/2, 3.2], \n — v A \ > 
3.2, |«3 — v 5 \ > 2\f2. Under these conditions the 
following inequalities hold. 

vor < -0.212. (368258024) 
t > 0.54525. (564618342) 

18.6 Groupf 

Let R be a hexagonal region with parameters 

(ai, 2, a 2 , 2, a 3 , 2, a 4 , 2, a 5 , 2, a 6 , 2), a 4 e {2, 2f }- 

Since this is a hexagonal region, we may discard 
any edge-length combinations that produce A < 0. 
Assume \v 2 - v 6 \, \v 3 - v 5 \ e [2\/2,3.2], |u 2 - v 5 | € 
[3.2,3.78], |«3 — «6| > 3.2. Under these conditions 
the following inequalities hold. 

vor < -0.212. (498774382) 
T > 0.54525. (544865225) 

18.7 Groupf 

vor < -0.221 - 2(0.009), if y 4 = 2y/2, y 5 € 
[2y/2, 3.2] , and y 6 e [2t , 2\/2] ■ (234734606f ) 

t > 0.486 - 2(0.05925), if y A = 2V2, y 5 e 
[2\/2,3.2], ye e [2i ,2\/2]. (791682321f) 

18.8 Groupf 

Let Q be a quadrilateral region with parameters 

(ai, &i, a 2 , 2, a 3 , 2, a 4 , 2t ), a 2 , a 3 , a 4 e {2, 2t }. 

Since this is a quadrilateral region, we may discard 
any edge-length combinations that produce A < 0. 
Assume \v± — v 3 \ > 3.2, \v 2 — w 4 | > 2^/2. bi = 
\vi — v 2 \ € [2\/2, 3.2]. Under these conditions the 
following inequalities hold. 

vor (Q) < -0.221 - 0.009. (995351614) 
t (Q) > 0.486 - 0.05925. (321843503) 

18.9 Groupf 

vor < -0.19 - (2/5 - V2)0.14, if y A = 2, y 5 e 
[2^,3.2], and y 6 e [3.2,3.47]. (354217730) 

T > 0.281, if i/4 = 2, y 5 e [2v^,3.2], y 6 G 
[3.2,3.23]. (595674181) 

vor < -0.11, if y 4 = 2, y 5 = 2t , y 6 = 3.2. 

(547486831f) 
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t > 0.205, if y 4 = 2, y 5 = 2* , Ve = 3.2. 18.13 Groupf 

(683897354f) 

vor < 0.009 + (y 5 272)0.14, if y 5 e [2^2,3.2], Let R be a P cnta § onal rc g ion with Parameters 

y4 = y6 = 2 ' (938003786) ( 2 , 2 , o 2 , 2, 2, 2, o 4 , 2, o 5 , 2* ), a 2 , o 4 , a 5 € {2, 2* }. 

Since this is a pentagonal region, we may discard 

18.10 Groupf an y edge-length combinations that produce A < 

0. Assume \vi - v 3 \ e [2^2,3.2], |«i - v 4 \ > 3.2, 
Let Q be a quadrilateral region with parameters | U3 _ „ 5 1 > 3 .2. Under these conditions the following 

~ /x „ ~ r„ „ i inequalities hold, 

(oi, 2, a 2 , 2V2, a 3 , 2, a 4 , 2* ), a, e {2, 2* }. 

vor (i?) < -0.221. (821730621) 

Since this is a quadrilateral region, we may discard to(R) > 486 (890642961) 
any edge-length combinations that produce A < 0. 
Assume \v\— v 3 \ > 3.2, |v 2 — v 4 \ > 3.2. Under these 

conditions the following inequalities hold. 18.14 Groupf 

vor (<2) < -0.221- 0.0461. (109046923) x nifi(J nmn ^ 

' vor < -0.168 - 0.009, y 4 = 2V2, y 5 ,y e e 

t (Q)> 0.486. (642590101) [2* ,2\/2]. (341667126+) 

T > 0.352- 0.05925, y 4 = 2y/2, y 5 ,y 6 e [2* ,2\/2]. 

18.11 Groupf (535906363 t ) 

Let R be a pentagonal region with parameters -^g ^5 Group 

(2,2,a2,2,2,2,a 4 ,2,2,2to), a 2 ,a 4 e {2,2* }. voro < -0.146, if y 5 , y 6 G [2y/2, 3.2]. (516537931) 



Since this is a pentagonal region, we may discard 



^(^(j/i, . . . ,y e )) + T (S(y 1 ,y' 2 ,y 3 ,2,y 5 ,2)) > 0.31, 



any edge-length combinations that produce A < 0. 11 2/5,2/6 e [2^2, 3.2], and y' 2 G {2,2* }- 
Assume |wi-w 3 | G [2\/2, 3.2], |i»3 — G [2V2,3.2]. (130008809) 



Under these conditions the following inequalities 
hold. 



18.16 Group 



vor (i?) < -0.221. (160800042) 

MR) > 0-486. (690272881) VOr ° < -°' 084 ' if 2/5 G [ 2i o,2V2], y 6 G [272,3.2]. 

y ' (531861442) 

vor < -0.084 - {y 5 - 272)0.1, if yi = y 3 = y 4 = 2, 
18.12 Groupf y 6 = 2t , y 5 G [2-3,3.2]. (292827481) 

t J. d v, + , • ♦ vor o < 0.009 + (j/5 - 2v/2)0.1, if j/5 G [2V2, 3.2], 

Let i? be a pentagonal region with parameters = = = = 2 (710875528) 

(ai,2,2,2,2,2,2,2,a 5 ,2t ), 01, o B € {2,2* }. t > 0.176, if y 5 G [2* ,2V2], y 6 e [272,3.2]. 

(286122364) 

Since this is a pentagonal region, we may discard 
any edge-length combinations that produce A < 0. 
Assume \v\ — v 3 \ > 3.2, \v 3 — v$\ > 3.2. Assume 

dih(0, v 3 , v 4 , v 5 ) + dih(0, v 3 , v 5 , vi) + dih(0, v 3 , v 1 ,v 2 ) 
> dih(5(2,2,2,3.2,2,2)) 
= arccos(— 53/75). 

Under these conditions the following inequalities 
hold. 

vor (i?) < -0.221. (713930036) 
t (R) > 0.486. (724922588) 
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